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1 - SUMMARY 



The impact of acid deposition on the aquatic ecosystems ot 
the Canadian shield region of central Ontario has, in recent 
years, become the subject of increased study. An important 
component in this study is the evaluation of acidity trans- 
mission by surface water. 

To this end, a model was required to estimate local runoff 
and to simulate the movement of runoff through a stream 
network. Since stream networks in the region of interest are 
dominated by interconnected lakes, these networks are repre- 
sented in the model as interconnected chains-of-lakes. The 
model is divided into three parts. The first part is used to 
calculate subbasin runoff, the second part is used to route 
the subbasin runoffs through the system and the third 
optional part can be used to calibrate the model for a 
specific basin. 

The runoff prediction utilizes equations to estimate average 
monthly runoff for any basin of interest. The equations were 
derived from a multiple linear regression analysis and are 
applicable only to basins similar to those used in deriving 
the equations. This limits the use of the equation to basins 
up to about 750 km2 in area located in the Canadian Shield 
region of central Ontario (from Lake Simcoe to Sudbury and 
east from Georgian Bay to Pembroke on the Ottawa River) . The 
variables incorporated in the runoff prediction equations are 
either physiographically or climatologically based. 

Two separate sets of runoff prediction equations were devel- 
oped. One set of equations is based on correlations of 
historic sequential aata while the other set is based on 
correlations of ranked data. 
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The relationships which were derived using historic sequen- 
tial data are applicable to 

- replicating or estimating runoff in a particular historic 
month in an ungauged catchment 

- reproducing the statistical behavior of monthly runoff 
(i.e., mean, standard deviation, skew of the monthly 
distribution) . These estimates will also replicate data 
for a specific historic month. 

The relationships which were developed using ranked data are 
applicable to 

- reproducing the statistical behavior of monthly runoff 
(i.e., mean, standard deviation, skew of the monthly 
distribution) . These estimates will not replicate data for 
a specific historic month. 

The second part of the model, the routing model, is determin- 
istic. This part routes the individual subbasin runoffs 
through a stream network. An estimate of discharge is 
provided at selected points, resulting from both local runoff 
and from inflows originating in upstream basins. In addition 
to accumulating contributions from upstream basins, the 
routing model accounts for natural regulation from major 
unregulated lakes and controlled regulation from the opera- 
tion of river or lake structures. 

The third portion of the model may optionally be used for 
calibrating the model to specific drainage basins where 
historic discharge data for the basin is available. The two 
calibration techniques which can be optionally selected 
within this part of the model are 



- historic replication of a specific basin discharge record 

- improvement of the runoff regression equations for a 
specific data record. 

The historic replication technique can be used tor determin- 
ing the distribution of runoff within a basin which has a 
long historic record of measured discharges. The technique 
cannot be used for extending data beyond the perioa of 
measurement at a downstream gauging station. This replica- 
tion adjusts the basin time series response to perfectly 
match a single downstream gauge. The adjustment is then 
distributed over the subbasin on the basis of subbasin 
runoff, as computed by the regression equations, as opposed 
to using the common technique of prorating based on subbasin 
area alone. Although the selected approach is physically 
realistic, it is not possible to judge its accuracy using the 
available data base. 

The runoff regression equation improvement technique can be 
used for extending data to periods when discharge records are 
not available but climatic data is available. 

The South Branch Muskoka River was selected for testing the 
chain-of-lakes model. Rule curves and historic lake level 
data were used for operating and verifying the model. The 
results of three test runs are presented. These runs are 

- historic replication (including a month-by-month calibra- 
tion) 

- statistical simulations using 

- historic sequential runoff prediction equations 
(including calibration to improve the runoff prediction 
equations) 

- ranked runoff prediction equations (including calibration 
to improve the runoff prediction equations). 



The results shows that considerable adjustment is required to 
calibrate the model on a month-by-month basis. 

Similarly, calibration using an improvement to the runoff 
regression equations resulted in significant changes to the 
relationships. 

A comparison of the two statistical simulations shows that 
the ranked approach is better suited for data extensions 
where the statistical behavior of runoff is of more concern 
than particular historic runoffs. 

From the results of the test cases on the South Branch 
Muskoka River, and from the natural runoff prediction scatter 
diagrams, it can be concluded that flow estimates for 
ungauged uncalibrated catchments are only reliable within a 
factor of approximately two. Flow estimates using the 
calibration routine of the model will be very much more 
reliable, although the accuracy of flow estimates for 
ungauged subbasins in a gauged basin can not be evaluated 
using the current data base. It is therefore recommended 
that the model be tested on a watershed for which subbasin 
discharge data is available. 

There are several appendixes to this report which give 
additional details of the study as follows: 

- discharge calculations for the routing model 

- details of the calibration procedures 

- user's guide to the input data 

- example input and output for a historic replication run for 
the South Branch Muskoka River. 
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In recent years, the phenomenon of acid deposition has become 
a serious environmental problem. The most widely publicized 
result of acid rain is the destruction of many forms of life 
within any given lake. This destruction of life occurs once 
the pH level has been reduced to a critical level through 
accumulation of acid precipitation. 

In order to assist in the monitoring and evaluation of the 
effects of acid deposition on the ecosystem, there is a need 
to predict changes in pH levels of lakes tor a given level of 
acid deposition. As part of the response to this need, the 
object of this study is the development of a model for 
predicting or simulating flow at any point within a stream 
network. Ultimately, the model is to be used to calculate 
the hydrologic input for a water quality model used to 
simulate water quality effects of acid rain deposition. 

Acid deposition exists throughout Canada, however, the 
effects are most acute in regions which do not have a natural 
buffer, such as limestone. Of particular concern is the 
Canadian Shield region of central Ontario. The thin soil and 
bedrock offer little buffering. Because of its significance, 
this region has been selected for detailed study. The region 
is characterized by thousands of lakes with a wide range of 
sizes and configurations. On a monthly basis, the flow 
regime of a stream network in this region is usually con- 
trolled by precipitation and lake storage, with in-channel 
storage having little effect. Consequently, a drainage basin 
or stream network may be modeled as a chain of lakes. 
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The development of a chain-of-lakes model is potentially very 
complex. However, the model presented herein is simple in 
both its theory and in its application. 

The model can be considered in three parts. The first part 
is statistically based and involves the development of runoff 
prediction equations from physiographic and climatologic 
data. The second part of the model is physically based and 
is primarily a routing model with the focus on routing runoff 
from individual basins through a stream network, or a network 
of basins. The third portion of the model may optionally be 
used for calibrating the model to a specific drainage basin. 

In this report, the development of the runoff prediction 
equations is discussed first. Subsequently, the flow predic- 
tion model is described. It includes the results of the 
runoff prediction equations and routing and calibration 
techniques. Finally, application of the model to the South 

Branch Muskoka River is discussed. The south Muskoka drain- 
age basin was selected for verification purposes because of 
the availability of a comprehensive data base. Appendix E 
provides a user's guide to the input data which is required 
to run the model. Appendix F is a test run for the South 
Branch Muskoka River, including input data files and an out- 
put listing. Other appendixes list data used for the runoff 
prediction correlations and provide greater detail of the 
routing model than appears in the main text. 
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3 - NATURAL RUNOFF 
PREDICTIONS 



Multiple regression analyses were used to develop monthly 
natural runoff prediction equations for the Canadian Shield 
region of central Ontario. The relationships are based on 
readily available climatologic and physiographic data. 
Special care was taken to ensure that the resultant relation- 
ships were physically meaningful. However, the reliability 
of the prediction equations is influenced by the meagre data 
base which is available. 

In any application of statistics to hydrology, particularly 
for analyses using limited data, there are certain dangers, 
misuses and misconceptions which must be recognized. In this 
study, these possible problems are especially relevant and 
should be mentioned. Yevjevich (Reference 1) has stated them 
as follows. 

(a) Expectations are often greater than the available data 
and the best results which statistical methods can 
produce. 

(b) False expectations that statistical methods may in some 
way replace the scarcity or unreliability of data, which 
is impossible. 

(c) The tendency to replace the necessary judgment and 
interpretation of physical facts by statistical computa- 
tions or statements. 

Within this section, techniques for developing a long period 
of reference hydrology are discussed. This long period 
reference hydrology is then available for analyzing the 
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long-term effects of acid precipitation fallout. Data 
extension is frequently required for developing long period 
reference hydrology and possible extension techniques are 
presented along with the two techniques which were selected 
for detailed investigation. 

The remainder of this section describes these extension tech- 
niques and the development of the statistical runoff predic- 
tion equations. A discussion of each of the following is 
included in this section: 

- techniques for data extension 

- data availability 

- methods of analysis 

- results. 



3.1 - Techniques for 
Data Extension 



There are various techniques that can be used for extending 
streamflow records to obtain a larger synthetic record for 
use as reference hydrology. This long period reference 
hydrology might be of interest for analyzing the long-term 
effects of acid precipitation fallout. Some of the available 
techniques are listed and discussed briefly as follows. 



(a) Runoff/Runoff 
Correlations 



Correlations of this type can be made between short- and 
long-term gauging stations in the same river basin and/ 
or in nearby basins. This technique is applicable to 
many circumstances, such as the original Serpent River 
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model, but could not be easily used in a general manner 
for estimates on a totally ungauged catchment as 
required in the current study. 



(b) Rainfall/Runoff 
Correlations 



Correlations of this stype can be made between long-term 
rainfall stations and short-term streamflow gauging 
stations, provided both types of stations are located in 
a homogenous climatologic and physiographic region. 
Unfortunately, in many circumstances rainfall stations 
are not well located for obtaining consistent estimates 
of basinwide rainfall. Consequently, the correlations 
are weak and the standard error of estimate of runoff 
from rainfall observations is large. 

The principle analyses for the current study use 
rainfall/runoff correlations extensively. However, in 
the current case the errors associated with the weak 
correlations are partially corrected by calibration 
against downstream flow data. This calibration is 
performed for an ungauged subcatchinent by adjusting the 
general runoff relationships, so that the total dis- 
charge from the catchment matches the records for the 
downstream gauge, A detailed description of this 
calibration technique is presented in Section 4.3.2. 

( c) Markov Modeling 

Synthetic streamflow generation can be done using a 
Markov type model. Such models view streamflow as a 
time series of serially dependent values, where each 
value has a deterministic component and a random 
component. 



The deterministic component can be comprised of a trans- 
fer function to convert precipitation to runoff. To 
apply Markov type modeling to the current study would 
require an appropriate transfer function. The limita- 
tions imposed by the current data set prevent applica- 
tion of this highly specialized hydrologic tool. 



(d) Ranked Rainfall/Runoff 
Correlation 



Ranked rainfall/runoff correlation is a correlation of 
rainfall in a basin with runoff, on the basis of 
frequency of occurrence. This method attempts to 
compensate for inadequate disposition of rainfall 
stations in or near a basin by bypassing the determin- 
istic relationship between rainfall and runoff to 
directly consider the probabilistic nature of the hydro- 
logic cycle. The ranked correlation approach assumes 
that both the basin and the meteorologic station are 
located in the same climatic zone and that their long- 
term precipitation density functions are similar. The 
ranked correlation approach involves reproducing the 
statistical behavior of a basin without regard for 
reproducing or predicting specific historic runoff. The 
method is therefore most applicable when analyzing 
long-term trends and overall basin behavior. 

The two techniques which have been adopted in this study are 

- rainfall/runoff correlation 

- ranked rainfall/runoff correlation. 

Both methods can be used for extending the discharge record 
beyond the period of available discharge data. However, if 
an estimate of runoff for a particular historic month is 
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required, then only the historic sequential rainfall/runof t 
correlation is applicable. For general extensions of data, 
either technique can be used. However, Section 5 (Test 
Results) demonstrates that the statistical behavior of the 
basin is reproduced more closely using the ranked rainfall/ 
runoff correlation technique. 



3.2 - Data Availability 

In total, 13 drainage basins ranging in size from 22 to 
741 kin2 were selected for analysis. These drainage basins 
were selected based on the availability of natural discharge 
records. The selected basins are illustrated on the location 
map in the pocket at the back of this report. 

The rest of this section presents information concerning the 
discharge data, the physiographic characteristics, and the 
climatologic data. 

3.2.1 - Streamf low 

Monthly discharge data was obtained from Water Survey 
of Canada (WSC) publications (References 2, 3 and 4) 
and is tabulated in Appendix A. The period 1970 to 
1981 was selected for analysis since few gauges 
predated 1970. The available period of record for each 
of the 13 gauges is illustrated on the bar chart in 
Figure 1. 



3.2.2 - Physiographic 

Characteristics 



All physiographic data were obtained from 1:50 000 
scale topographic maps published by the Surveys and 
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Mapping Branch of the Department of Energy, Mines and 
Resources. The data were obtained by measurements 
using a digitizer linked directly to Acres VAX com- 
puter. The map sheet reference numbers, which are 
applicable to each of the 13 selected basins, are 
tabulated in Table 1 along with other pertinent physio- 
graphic data. 

The particular physiographic data which were measured 
are discussed below. 



(a) Total Basin Area - A„ 



The total basin areas tributary to each stream 
gauge are published by WSC. These areas were 
checked to ensure that the drainage basins which 
were drawn on the topographic maps corresponded to 
the basins as defined by WSC. 



(b) Area of Lakes - A 



The area of lakes in each of three size categories 
was measured for each basin. These sizes are 

- less than 0.01 km2 

- from 0.01 to 0.5 km2 

- greater than 0.5 km2. 

As discussed further on in this report, the avail- 
able data base was not suitable for regression on 
a large number of independent variables. Conse- 
quently, these three categories were later lumped 
together into a total area of lakes for each 
basin. 
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TABLE 1 



SUMMARY OF PHYSIOGRAPHIC 


DATA 










USED IN THE REGRESSION ANALYSIS 






















1:50000 


River 


Gauge 
Number 


Area of 
Basin 


Area of 
Swamps 


Area of 
Lakes 


Climatological 
Station 


Topographic 
Map Index 


N. Magnetawan 


2EA005 


172.0 


10.1 


16.5 


Bark's Falls 


31E/11, 14 


N. Magnetawan 


2EA010 


149.0 


3.8 


4.7 


Burk's Falls 


31E/11, 14 


Shawanaga River 


2EA012 


235.0 


37.6 


24.2 


Burk's Falls 


41H/8, 9 
31E/12 


Whitson River 


2CF007 


117.0 


0.4 


0.0 


Sudbury 


411/10, HE 


Whitson River 


2CF008 


155.5 


29.1 


8.7 


Sudbury 


411/10, HE 


Vermilion River 


2CF100 


635.0 


45.1 


61.4 


Sudbury 


41I/10W, HE, 
14E, 15W 
41P/3E 


Duchesnay River 


2DD00 8 


90.0 


2.6 


0.9 


North Bay 


31L/5, 6W 


Veuve River 


2DD012 


741.0 


66.5 


23.4 


North Bay 


411/7, 8, 9, 10 


La Vase River 


2DD013 


70.4 


3.6 


0.2 


North Bay 


31L/3, 6 


Commanda Creek 


2DD015 


106.0 


10.8 


5.6 


North Bay 


31E/13, 14 


Rosseau River 


2EB103 


130.0 


10.1 


4.9 


Huntsville 


31E/4, 5, 6 


Bob Creek 


2HF004 


22.0 


2.9 


0.5 


Minden 


31D/15 
31E/2 


Indian River 


2KC014 


443.0 


21.0 


11.0 


Chalk River 


31F/H, 12E, 
13E, 14 
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(c) Area of Swamps - A 



The area of swamps was measured for each basin. 
The definition of swampland is dependent on topo- 
graphic mapping detail because the swamp area was 
measured from areas depicted symbolically on the 
maps as swampland. The area of swamps normally 
varies throughout the year and from year to year. 

3.2.3 - Climatologic Data 

Cliraatologic stations with sufficient record, that are 
within or near a basin, were selected from the 
Cliroatoiogical Station Data Catalogue for Ontario 
(Reference 5). The selected stations are 

- Sudbury 

- North Bay 

- Burk's Falls 

- Huntsville 

- Minden 

- Chalk River. 

The particular station which applies to each of the 
drainage basins is listed in Table 1. The stations are 
also shown on the location map in the pocket of this 
report. Figure 1 shows the available period of record 
for each station. Appendix B is a compilation of the 
pertinent data for each of the 6 stations. This data 
was compiled from Monthly Records of Meteorological 
Observations in Eastern Canada by the Atmospheric 
Environment Service (AES) (Reference 6). 

The particular climatologic data which was tabulated is 
discussed below. 
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{ a) Precipitation - P 

The total monthly precipitation for each month was 
tabulated. 



(b) Mean Monthly Temperature - 
Tmax or Tmin 



Three individual measures of the mean temperature 
for each month are available based on 

- raaxirauni daily temperature 

- minimum daily temperature 

- mean daily temperature. 

The mean daily temperature which is tabulated by 
AES is just the average of maximum and minimum 
daily temperatures. Consequently, only the mean 
monthly temperatures based on maximum and minimum 
daily temperature were extracted from the records 
for the purpose of this study, 

(c) Degree Days Above 

Freezing - Dmax or Dmin 

The number of degree days above freezing (CC) for 
each month was derived from AES published records. 
Two measures of this parameter were used — based on 
either the maximum daily or minimum daily 
temperature. 

In the summer, this variable is equivalent to the 
average temperature and consequently was not tabu- 
lated. However, in the winter and spring, this 
variable reflects the importance of the number of 
melt degrees rather than the overall temperature. 
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3.3 - Method of Analysis 

The objective of this analysis was to develop monthly natural 
runoff prediction equations for ungauged catchments in the 
Canadian Shield region of central Ontario. The previously 
discussed independent variables (Ag, a,, p^ p^ Tmax, 
Tmin, Etaax, Dmin) were used in a multivariate regression 
analysis to derive these prediction equations. A stepwise 
multiple linear regression computer program was used for the 
analysis. The general characteristics of this program and 
the actual analysis which was undertaken are discussed in the 
following two subsections. 

3.3.1 - Regression Analysis 
Computer Program 

The program performs the stepwise regression analysis 
using a simple correlation coefficients matrix. Any 
combination of one dependent variable and one or more 
independent variables may be selected for analysis, 
using the input either directly or transformed. 

In the stepwise procedure, independent variables are 
added one at a time to the regression, giving the 
following intermediate equations 



Y = B' + b;x. tU 
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Y = B^ + B^X^ + B'^X^ 12} 



Y = B^' + B^'Xj^ + ^y^2 ^ ^3*^3 ^^^ 



In this manner, it is possible to obtain valuable sta- 
tistical information on the improvement in predictive 
capability at each step of the calculation. 
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At each step, those independent variables not included 
in the regression are inspected to find the one that 
will give the greatest reduction in the variation of Y. 
This variable is then tested for significance; i.e., 
the computed F-ratio of the variable is compared to the 
supplied critical F-value. If the computed value is 
greater than the critical value, the variable is con- 
sidered significant and is added to the regression 
solution. Where appropriate, "forced independent vari- 
ables" are defined to always be included in the 
equations and are not tested tor significance. 

After each new variable has been added, those variables 
already in the regression are inspected to see if any 
of them can now ba deleted because their contribution 
to the reduction in the variation of Y is no longer 
significant. Those variables which have a computed 
F-ratio less than the critical value are considered 

insignificant and are deleted from the regression 
solution. 

This process is continued until no more variables are 
added or deleted. Thus, the final regression solution 
contains only those variables that are statistically 
significant (except possibly forced variables). The 
general equation which results is thus 



Y = B^ + B.X, + B_X. + ...B„X„ (4) 

U 1 i 2 2 n n 



Under user control, it is possible to force one or more 
of the independent variables into the regression solu- 
tion. These variables are added before any other 
variables, with no regard for statistical significance. 
The remaining independent variables are then added to 



I 
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the equation in the normal stepwise procedure as out- 
lined in the preceding paragraphs. Once a variable has 
been forced into the regression, it will be kept 
regardless of the statistical significance. 

The application of transformed variables can alter the 
linear form of Equations 1 through 4 and, the multiple 
linear program can then be used for general multivari- 
ate regression analysis. Two sets of transformations 
were investigated for application in the current study. 
These transformations are the following. 



(a) Log Transform of 
the Dependent and 
Independent Variables 



Substituting log Y, log X;^ , log X2 , etc, for 

Y, Xi , X2 f etc, in Equations 1 through 4 
results in the following general equation 



log y = B„ + B.log X, + B„log X- + ,_, 

11 22 (5) 

. . .B„log X 
n ^ n 



This equation simplifies to the following Power 
Law Relation 



Bp. B, B- B 

Y - 10 ^ X, X- '^...X " (6) 

1 2 n 



(b) Log Transform of 

the Dependent Variable 

Substituting log Y for Y in Equations 1 through 4 
results in the following general equation 



log Y = B„ + B,X, + B-X^ + ...+B X (7 

u 1 i 2 2 n n 
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This equation reduces to the Exponential Relation 
which follows 



A 1 O 

Y ^ 10 " b. ^ b^ '^.. .b " (8) 

12 n 



where 



bj^ = 10 ■^, b^ = 10 '', etc 



3.3.2 - Alternative Analytic Approaches 

There are two alternative approaches which were inves- 
tigated for the regression analyses. These two 
approaches can be simply described as a one-step 
analysis and a two-step analysis. Each technique is 
described briefly below, followed by a general discus- 
sion of the relative merits and limitations of each. 

(a) One-Step Analysis 

Analyses using the one-step technique involve 
determining monthly relationships as follows. 

R = Function{ Physiographic Data, 

Climatologic Data) (9) 

where R is the monthly average discharge for a 
natural catchment. 

Both dependent and independent variables are 
lumped together into one long record containing 
all of the years of record for all of the basins 
for any particular month of the year. This tech- 
nique is commonly used in hydrologic regressions 
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and is called the station-year method by Yevjevich 
(Reference 1). It must be recognized that the 
effective number of observations is much smaller 
than the total number of observations, physio- 
graphic data (A3, Al, Ag) for the indepen- 
dent observations are repeated in blocks for each 
basin and are not variable throughout the entire 
set of observations. Similarly, the same meteoro- 
logic station is used for several basins and the 
precipitation and temperature variables are also 
repeated in blocks. Because of the obvious mutual 
correlation of these variables, it is readily 
apparent that the total sample size can only 
provide as much information as some lesser number 
of totally independent observations. In this 
study the total data set contains approximately 
100 observations for each month. However, the 
equivalent number of independent sets of observa- 
tions would be much less. 

( b) Two-Step Analysis 

Analyses using the two-step technique involve 
separately determining the following monthly 
relationship for each basin 

R = Function( Climatologic Data) (10) 

This analysis results in a set of coefficients 
which describe the runoff in each particular basin 

(i.e., Bq , B^ , B2, etc, from Equation 4 for 
each basin) . The second half of the two-step 
analysis involves explaining the variation in 
these coefficients as follows 
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B„, B, ...B = Functions{ Physiographic Data) (11) 



The derivation of natural runoff for an ungauged 
catchment therefore requires two steps as follows. 

- Determine the coefficients Bq , Bi...Bn 

based on the basin-specific physiographic data 
(from Equation 11). 

- Determine the basin runoff based on the coeffi- 
cients and climatologic data (from Equation 10). 

The two-step technique has the following advantages. 

- The data limitations are readily apparent at all 
times. It is always immediately apparent that each 
regression of the form of Equation 10 is based on a 
relatively short record (4 to 12 years in the current 
study) , 

- There is no compromise of formal statistical rules 
concerning dependency and independency. 

The one-step technique has the following advantage, 

- The total number of regression analyses is greatly 
reduced. In the current study, 15 or 16 runs would 
be required for the two-step approach as opposed to 1 
run for the one-step approach. When compounded by 
the requirement for 12 monthly regressions, it is 
apparent that the one-step technique is much simpler 
and more straightforward. 



I 
I 
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A test was performed using both techniques for the 
month of May. An evaluation of the final results was 
done using comparisons of observed and estimated 
discharges. It was discovered that the final estimates 
using both techniques were quite similar. The standard 
error of estimate was marginally lower for the one-step 
approach on 10 of the 13 basins and marginally lower 
for the two-step approach on the other 3. No distinct 
pattern was observed for this variation between the 
techniques. Because of the computational and data 
handling advantages of the one-step technique, it was 
adopted for all of the remaining regression analyses in 
this study. 



3.4 - Results of 

Regression Analysis 

The analyses were done using the multiple linear regression 
computer program described in Section 3.3.1. 

The independent variables which were used in the initial 
regressions are the following (where i refers to the current 
month) . 

- Monthly precipitation - P. , P. , , P. ^, P. -, , P- . (m) 

1 1—1 1 — *■ 1— J 1 — 4 

- Mean monthly temperature based on the maximum daily 
temperature Tmaxj^, Tmaxi_i 

- Mean monthly temperature based on the minimum daily 
temperature Tmin^^, Tmini_i 

- Degree days above 0°C based on the daily maximum tempera- 
ture - Dmax^, Dmaxi_i, Dmaxi_2, Dmaxi_3 
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The multiple linear regression computer program described 
I in Section 3.3.1 was used for general raultivariant regression 

analyses using the following transformations. 

I 
I 
I 
I 
I 

It was found that satisfactory results could only be obtained 
H for a linear relation. The power law relation created 
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- Degree days above O^C based on the daily minimum tempera- 
ture - DmiHj^^ Draini_i, Dmini_2r Dmini_3 

- Area of basin - Ag (km2) 

- Area of lakes - A^ (km2) 

- Area of swamps - Ag (kra2) 

The variables Dmax and Dmin were only included in the analy- 
sis for the period from November to April when the impact of 
melt degrees could be significant. 

The dependent variable was the current monthly recorded 
natural discharge, Rj^. 



- Direct independent and dependent variables (linear 
relation) . 

- Logarithm of both independent and dependent variables 
(power law relation). 

- Logarithm of dependent variable and direct independent 
variables (exponential relation) . 



unacceptable anomalies if one or more of the independent 
variables was very small. The exponential relation did not 
have this deficiency, however, the results were not as good 
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as for the simple linear relation. Consequently, all of the 
final regression equations are multiple linear equations. 

Many alternative combinations of variables were tested before 
the final equations were derived. The F-statistic at the 90 
percent confidence level was used as the criterion of accep- 
tance of independent variables into the predictive equation. 
However, this was tempered with judgment, when considering 
acceptance or rejection of variables, based on the physical 
significance of the variables. 

Precipitation and area of basin are always present in the 
relationships because of the physical relation between the 
runoff and the volume of water (P x Ag) falling on the 
basin. Precipitation was initially included for the current 
month and also lagged by up to 4 months (Pj^ to Pi_4). 
It was found that in some months the precipitation occurring 
several months prior to the current month was statistically 
more significant than the current months precipitation. For 
basins ranging in size from 22 to 741 km2^ this anomally 
seems physically unreasonable. Consequently, the independent 
variable describing precipitation was changed to the average 
precipitation for a period which varied in duration from 1 to 
4 months. There was also a secondary advantage to this 
choice because the number of independent variables was 
reduced, thereby providing more degrees of freedom to the 
final relationship. 

The general regression equations which resulted for both the 
ranked and historic sequential rainfall/runoff correlations 
are as follows. 

April to November 



R = a, + A^(a_P + a.) + A^(a.P + a.) + A.a^ (12) 

1B2 3 24 5 56 
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where 

^1 to ag are constant coefficients which vary from 
month to month 

P is either i?^ + Pi-i}/2 or (Pi + Pi_i + 
Pi-2)/3 or [Pi + (Pi-i + Pi-2 + 
^i-3)/3]/2, depending on the specific month 

December to March 



R = a^ + A^a^ + AgU3 + P(a^ + a^D)] (13) 



where 

D is the monthly degree days above O^C. D is either 
(Droaxi + Dmaxi-i)/2, (Dmini + Dmini-i)/2 or 
Dmini, depending on the specific month. 

The general relationship for April to November (Equation 
12) can be reduced to the following basin-specific runoff 
equation. 



R s ttj^ + B^p (14) 



where 



a = a, + a-A„ + a_A^ + a^A„ 
1 1 3B 5L 6s 



^1 = ^2^B "■ ^A 



Similarly the general relationship for the December to 
March period (Equation 13) can be reduced to the following 
basin-specific runoff equation 



I 
I 
I 
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R ^ a^ + ^2? + YjPD (15) 



where 



°'2 = ai + a2AL + asAs 



^2 = a4AB 



^2 = asAB 

These reduced equations are simple linear rainfall/runoff 
prediction equations. The comparison o£ the one-step and 
two-step regression techniques which was discussed in 
Section 3.3.2 is analogous to a comparison between 
Equations 12 and 14 and Equations 13 and 15. In the two- 
step analysis, the coefficients a^, 0-2, ?>i, 62 and 

"^2 would be determined for each of 13 basins. Prediction 
equations would then be derived for a^ , 0.2, 3i, 60 

and Y2, based on physiographic characteristics for each 
basin. In the one-step analysis, a^^, a2, Bi, 62 

and Y2 are derived implicitly within Equations 12 and 13 
using an overall lumped data set. As can be seen, the 

final prediction equations, using either the one-step or 
the two-step technique, would have similar forms. 



3,4.1 - Historic Sequential 
Correlations 



The regression analyses produced a set of coefficients 
(aj^ to 35) for both general equations (Equations 12 
and 13). These coefficients are tabulated in Table 2. 
The scatter diagrams relating measured and calculated 
runoff for each month are plotted in Figures 2 to 7 . 



TABI£ 2 

OOEFFICIEhrrS FOR HISTORIC SEQUEMTIAL OQRPEIia'IOtiS 

General Bquatican: (^ril to Novenber) 

R = a^ + Ag (a^P + 33) + A^ (a^P + a^) + a^A^ 





Precipitation 

m 


Apr 


(Pi + Pi_i)/2 


May 


(Pi + Pi_i)/2 


Jun 


(Pi ^ Pi_l)/2 


Jul 


<^i ^ Pi-1 ^ Pi-2'/3 


Aug 


<^ ^ "^i-l ^ **i-2)/3 


Sep 


t^ ^ Pi-1 ^ '?i^2>/3 


Oct 


<Pi ' h-1 ^ ^i-2>/3 


ttov 


(Pi ^^ (Pi_i + Pi_2 + 




P._3)/3.0)/2 



h 


^2 


^3 


*t 


*5 


"6 


r 


Siqnif 
(90 %) 


leant 
(95%) 


-0.926 


0.529 


0.020 


0.0 


-0.267 


0.267 


0.83 


0.64 


0.80 


-0.384 


0.051 


0.023 


0.0 


0.213 


-0.108 


0.82 


0.64 


0.80 


0.071 


0.119 


0.0 


0.0 


0.114 


-0.074 


0.80 


0.64 


0.80 


0.399 


0.064 


0.0 


0.927 


0.0 


-0.066 


0.67 


0.64 


0.80 


0.080 


0.051 


0.0 


0.522 


0.0 


-0.038 


0.64 


0.64 


0.80 


0.050 


0.170 


-O.O086 


0.0 


-0.022 


0.0 


0.65 


0.62 


0.70 


0.064 


0.313 


-0.017 


0.0 


-0.024 


0.0 


0.78 


0.62 


0.70 


-0.060 


0.296 


-0.012 


0.0 


0.0 


0.0 


0.66 


0.48 


0.57 



SEE 
(in3/s)* 


Q 3 
(mVs)* 


6.48 


13.8 


3.35 


6.91 


1.33 


2.42 


1.10 


1.52 


0.86 


1.13 


1.58 


1.76 


1.87 


2.94 


2.21 


3.30 



General Equation: (Deoaiter to Mardi) 

R = a^ + a2Aj^ + Ag (a^ + P (a^ + a^D) ) 

Precipitation Degree Days lu 

(m) (OQC) * 

Dec (P. + P. J/2 (Dnax. + Dmax. ,)/2 0.161 

X i— J, 1 Jt— i 

Jan (P^ + Pi_i)/2 (Qnax^ + DinaXj^_j^) /2 0.244 

Feb (P. + P. ,)/2 (ftidn. + Dnln. ,)/2 0.312 

Mar (P. + P._, + P _J/3 Dndn. 1.033 

1 X^X X 4b X 



* r ^ is the coefficient of determinaticffi 

* PFR (m /s is the stanc3ard error estimate 

* Q (mVs) is the mean ironthly discharge for all stations used in the analysis 
**Assiining 10 effective independent c3bservations 



'^ 


*3 


^ 


«5 


r^ 


Slgnif 
(90%) 


icant r2** 
(95%) 


SEE 
(m3/s) 


(mVs) 


0.058 


0.0 


0.036 


0.0013 


0.84 


0.64 


0.80 


1.13 


2.71 


0.048 


0.0 


0.025 


0.00066 


0.76 


0.64 


0.80 


0.59 


1.33 


0.042 


0.0 


0.018 


0.049 


0.82 


0.64 


0.80 


0.90 


1.40 


-0.039 


0.0071 


0.0 


0.043 


0.76 


0.64 


0.80 


2.80 
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FIG. 2 
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SCATTER DIAGRAM OF RUNOFF SEQUENTIAL CORRELATION 
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The coefficients which are listed in Table 2 respond as 
expected and are physically reasonable. In all months, 
one or more of the coefficients is zero. These zero 
coefficients occur when the particular independent 
variable by which the coefficient is multiplied is not 
significant for that particular month. 

Table 2 lists values of the coefficient of determina- 
tion (r2) for each monthly relationship. Significant 
values of r2 at the 90 and 95 percent confidence 
levels are also tabulated. These significant r2 
values are based on the assumption of 10 effective 
independent observations (see discussion of Section 
3.3.2). The number of independent variables included 
in the relations varies from 2 to 4 and, hence, based 
on an assumption of 10 effective independent observa- 
tions, the number of degrees of freedom varies from 5 
to 7. The monthly relations for February, April, May, 
June, October, November and December are significant at 
the 95 percent level. The remaining months are all 
significant at the 90 percent level. The relations for 
July and August are only marginally acceptable at the 
90 percent confidence level. 

The scatter diagrams (Figures 2 to 7 ) indicate consid- 
erable variability between measured and calculated 
discharges. While reasonable relationships are indi- 
cated for most months, a poor predictive capability is 
indicated for July and August. 

The derived equations represent the best possible 
relationships for statistically estimating the basin 
runoff in an ungauged catchment for a particular 
historic month. 



36 



3.4.2 - Ranked Correlations 

In general, it is preferable to develop prediction 
equations for runoff based on the historic correspon- 
dence of runoff and climatic behavior. However, in 
cases where monthly basin precipitation estimates are 
not reliable and the objective is to reproduce or 
simulate the statistical behavior of basin runoff over 
a period of time, a ranked correlation analysis can 
sometimes give superior results. This problem of 
estimating basin precipitation is particularly acute 
where the nearest precipitation station is a long way 
from the basin. In this case, the statistical 
distribution of monthly rainfall at the gauge and over 
the basin may be similar but in any given month the two 
may differ widely. 

To utilize the ranked approach, both precipitation and 
runoff are ranked from maximum to minimum within a 
specific basin. These two ranked vectors, along with 
the physiographic data, are then incorporated into the 
regression analysis. 

During some months, particularly during spring runoff, 
the degree days of melt are more significant than pre- 
cipitation. In such cases, the ranking procedure is 
based on degree days and runoff. When ranking the 
degree days, the precipitation is arranged to histor- 
ically correspond with the ranked degree days. 
Similarly, if precipitation dominates, but degree days 
are still a contributing variable, the degree days are 
arranged to historically correspond with the ranked 
precipitation. 
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The analyses for the ranked approach resulted in the 
same general equations as were derived for the historic 
sequential correlations. The coefficients which were 
derived are tabulated in Table 3. The scatter diagrams 
relating measured and calculated runoff for each month 
are plotted in Figures 8 to 13. 

There is considerably less scatter associated with this 
set of plots than the previous scatter plots of Figures 
2 to 7. However, a direct comparison is somewhat 
deceptive. The scatter in the plots is a measure of 
both the basin-to-basin variability and the difference 
between the calculated and actual probability distribu- 
tions of flow. 

The scatter associated with Figures 8 to 13 is primar- 
ily due to the basin-to-basin variability. Only the 
data within each basin was ranked. When data for all 

13 basins are lumped together, any bias is reduced 
because of the basin-to-basin variability. The coeffi- 
cient of determination {r2) ig listed in Table 3 for 
each relation. However, because of the procedure 
adopted, the values cannot be interpreted in the normal 
way. Therefore, significant values of r2 for 90 and 
95 percent confidence levels are not tabulated. The 
relations are expected to have a similar significance 
to those presented in Table 2 for the historic sequen- 
tial correlations. 



3.5 - Summary 

Physically acceptable equations have been developed for esti- 
mating natural monthly runoff in ungauged catchments of the 
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TABI£ 3 

OOEFFICIEWrS POR RANKED OORRELftTIONS 

General Equation: (April to Novenfcer) 





Precipitation 
(m) 


Apr 


<Pi 


+ 


Pi_l)/2 


May 


(Pi 


+ 


Pi_l)/2 


Jun 


<Pi 


+ 


Pi_l)/2 


Jul 


<^i 


+ 


Pi_l)/2 


Aug 


<^i 


+ 


P,.,)/2 


Sep 


<^ 


+ 


Pi.l)/2 


Oct 


'^ 


+ 


Pi=l)/2 


Nov 


«Pi 


+ 


Pi_l)/2 



2 — 

a, a_ a^ a< ^E ^c r * SEE Q , 

^ ^ ^ 4 5 6 ^^3/3)* (^3/g,. 

-1.378 0.909 0.0 0.0 -0.249 0.192 0.92 4.30 13.75 

-0.587 0.374 0.0 0.0 0.247 -0.118 0.94 1.85 6.86 

-0.001 0.171 -0.004 0.0 0.118 -0.082 0.93 0.79 2.40 

0.318 0.107 -0.003 0.0 0.079 -0.072 0.78 0.89 1.49 

0.096 0.098 -0.004 0.0 0.043 -0.043 0.84 0.58 1.13 

0.019 0.166 -0.008 0.0 -0.025 -0.025 0.79 1.22 1.75 

0.194 0.278 -0.013 0.0 -O.049 0.0 0.79 1.83 2.94 

-0.036 0.394 -0.018 0.0 0.0 0.0 0.88 1.29 3.31 



General Equation: (Deo^rber to March) 
^ = ^ + ^\*\ (aj + PCa^ + a^D)) 





Precipitation 
(ra) 


Deyiue Days 

(ooq 


^1 


^2 


^3 


^4 


^5 


r2 


SEE 
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*r , is the coefficient of determination 

*SEE (m /s) is the standard error of estimate 

*Q (m3/s) is the mean monthly discharge for all statlor^ used in the analysis 
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Canadian Shield region of central Ontario. Although equa- 
tions were developed for all 12 months, during July and 
August the relationships have limited predictive capability. 

Two types of correlation analyses were performed. These are 

- historic sequential 

- ranked. 

The relationships which were derived using historic sequen- 
tial data are applicable to 

- replicating or estimating runoff in a particular historic 
month 

- reproducing the statistical behavior of monthly runoff 
(i.e., mean, standard deviation, skew of the monthly 
distribution) . These estimates will also replicate data 
for a specific historic month. 

The relationships which were derived using ranked data are 
applicable to 

- reproducing the statistical behavior of runoff (i.e., mean, 
standard deviation, skew of the monthly distribution) . 
These estimates will not replicate data for a specific 
historic month. 

Initially it was not apparent whether the ranked or the his- 
toric sequential approach would best reproduce the statis- 
tical behavior of runoff. The test case of Section 5, which 
includes routing effects, indicates that the ranked approach 
is superior. The details of this test result are discussed 
in Section 5.4. 
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4 - FLOW PREDICTION MODEL 



The development of a hydrologic moael describing surtace 
water movement depends on both the calculation ot local 
runoff as described in section 3 and on a procedure tor 
calculating interbasin flow. The overall chain-ot-lakes flow 
prediction model therefore contains two major components. It 
includes equations for predicting local runoff and it 
includes a routing submodel. The purpose ot the routing 

submodel is to route estimates of local runoff through a 
chain of lakes which is subject to alternative types of 
regulation. A third component of the model incorporates 
calibration procedures so that the calculated runoffs and 
routed discharges can be adjusted to provide reliable 
replication of measured discharge at a basin outlet. 

The overall model may be used for two distinct applications — 
for simulating historic monthly discharges and for simulating 
the statistical behavior of the basin in order that data 
extension can be performed. The primary difference between 
these two applications is the input aata which is required 
and the calibration technique which is undertaken. A 
detailed discussion of the requirements for both is included 
later in this section. 

The remainder of this section describes the following three 
components of the model 

- natural runoff predictions 

- lake routing 

- calibration. 
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4.1 - Natural Runoff Predictions 

Local runoff is calculated using the natural runoff predic- 
tion equations which were discussed in Section 3. The equa- 
tions which are currently in the model are Equations 12 and 
13. The coefficients a^ to a^ of these linear relations 
are read in as input data. These coefficients are listed in 
Tables 2 and 3 for the historic sequential and the ranked 
rainfall runoff correlations. If, in the future, the regres- 
sion analyses of Section 3 are redone it is possible that 
equations with an alternative form, such as exponential or 
logorithmic will be derived. In this case it will be 
necessary to alter the Fortran coding for one subroutine of 
the model. 



4.2 - Lake Routing 

Interbasin flows are estimated by the use of the routing 
submodel presented in this subsection. Assumptions and 
simplifications of the routing submodel are discussed 
followed by a presentation of the subbasin discharge calcula- 
tions. The subsection on subbasin discharge calculations 
includes a discussion of the two methods which are available 
within the model for estimating lake storage effects. 



4.2.1 - Assumptions and 
Simplifications 



The assumptions and simplifications which form the 
basis of the routing submodel are listed below. 

- Channel routing and channel storage are neglected. 
This is based on the assumption that lake storage ana 
lake routing predominate in regulating flows over a 
1-month period. 
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- Lakes which do not have signiticant storage effects 
are not included in the routing procedure, but rather 
are handled implicitly through their inclusion in the 
runoff prediction equations. 

- The model allows for one regulating lake per sub- 
basin. This lake must be located at the downstream 
end of the subbasin. 

- Outlet flows from lakes are modeled as either natural 
outlets or as regulated outlet structures. 



4.2.2 - Subbasin Discharge 
Calculations 



The subbasin discharge is calculated using basic 
continuity principles. The discharge from a subbasin 
is simply the sum of the following components. 

- Local runoff estimated using the runoff prediction 
equations. When calculating the local runoff, the 
area of the regulated lake is excluded from both the 
area of basin and the area of lakes, 

- Discharge entering the subbasin from one or more 
upstreaun subbasins. 

- Net precipitation falling directly on the regulated 
lake. Net precipitation is the difference between 
monthly precipitation and monthly evaporation where 
evaporation is assumed to equal the long-term average 
for that month. 

- Storage changes in the regulating lake. 
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It the subbasin outlet is selected such that no lake is 
at the outlet, the continuity equations reduce to local 
runoff estimates plus upstream inflows. Storage 
changes within each of the regulating lakes are esti- 
mated using either one of two techniques. The two 
techniques for estimating storage changes are outlined 
below. 

(a) Regulation Subject to 
Rule Curve Operation 

A rule curve for a particular lake is an operating 
policy which specifies an objective time series of 
lake elevations (normally a 1-yr cycle). Storage 
changes for rule curve operation assume that the 
lake is operated to match objective lake eleva- 
tions. After accounting for inflows, water is 
discharged from the subbasin so as to ensure that 
each month-end lake level does not exceed the 
level specified by the rule curve, Vvhere insuffi- 
cient water supply exists, the end of month lake 
level is allowed to drop below the rule curve 
elevation. In this case the starting water level 
for the next month's storage calculation is the 
previous month's actual ending water level and not 
the rule curve level. 

Under some circumstances historic monthly reser- 
voir levels are available. These actual lake 
levels should be used instead of rule curve levels 
when they are available. The model is designed to 
accept these values as input, however, tor actual 
operation of the model the historic lake levels 
are treated as if they were rule curve values. It 
is always preferable to use historic lake levels 
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if they are available instead of rule curves for 
historic simulations. 

When using the model for statistical simulation 
purposes, rule curves are adequate, providing they 
reflect historic reservoir operation, simple 
calculations can be made outside of the moael to 
derive average "historic" rule curves rather than 
"objective" rule curves. The two do not always 
coincide. 

(b) Regulation Subject to 
a Rating Curve 

A rating curve for a particular lake is an eleva- 
tion discharge relation tor the lake outlet. 
Normally a rating curve will apply only to natural 
lakes but it will also apply to lakes with non- 
operated man-made structures. Rating curve stor- 
age estimates are made in the following manner. 

- Inflow is assumed constant throughout the month 
and is divided into 10 equal segments. 

- Each equal flow increment is routed through the 
lake, with the outflow during each time step 
governed by the lake elevation of that time 
step. 

- The resulting 10 outlet aischarges are averaged 
to give the mean monthly discharge. 

A monthly subbasin discharge is produced by one of the 
above techniques in conjunction with the general 
continuity equation. This flow is then applied as an 
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inflow to the next immediate downstream subbasin. The 
procedure is repeated for each subbasin until monthly 
network discharges are calculated at all points. 

Details of the discharge calculation procedure can be 
found in Appendix C. 



4.3 - Calibration Procedures 

The chain-of-lakes model consists of subbasin runoff calcula- 
tions, followed by flow routing from subbasin to subbasin. 
The discharge at the downstream end of the drainage network 
can optionally be compared within the model to measured 
discharge records. Typically, differences will exist between 
the measured and calculated discharges. 

Calibration is required to make the model basin specific. 
The following two calibation techniques have been incorpo- 
rated into the model. 

- Historic replication of a specific basin discharge record. 

- Improvement of the runoff regression equations for a 
specific data record. 

Each calibration technique is intended for a different pur- 
pose and consequently there are application limitations to 
each technique. The techniques and their limitations are 
discussed in the following subsections. Details of each are 
given in Appendix D. 
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4.3.1 - Historic Replication 

Calibration using historic replication of a specific 
basin discharge record adjusts the subbasin runoff 
values on a raonth-by-month basis so that the total 
calculated basin discharge is equal to the measured 
basin discharge. This method Is event specific and 
therefore is limited to application only for periods 
where historic discharge data is available. The sub- 
basin runoff values are adjusted by prorating on the 
basis of subbasin runoff, as computed by the regression 
techniques. This approach is simply an extension and 
improvement of the common technique of prorating dis- 
charges based on subbasin area alone. Although this 
approach is physically realistic, it is not possible to 
judge its accuracy using the available data base. This 
calibration technique is intended for application on a 
basin which has a long historic record of measured 
discharges at a arainage basin outlet. The calibration 
will ensure that calculated and measured discharge 
match perfectly. In addition, it will be possible to 
estimate the origin of these flows from upstream sub- 
basins based on prorating subbasin calculated runoffs 
(from regression equations) rather than simply from 
subbasin drainage areas. 

To replicate the historic basin discharge, the follow- 
ing steps are followed by the model. 

1 - The subbasin runoff is computed on a month-by-month 

basis for all subbasins. These values are summed 
to provide an overall basin runoff estimate. 

2 - The subbasin runoffs are routed through the chain- 

of-lakes network and a calculated discharge at the 



basin outlet is derived. For comparative purposes 
this initial calculated basin discharge is output 
along with the measured basin discharge. 

3 - The total measured basin runoff is approximated by 

backrouting measured discharge with respect to 
changes in lake storage. Adjustments are also made 
for net lake evaporation. 

4 - A ratio is computed by dividing the total measured 

basin runoff (Step 3) by the sum of the subbasin 
runoff (Step 1) . 

5 - The calculated runoff for each subbasin (Step 1) is 

then multiplied by the ratio which was computed in 
Step 4. 

6 - The adjusted subbasin runoffs are routed through 

the drainage network, 

7 - If additional iteration is necessary, Steps 3 to 6 

are repeated. 

The calculated basin discharge for each month will 
implicitly converge to the measured basin discharge. 
More than one iteration may be required because the 
total volume of stored water will normally vary from 
iteration to iteration. 

4.3.2 - Improvement to Runoff 
Prediction Equations 

Instead of replicating specific runoff on a month-by- 
month basis, an improvement to the runoff prediction 
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equations can be made within the model. This improve- 
ment to the runoff prediction equations has a broad 
basis for application and can be used for both historic 
simulation and for extension of discharge records 
beyond the period of available data. 

Revised runoff prediction equations are derived for the 
overall basin based on actual measured basin discharge. 
The measured discharge is backrouted through the 
regulated lakes to remove the effect of lake storage. 
Monthly correlations are then made between these 
naturalized flow records and the appropriate climato- 
logic data. The form of these equations is identical 
to Equations 14 and 15 of Section 3.4. The result is a 
set of 12 revised runoff prediction equations which are 
calibrated for the specific basin under study. These 
revised equations are then used to calculate an overall 
basin runoff on a month-by-month basis. In the same 
manner as was done for Section 4.3.1, ratios are 
calculated between revised basin runoff and initial 
basin runoff. These ratios are then multiplied by the 
initial individual subbasin runoff estimates. A 
detailed description of this calibration technique is 
presented in Appendix D. The general methodology is as 
follows. 

1 - For the entire simulation period, subbasin runoffs 

are calculated by using the runoff prediction 
equations without any adjustment factor. These 
values are summed to provide an overall basin 
runoff estimate. 

2 - For the entire simulation period, the subbasin 

runoffs are routed through the drainage network and 
calculated discharges are derived for the basin 
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outlet. For comparative purposes these calculated 
basin discharges are output along with the measured 
basin discharge. 

3 - The total measured basin runoff is approximated by 

backrouting the measured discharge with respect to 
changes in lake storage. Adjustments are also made 
for net lake evaporation. 

4 - For each of the 12 months, linear equations are 

developed which describe the total measured basin 
runoff in terms of the climatologic data. The 
equations are used only to calculate the total 
basin runoff. 

5 - On a month-by-month basis, the total basin runoff 

is calculated by using the revised regression 
equations, as developed in Step 4. 

6 - On a month-by-month basis, adjustment ratios are 

calculated by dividing the initial predicted basin 
runoff (Step 1) into the adjusted basin runoff 
(Step 5) . 

7 - The individual subbasin runoffs (Step 1) are then 

multiplied by the appropriate ratio as determined 
in Step 6. 

8 - For the entire simulation period, the subbasin 

runoffs are routed through the drainage network. 

9 - If additional iteration is necessary. Steps 3 to 8 

are repeated. 
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The total basin average monthly discharge, as computea 
in the above calibration technique, will implicitly 
converge toward the corresponding average for measured 
discharge. More than one iteration may be required 
because the total volume of stored water will normally 
vary from iteration to iteration. This variation in 
volume of stored water from iteration to iteration will 
only occur for the rating curve storage assumptions or 

when lake levels do not reach rule curve levels. If 
all lakes are operated to rule curve and all lakes 
achieve this objective level then only one iteration 
will be required for perfect closure. 

The results of the calibration can be used for subse- 
quent data extension runs of the model. The results of 
the linear runoff equations (Step 4) can be read into 

the model as input data. Steps 1, 2 ana 5 to 8 are 
then repeated within the model to produce basin 
specific subbasin and total basin discharge values. 



4.4 - Model Application 

This section presents an overview of necessary steps and data 
requirements for running the flow prediction model. A 
description of output options is included. 

4.4.1 - Required Input Data 

Details of the data input requirements are given in 
Appendix E. 

Prior to collection of pertinent data it is necessary 
to delineate the basin and subdivide it into the 
desired subbasins. Guidelines for this task are as 
follows. 
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- Any stream network to be simulated should preterably 
incorporate the drainage area above a streamflow 

gauge. It is desirable, although not necessary, to 
include a known gauge as the downstream node of the 
network in order to take advantage of the calibration 
procedures within the model. If additional stream- 
flow gauges exist within the network, then the drain- 
age area should be divided into separate networks tor 
each gauge. 

- Individual subbasins should preferably not exceed 
about 750 km^. use of larger basins requires using 
the runoff prediction equations beyond the range for 
which they were developed. 

- Subbasin nodes should be located at the outlet of all 
regulated lakes, or of all lakes which have signifi- 
cant natural regulation. 

- Subbasin nodes should also be located at any point 
where a flow prediction is required. 

- The subbasins must be assigned a nodal number. Any 
numbering system can be assigned because the model 
does not make assumptions on the numbering hierarchy. 
The interconnection of nodes is defined through the 
input data. 

The required data for each of the three components of 
the model is discussed below. 

(a) Natural Runoff Predictions 

- Climatologic data such as precipitation and 

degree days above O'C for each subbasin based on 
the nearest or most appropriate station. 



- Evaporation data for each month. 

- Physiographic data such as areas of subbasins, 
lakes and swamps, 

- Coefficients of the runoff prediction equations 

Oj to ag) for either the historic sequen- 
tial or ranked rainfall/runoff correlations. 

( b) Lake Routing 

- Reservoir rule curves (either objective or 
historic) . 

- Historic month-end lake elevations. 

- Outlet descriptions for natural lake outlets or 
nonoperated lake structure outlets. 

( c) Calibration 

- Historic monthly discharges at the outlet to the 
overall basin. 

In addition there are numerous flags and codes that 
must be input to control the selection of computational 
options. Details of these codes and flags can be found 
in Appendix E. 

4.4.2 - Output Selection 

There are three basic output options available within 
the model. They are described briefly below. 
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(a) Detailed Month-by- 
Month Output 



This output lists details for all subbasins, one 
page for each month of simulation. Details of 
input data, lake level changes and calculated dis- 
charge are provided tor each step of the computa- 
tions. This output can be used for assessing the 
relative significance of various components of the 
model or of particular subbasins within the 
network. 

( b) Discharge Summary 

This output provides a discharge summary for the 
simulation, one page for each subbasin. In 
addition, for calibration runs, a summary table is 
provided which lists measured and calculated 
discharge at the network outlet. This summary is 
provided for each iteration of the calibration 
procedure (normally 2 or 3 iterations). 

( c) Lake Level Summary 

This output provides a lake level summary for all 
months of simulation, one page tor each subbasin 
(i.e., each major lake). This output provides 
infoinnation on lake level variations and indicates 
whether or not objective rule curve levels were 
achieved. 

There are two other forms of output which are always 
produced and which are intended for use as input data 
for future computations. These are 
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- output in the correct format to be read by the 
existing serpent River water quality model 

- output in the correct format to be read by this 
chain-of-lakes flow prediction model as an input to a 
future downstream network. 

Both of these output forms are sent to separate 
computer output units. Details of the output units can 
be found in Appendix E and example output can be found 

in Tables F-8 and F-9 of Appendix F. 



5 - TEST RESULTS 
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5 - TEST RESULTS 



The results of calibrating and running the model are pre- 
sented in this section in order to verify correct operation 
and to demonstrate an application of the model. 

The South Branch Muskoka River was selected tor testing the 
chain-of-lakes model. This basin was selected because of 
Acres' familiarity with the region, and because the required 
data was readily available. 

A description of the basin and its operation, including data 
availability and limitations, is presented in this section. 

The test results are presented for the following two cases 

- month-by-month historical replication 

- statistical simulation. 

A discussion of the accuracy and acceptability of the results 
for each test is included in Sections 5.3 and 5.4. 
Sections 5,1 and 5.2 are a description of the South Branch 
Muskoka River and the lake levels used in the test runs. 



5.1 - Description of the 
South Branch 
Muskoka River 



The South Branch Muskoka River is a subbasin ot the overall 
Muskoka drainage basin which drains an area of 4,770 kni2 ^ 
stretching from Algonquin Provincial Park to the outlet ot 
Lake Muskoka near Georgian Bay. The south branch drains an 
area of approximately 1,500 km2. This study only used the 
drainage area of 1,390 km2^ located upstream from the WSC 
streamflow gauge located at Baysville (2EB008). 
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The South Branch Muskoka River is characterized by six major 
lakes subject to regulation policy, as well as several 

smaller lakes regulated by nonoperative structures. Follow- 
ing the guidelines of Section 4.1.1 the overall network was 
divided into six subbasins. The outlet of each of the six 
major lakes represents the outlet of each of the six sub- 
basins. Figure 14 is a schematic of the network. 

The water levels of all six lakes are subject to rule curve 
operation. The water levels of Lake of Bays and of Lake 

Kawagama are monitored daily. Adjustments to log settings 
are made as necessary. The water levels of the remaining 
four lakes are not monitored on a daily basis. As well, the 
log settings are only changed a few times a year. 

Daily water levels are available for Lake Kawagama dating 
from 1975. For Lake of Bays, data for the period dating from 
1973 is presently on hand. 

Discharge records for the streamflow gauge at Baysville date 
from 1942; however, because of the limitations of the other 
data bases, the period from 1970 to 1981 was selected. The 
discharge records are presented in Table 4. 

Meteorologic data for Huntsville (1970 to 1981) was selected 
as representative of the South Branch Muskoka River. The 
data are presented in Table B-3 of Appendix B. Occasionally, 
data for Huntsville is missing. For the missing months, 
meteorological data for Dwight is used as a substitute for 
the missing Huntsville data. 

The physiographic data were taken from 1:50 000-topographical 
maps and from records supplied by the Ontario Ministry of 
Natural Resources. 
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S8.S 


34.7 


31.9 


17.6 


16.5 


15.9 


16.8 


10.7 


23.5 


29.6 


1974 


25.9 


33.3 


40.6 


57.1 


64.7 


19.7 


7.6 


6.8 


12.3 


26.1 


28.4 


22.5 


28.7 


1975 


20.3 


23.1 


40.1 


35.5 


31,5 


4.6 


3.4 


5.1 


13.8 


20.1 


10.4 


32.0 


20.0 


1976 


25.0 


27.6 


50.5 


69.1 


24.0 


8.5 


18,6 


6,3 


6.2 


11.8 


8.2 


12.9 


22.4 


1977 


15.4 


13.S 


41.5 


38.7 


13.9 


4,0 


J.O 


3.3 


14.1 


33.7 


24.1 


27.5 


19.5 


1978 


29.7 


19.5 


27,7 


25.1 


53.4 


11,3 


4.0 


6,3 


19.7 


25.6 


14.7 


21.4 


21.6 


1979 


21.5 


30.7 


50.1 


57.4 


44.7 


15.8 


S.l 


9.5 


22.1 


28.3 


43.1 


40.1 


30.6 


1980 


30.2 


25.1 


23.0 


78,5 


20.6 


22,5 


17.0 


22.8 


24.3 


33.9 


30.3 


24.3 


29,3 


1981 


20.9 


37.1 


52.6 


31.8 


22.3 


13,6 


5.2 


4.0 


61.3 


37.5 


14.6 


16.7 


26.4 
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5.2 - Lake Level Estimates 

When calibrating and running the model, four alternatives are 
generally available for accounting for changes in lake eleva- 
tions. One method utilizes rating curves. It is not 
discussed here since it is not applicable to the South Branch 
Muskoka River. The three remaining methods are 

- rule curves 

- historic rule curves based on average recorded lake levels 

- historic lake levels. 

An examination of historic lake levels indicate that the 
major lakes (Lake Kawagama and Lake of Bays) were not oper- 
ated exactly to rule curve. Consequently, average historic 
lake levels were used to derive average historic rule curves. 
When using the model for statistical simulation, these 
average historic rule curves for Lake Kawagama and Lake of 
Bays were used in place of the objective rule curves. The 
objective rule curves were used for the four remaining 
regulated lakes. However, when calibrating and running the 
model on a month-by-month basis, historic lake levels were 
used for Lake Kawagama and Lake of Bays. These levels were 
used in order to take advantage of the best available data. 



5.3 - Historic Replication 

The model was run for the period 1975 to 1980. The historic 
replication calibration procedure was selected tor use in the 
model to ensure that measured and calculated discharge at the 
downstream node of the basin matched perfectly. Table 5 
lists the basin runoff provided initially by the prediction 
equations (sura of all six subbasins) and the final basin 
runoff after calibration. 
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TABLE 5 COMPARISON OF CALIBRATED AND UHCALIERATED 
BASIN RUNOFF FOR HISTORIC REPLICATION OF 
THE SOUTH BRANCH «USKOKA RIVER 

Discharae in Cubic Meters per Second 



YEAR CALIBRATION JAN FEB MAR APR HAY JUN JUL AUG SEP OCT NOV DEC 



1975 m 10.4 20.4 21.0 44.1 46.7 16.8 11.1 8.0 7.6 ?.? 19.0 36.9 
YES 14.2 11.4 18.9 40.7 41.7 10.9 6,6 4.1 4.6 S.t 12.1 30.3 

1976 NO 11.2 7.5 133.7 49.5 45,7 19.2 10,3 6.4 5.9 12.5 24.7 19.0 
YES 13.5 9.6 51.7 93.0 26.0 11.9 20.0 5.1 -3.4 2.1 2.8 9.9 

1977 NO 9.4 7.2 35.3 42,1 45,6 15.1 8.8 8.6 9.2 18.4 32.1 15.8 
YES 2.6 5.0 43.3 68.3 17.6 6.6 5.3 3.5 4.1 23.0 26.1 23.0 

1978 NO 9.4 3.3 11.4 32.1 46.0 18.3 12,1 11,3 11,7 18.4 20.2 13.3 
YES 14.2 8.7 10.8 55.9 69.3 14,4 6.4 4.5 6.6 17.0 16.7 18.9 

1979 NO 11.1 7.6 45.4 67.3 47.6 21.2 14.4 8.0 6.6 19.2 21.7 22.4 
YES 13.5 9.7 45.5 91.1 41,9 10.2 15.7 10.8 5.1 22.8 43.1 33.5 

1980 NO 9.2 6.6 62.2 56.7 47.8 23.4 17.7 13.6 11.5 21.S 1^3 18.1 
YES 23.7 9.2 24.4 99.9 22.9 26.8 20.4 14.3 16.9 22.0 31.8 16.0 
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The results of the simulation illustrate that the model can 
be run to produce historic discharges at the downstream node. 
However, the initial runoff prediction may have to be 
adjusted considerably to produce the desired result (see 
Table 5) . The magnitude of the adjustments may be attributed 
to the following causes: 

- the runoff prediction equations are not basin specific 

- climatological data from Huntsville may not adequately 
represent precipition and temperature in the basin of the 
south branch of the Muskoka River on a month-by-month 
basis 

- there may be errors in the historic lake levels recorded 
for Lake Kawagama and Lake of Bays 

- the rule curves may not always represent actual lake levels 
in each of the four smaller regulated lakes. 

Detailed output for this test run is presented in Appendix F, 
including input data files and output listings. 



5.4 - Statistical Simulation 

Two separate statistical simulation test runs were made. The 
runs were made using both the historic sequential and ranked 
rainfall/runoff correlations. In both cases the model was 
run for the period 1970 to 1981 using historic rule curves 
for Lake Kawagama and Lake of Bays. The calibration proce- 
dure of Section 4.3.2 was applied to improve the runoff 
prediction equations and make them basin specific. 
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Table 6 is a presentation of the results using the ranked 
runoff prediction equations. The measured basin discharges 
and the calculated basin discharges are listed in the table. 
Since the ranked approach was used, each measured discharge 
is paired with the calculated discharge of similar rank. At 
the bottom of the table the means, standard deviations and 
coefficients of skew are listed for both the measured and 
calculated data sets. The calculated and measured means are 
equal as required by the calibration. The standard aevia- 
tions are virtually identical in all months while the coeffi- 
cients of skew are not at all similar. However, 12 data 
points are not sufficient to adequately evaluate skewness. 
The standard error between ranked measured and ranked 
calculated discharge is also given in Table 6. 

An identical run was made using the historic sequential run- 
off prediction equations. The calculated and measured means 
are again equal. Table 7 lists the calculated standard 
deviations. In addition. Table 7 lists the measured standard 
deviations and the standard deviations which were calculated 
using the ranked runoff prediction equations. 

The results which are presented in Table 7 verify that the 
ranked approach produces better statistical results than the 
historic sequential approach. Therefore, it can be concluded 
that the ranked approach is probably better suitea than the 
historic sequential approach for extending discharge data 
beyond the period of available data. However, this conclu- 
sion is valid provided that only the statistical behavior of 
runoff is required and that reproduction of specific historic 
monthly discharge is not required. If a reproduction of 
historic monthly discharge records is required, then only the 
historic sequential approach can be used tor data extensions. 



TABLE 6 COttPARISOf* OF MEASURED AND CALCULATED DISCHARGE AT 
BAYSVILLE USI«G THE RANKED CORRELATION TECHNIQUE 



YEAR 


TYPE 
CALC 


JAN 
17.2 


FEB 

22.2 


MAR 
32.0 


APR 
46.6 


HAY 
55.0 


JUN 
17,4 


JUL 
37.5 


AUG 
18.4 


SEP 

11.1 


OCT 


NOV 


DEC 
28.7 


MEAN 


1 


26.4 


16.2 


27.4 




MEAS 


13.4 


13.8 


27.7 


47.5 


64,3 


19.7 


49,? 


20.4 


12.3 


28.3 


14.6 


32.0 


28.8 


2 


CALC 


21,9 


27.1 


31.8 


26.3 


17.2 


10.2 


6.0 


6,3 


1.6 


6.7 


15.4 


28.6 


16.6 




MEAS 


20.9 


30.7 


23.0 


28.3 


13.9 


8.5 


5.1 


5.8 


11.2 


7.6 


8.1 


27.5 


15.9 


3 


CALC 


31.2 


24.6 


41.6 


46.4 


44.6 


14.2 


19.3 


18.1 


15.9 


25.2 


15.9 


20.4 


26.5 




MEAS 


30.2 


27,3 


50.1 


40.1 


34.7 


18.5 


17.6 


16.5 


13.8 


24.7 


10.7 


21.4 


25.5 


4 


CALC 


28.2 


37.5 


57.3 


51.9 


50.6 


35,3 


33,2 


0.0 


0.0 


13.6 


17.5 


25.3 


1-9.2 




KEAS 


29.5 


37.1 


52.6 


57.4 


53.1 


31.9 


26.1 


4,0 


6.2 


11.8 


21.3 


25.5 


29.7 


5 


CALC 


29.9 


22.8 


32.4 


68.5 


62.4 


11.2 


0.0 


6.2 


26.6 


32.2 


19,3 


23.5 


27.9 




HEAS 


29.7 


24.0 


31.6 


69.1 


64.7 


11.3 


3.0 


5.1 


22.1 


33.7 


30.3 


24.3 


29.1 


6 


CALC 


24.1 


36.7 


33.8 


41,7 


31.2 


8.6 


6.7 


3.8 


17.3 


18.1 


17.1 


11.0 


23.7 




HEAS 


25.0 


33.3 


36.1 


38,7 


31.5 


4.6 


5.2 


6.3 


14.4 


16.8 


14.7 


40.1 


22.2 


7 


CALC 


23.2 


25.4 


61.8 


51.5 


19.8 


12.6 


2.1 


0.0 


21.3 


25.3 


18.8 


18.2 


23.3 




MEAS 


21.5 


27.4 


57.3 


57.1 


24.0 


13.6 


3.1 


3.3 


16.0 


25.6 


24,1 


16.7 


24.2 


8 


CALC 


20.4 


22.6 


35.1 


37.4 


17.3 


7.2 


4.5 


10.7 


23.7 


26.2 


22,4 


9.6 


19.7 




HEAS 


20.3 


23.1 


40.6 


31.8 


20,6 


4.0 


4,0 


6,8 


19.7 


26,1 


31.9 


10,4 


19.9 


9 


CALC 


18.3 


22.8 


29.9 


20.2 


19,7 


11.2 


12.3 


17,9 


26.9 


20.3 


15.6 


15.0 


19.2 




MEAS 


19.;. 


23.1 


22.3 


25.1 


22.3 


11,1 


17.0 


9.5 


24.3 


20.1 


10.4 


12.9 


18.3 


10 


CALC 


24.7 


22.3 


40.0 


38.8 


44,8 


13.1 


6.3 


8.8 


14.3 


33.0 


19.2 


26.7 


24.3 




MEAS 


25.9 


19.5 


41.5 


35.5 


43.6 


13.8 


5.2 


6.3 


13.2 


33.9 


28.4 


25.6 


24.4 


11 


CALC 


20.1 


22.6 


42.7 


70.7 


43.5 


19.8 


26,7 


18,6 


21.0 


30.2 


14.2 


21.9 


29.3 




HEAS 


20.2 


22.2 


50.5 


78.5 


33.6 


22.5 


18.6 


22.8 


15.9 


30.9 


8.2 


22,5 


28.9 


12 


CALC 


18.4 


38.0 


35.1 


67.9 


45.1 


13.9 


8.2 


25.2 


51.1 


39.6 


54.4 


21.7 


34.9 




HEAS 


19.7 


40.9 


40.1 


58.8 


44.7 


15,8 


7.6 


31.9 


61.3 


37.5 


43.1 


21.8 


35.3 


MEAN 


CALC 


23.1 


27.1 


39.5 


17.3 


37.6 


14.6 


13.6 


11,6 


19.2 


24.7 


20.5 


23.4 






HEAS 


23.1 


27.1 


39.5 


47.3 


37.6 


14.6 


13.6 


11.6 


19.2 


24.7 


20.5 


23.4 




STANDARD 


CALC 


4.7 


6.1 


10.3 


16.0 


15.9 


7.5 


12 „S 


3.0 


13.3 


9.0 


10.9 


7.9 




DEVIATION 


HEAS 


4.8 


7.6 


11.7 


16.9 


16.3 


7.8 


13.8 


9,2 


14.3 


9.2 


11.2 


8.0 




COEFFICIENT CALC 


0.5 


1.2 


1.5 


0.0 


-0.1 


2.3 


1.0 


0,1 


0,9 


-0.5 


3.2 


0.5 




OF SKE« 


HEAS 


0.3 


0.3 


-0.1 


0.4 


0.4 


0.8 


1.9 


1.2 


2.3 


■0.6 


0,7 


0.4 




STANDARD 




1.2 


3.6 


6.3 


5.7 


i.9 


2,8 


5.6 


4.4 


6.0 


1.3 


7.9 


1.6 




ERROR 































NOTE! All values are in units of cubic ielers per second with the 
exception oT the coefficient of skew which is diuensionless. 
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TABLE 7 COHPARISON OF MEASURED AND 

CALCULATED STANDARQ OEVIATIOMS 



MONTHLY STANDARD DEVIATION HONTHLY STANDARD DEVIATION 

ffONTH (Cubic Heters/Second) MONTH (Cubic Meters/Second) 

Mess Csk-l C3lc-2 Mess Csk-l C3k-2 



JAN 


4.S3 


4.65 


2.28 


JUL 


13.76 


12.61 


10.20 


FEB 


7.57 


6.43 


4.83 


m 


%& 


8.02 


5.72 


flnn 


11.71 


10.25 


5.1S 


SEP 


14.27 


13.29 


S.34 


APR 


14.89 


16.05 


9.27 


OCT 


hm 


9.03 


7.32 


mm 


16.83 


15.90 


4.16 


NOV 


11.20 


10,92 


5. 89 


m 


7.85 


7.55 


7.04 


m 


§M 


7,87 


2.65 



(^k-1 SUndard Deviations usind ranked regression eouations 
C8lc-2 SUridard Devl3tions usind seouentisl re<iression eauations 
Hess Standard Deviations ysind leasured discharge data at Bsvsville 
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APPENDIX A 
SUMMARY OF HYDROMETRIC DISCHARGE DATA 



A-1 



TftBLE fl-1 (IHITStW RIVER AT CHELMSFORD 
Station 2CF0O7 

Discharge in Cubic Metres per Second 



YEAR JAN FEB WAR APR MAY JUN JUL AUG SEP OCT NOV DEC 



1974 


— 


— 


— 


™ 


~ 


~ 


— 


— 


~ 


'- 


0.64 


0.53 


1975 


0.3A 


0.3S 


0.46 


6.89 


1.59 


0.98 


0.51 


0.42 


0.48 


0,65 


0.36 


0.72 


1976 


0.39 


0.28 


1.4? 


8.07 


2.30 


0.99 


1.02 


0.61 


0.48 


0.36 


0.34 


0,25 


1977 


0.17 


0.18 


1.25 


5.57 


0.76 


0.41 


0.40 


0.36 


0.53 


1.02 


1.08 


0.32 


1978 


0.26 


0.29 


0.43 


5.83 


2,54 


0.51 


0.50 


0.44 


0.45 


1.86 


0.81 


0,70 


1979 


0.39 


0.31 


2.00 


10.20 


3.27 


0.75 


0.81 


0.48 


0.11 


0.39 


2.04 


0.99 


1980 


1.00 


0.43 


0.90 


8.08 


3.20 


0.61 


0.40 


0.36 


0,94 


1.73 


0,79 


0.33 



1981 0.33 1.26 2.54 5.67 1.41 0.75 0.80 0.37 0.92 1.65 1.05 0.42 



TABLE A-2 UHITSON RIVER AT VAL CARON 
Station 2CF008 

Discharse in Cubic Metres per Second 
YEAR JAN FEB IttR APR iWY JUN JW. m SEP OCT NOV DEC 



1970 


~ 


— 


— 


7.82 


— 


~ 


— 


~ 


~ 


~ 


— 


-_ 


1971 


— 


— 


— 


11.30 


— 


~ 


— 


~ 


— 


— 


— 


— 


1972 
1973 
1974 


~ 


~ 


-- 


3.47 
7.27 


— 


— 


~ 


— 


~ 


— 


3.34 


0.83 


1975 


0.81 


0.56 


0,73 


8.01 


2.94 


0.73 


0.10 


0.10 


0.20 


0.75 


1.09 


1.78 


1976 


0.26 


0.25 


1.62 


8.23 


2.23 


0.69 


0.21 


0.18 


0.90 


1.04 


0.70 


0.34 


1977 


0.29 


0.27 


3.83 


6.63 


0.84 


0.35 


0.60 


0.67 


1.73 


3,10 


2.74 


1.65 


1978 


0.91 


0.62 


0.62 


6,67 


2.99 


0.63 


0.70 


0.98 


1.81 


4,20 


1,57 


1,15 


1979 


0.86 


0.62 


2.62 


15.10 


3.62 


0.89 


0.59 


0.35 


0.93 


1.31 


2,12 


1.85 


1980 


1.31 


0.54 


0.72 


9.32 


3.03 


0.86 


0.52 


0.56 


1.32 


2.51 


1.53 


0.75 


1981 


0.36 


1.68 


2.89 


7.13 


1.32 


0.79 


0.56 


0.23 


1.87 


2.63 


1.30 


0.80 



A-2 



TABLE A3 VERMILLION RIVER NEAR CAPREOL 
SUtion 2CF100 

Discharte in Cubic Metres per Second 



YEAR JAH FEB »AR APR HAY JUN JUL AUG SEP OCT NOV DEC 
1970 - - 29.10 38.90 11.50 6.12 4.9i 10.20 8.29 7.90 



iy/1 


"■■" 


"^ 


"" 


"■" 


iiB./g 


1«.W 


J. 3d 


l.HY 


1,13 


J.UD 


4.U.: 


t,n^ 


1972 


4.78 


3.52 


3.61 


15.90 


35.90 


7.33 


4,69 


7,08 


7.20 


6,83 


8.25 


5.38 


1973 


4,54 


4.38 


13.60 


21.50 


20.90 


12.60 


7.25 


6.14 


3.87 


4.92 


7.16 


6.05 


1974 


4.60 


3.92 


3.66 


21.30 


21.00 


11.30 


5.59 


3.33 


2.17 


3.82 


10.80 


6.05 


1975 


— 


— 


-- 


— 


24.90 


11.90 


5,35 


2.25 


1.67 


2,32 


4.58 


8.84 


1976 


— 


— 


— 


48.20 


26.S0 


9.05 


3.14 


1.49 


1,86 


2.44 


— 


~ 


1977 


— 


~ 


~ 


~ 


11.20 


3.75 


3.07 


2.06 


3.86 


8.21 


8.92 


7.27 



TABLE A-4 DUCHESNAY RIVER NEAR NORTH BAY 
SUilon 2BD008 

Biscfiarse in Cubic Metres per Second 



YEAR JAN FEB «Aft APR NAY JUN JUL AUG SEP OCT NOV DEC 
1970 0.20 0.18 0.25 6,22 4,25 3.26 5.86 1,19 1.98 1.98 1.61 0.74 



1971 


0.32 


0.33 


0.47 


6.58 


3.01 


0.30 


0.11 


0.16 


0.43 


0.66 


0.88 


1.84 


1972 


0.69 


0.37 


0.32 


5.31 


5.95 


1.14 


0.63 


3.71 


0.66 


1.47 


1.49 


0.68 


1973 


0.52 


0,52 


3.78 


6.27 


2.96 


1,52 


0.41 


0.80 


0.46 


1.39 


1.72 


1.11 


1974 


0.52 


1.06 


0.54 


7.33 


3.09 


2.70 


0.52 


0.47 


0.45 


1.49 


4.21 


1.00 


1975 


0.50 


0.29 


0.55 


4.47 


4.11 


i.74 


0.41 


0.11 


0.31 


1.25 


2,02 


2.84 


1976 


0.40 


0.19 


1.84 


6,79 


2.95 


0.64 


1,44 


0.39 


0.91 


0.89 


0,51 


0.37 


1977 


0.25 


0.23 


2.10 


7.17 


0,69 


0.38 


1,63 


1.55 


1.77 


1.50 


2,03 


0.72 


1978 


0.40 


0.43 


0.27 


4.21 


4.13 


0.57 


0.54 


1.31 


2.38 


2.49 


1.19 


0.40 


1979 


0.44 


0.43 


3.27 


8.91 


2.90 


0.74 


0,58 


0.58 


0.96 


1.71 


1.90 


1.53 


1980 


0.69 


0.29 


1.76 


8.50 


1.80 


1.08 


1.18 


0.63 


2,79 


2.46 


1.60 


0.65 



1981 0.32 1.50 1.56 7.32 2.15 1.34 0.31 0.30 1.93 3.55 1.16 0.67 
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A-3 



TABLE fl-5 VEUVE RIVER NEAR VERNER 
Ststion 2DD012 

Discharge in Cubic Metres per Second 



YEAR J« FEB BAR APR HAY JUN JUL AUG SEP OCT NOV DEC 



1973 


— 


— 


— 


— 


21.60 


5.75 


1,79 


5.92 


0.69 


2.28 


8.83 


10.30 


1974 


2.81 


0.88 


1.33 


36.10 


13.80 


11.20 


1.43 


0.74 


2.81 


12.90 


30,80 


5.70 


1975 


3.49 


1.64 


1.70 


49.10 


20.30 


5.55 


0.75 


0.41 


0.30 


0.91 


4,04 


13.60 


1976 


1.87 


1.25 


20.20 


68.30 


23.20 


3.26 


0.84 


0,92 


1.41 


2.39 


4.77 


2.42 


1977 


1.35 


1.18 


34.10 


40.40 


3.57 


0.52 


1.59 


5.70 


17.60 


19.00 


23.60 


13,20 


iv/a 


4.15 


2.31 


1.93 


50.60 


25.00 


1.98 


0.97 


1.84 


14,90 


25.50 


9.65 


6.10 


1979 


3.55 


2.78 


29.40 


89.10 


18.90 


5.73 


2.79 


2.08 


2.73 


5.15 


21.80 


11.00 


1980 


9.02 


1.97 


10.00 


68.90 


18.40 


3.22 


2.26 


1.82 


10.80 


18.40 


11.80 


4.74 



1981 2.65 16.10 17.20 61.00 12.20 8.73 4.21 1.11 8.99 17.60 8.86 5.71 



TABLE A-6 LA VASE RIVER AT NORTH BAY 
SUtion 2DD013 

iischsrae in Cubic Metres per Second 
YEAR JAN FEB MAR APR MAY JUN JUL AUG SEP OCT NOV DEC 



1974 


~ 


._ 


1.23 


3.83 


1.54 


0.59 


0.13 


0.22 


0.54 


0.68 


2.03 


0.38 


1975 


0.28 


0.21 


0.56 


4.03 


1.11 


0,43 


0.08 


0,06 


0.07 


0.33 


0.62 


1,13 


1976 


0.13 


0.13 


3.32 


2.68 


2,10 


0.18 


0.25 


0,09 


0.11 


0.19 


0.36 


0,22 


1977 


0.13 


0.11 


2.51 


3.27 


0.28 


0.13 


0.17 


0.38 


1.12 


0.81 


1.06 


0.69 


1978 


0.26 


0,17 


0.14 


4.19 


0,80 


0.14 


0.16 


0,55 


1.05 


1.20 


0.73 


0.42 


1979 


0.31 


0.30 


2.61 


4.98 


1.39 


0.39 


0.11 


0.29 


0.73 


0.99 


1.30 


0,97 


1990 


0.38 


0.11 


1.23 


4.65 


0.82 


0.59 


0,87 


0.58 


2.39 


2.49 


1,28 


0.57 


1981 


0.20 


2.05 


1.54 


2.36 


1.33 


0,94 


0,25 


0.16 


1.45 


1.63 


0.31 


0.47 



A-4 



TABLE A-7 COMMANDA CREEK NEAR COMHANOA 
SUtion 2DD01S 

Discharte in Cubic Metres per Second 



YEAR JAN FEB NAR APR MAY JUN JUL AUG SEP OCT NOV DEC 



1974 


— 


— 


— 


— 


-- 


— 


~ 


~ 


— 


— 


2.90 


1,35 


1975 


1.80 


1.32 


1.17 


5.26 


3.76 


1.06 


0,32 


0.11 


0.51 


0.76 


1.67 


4.74 


1976 


0.92 


0,82 


4.99 


6.44 


2.21 


0.83 


1.02 


0.30 


0.47 


0.52 


0.31 


0.94 


1977 


0.44 


0.45 


4.58 


5.86 


0.84 


0.35 


0.15 


0,32 


0.91 


1.96 


2.46 


1.50 


1978 


0.86 


0.54 


0.47 


4.74 


3.34 


0.83 


0.20 


0.42 


0.83 


1.79 


1.43 


1.37 



1979 0.78 0.56 3.82 7.34 2.71 1.12 0.44 0,65 0.32 1.78 2.37 1.19 

1980 1.02 0.43 1,83 7.72 1.63 2.19 1.24 0.99 1.43 2.92 2.17 1.08 

1981 0,64 2,63 2.54 6.31 2.56 2.00 0.67 0.44 1.78 2.59 1.26 1.27 



TABLE A-e NORTH HAGNETAUAN NEAR BURK'S FALLS 
Station 2EA005 (NET) 

DischsrSe in Cubic rtetres per Second 
YEAR JAN FEB NAR APR MAY JUN JUL AUG SEP OCT NOV DEC 



1970 


1.80 


1.52 


1,17 


8,49 


7,42 


2.93 


4.84 


2.07 


1.21 


3.19 


2.89 


3.44 


1971 


1.38 


0.95 


1.67 


11.40 


6,87 


2.42 


1.03 


0.98 


0.94 


1,19 


1.81 


3.35 


1972 


1.60 


0.63 


0.17 


5.74 


8.89 


2.02 


2.63 


2.23 


1.49 


1.56 


3.25 


2.02 


1973 


2.10 


1.80 


3.69 


7.99 


4.29 


2.12 


2,25 


2.11 


0.60 


1.47 


1.74 


2.87 


1974 


1.27 


1.53 


1.42 


8,80 


6.36 


2.78 


1.07 


0.70 


0.58 


2.67 


4.24 


2.74 


1975 


1.82 


1.33 


1.82 


6,35 


7.95 


1.53 


0.36 


0.55 


1.30 


2.70 


2.67 


4.08 


1976 


1.53 


1.36 


3.43 


11.89 


4.21 


1.49 


1.25 


0.25 


0.00 


0.37 


0.70 


1.85 


1977 


0.80 


0.48 


2.28 


10,76 


2.69 


0,41 


0.80 


0.14 


0.73 


5,05 


5.49 


4.05 


1978 


1.91 


1.09 


0.61 


4,64 


8.53 


1.94 


0.55 


0.46 


1.10 


2.69 


2.19 


2.10 


1979 


1.71 


1.22 


4.39 


10.00 


5.97 


2.19 


0,48 


1.14 


0.90 


2.23 


4.43 


4.29 


1980 


2.68 


1.32 


1.64 


13.20 


3.98 


3.20 


3.37 


3.03 


3.56 


3.87 


5.12 


3.49 



1981 1.32 2.52 5.01 11,10 4,57 2,70 1.98 1.76 6.70 6.05 2.34 2.23 



A-5 



TABLE A-9 NORTH MACHETAyAN RIMER ABOVE PICKERa LAKES 
Stsiion 2EA010 

DischsrSe in Cubic Metres per Second 



YEAR JAN FEB HAR APR HAY JUN JUL AUG SEP OCT NOV DEC 



1970 


0.72 


0.61 


0.82 


9.41 


3.98 


2,07 


3.98 


0.69 


1.68 


2.78 


2.39 


1,80 


1971 


1,14 


0.69 


1.00 


10.40 


4.13 


1.21 


0.57 


0.71 


0.67 


0.76 


1.60 


2.84 


1972 


1.07 


1.02 


1.26 


9.96 


9.31 


3.84 


1.90 


2.32 


1.25 


2.84 


2.83 


1,13 


1973 


1.19 


0.99 


5.51 


8.01 


3.82 


2.32 


1.41 


1,18 


1.03 


1.71 


2.76 


3,01 


1974 


2.32 


2.40 


2.79 


12.90 


5.94 


2.02 


1.27 


0.66 


0.90 


2.38 


4.82 


1.66 


1975 


1.76 


1.38 


1.88 


3.85 


5.55 


0.98 


0.94 


0.53 


2.52 


2.68 


2.79 


5.51 


1976 


1.11 


0,98 


9,37 


9.01 


4.63 


1.56 


1.43 


0.55 


0.99 


1.16 


1.63 


1.43 


1977 


0.74 


0.52 


6,11 


9.04 


1.56 


0.70 


0.61 


0.72 


2.34 


3.78 


4.19 


2.55 


1978 


1.03 


0.65 


0,62 


6,96 


5.87 


1.26 


0,51 


0.98 


1.43 


3.12 


2.05 


1.66 


1979 


1.04 


0.87 


6.81 


10,70 


3.83 


1.30 


0.88 


0.71 


0.92 


2.92 


4.63 


2.71 


1980 


2.79 


0,96 


2,12 


11,90 


2.63 


3.87 


3.08 


1.82 


3.14 


4,09 


3,70 


1,54 


1981 


1.02 


5.57 


3.14 


9.70 


3.96 


2.71 


2.67 


1.66 


5.70 


5.55 


2.26 


1.94 



TABLE A-10 SHAVANABA RIVER AT NO, 69 
Station 2EA012 

Hischsrie in Cubic Metres nt Sectmd 



YEAR JAN FEB iWR APR 



tMY 



JUN JUL AUG SEP OCT MOV DEC 



1974 


~ 


~ 


— 


20,40 


7.77 


3.10 


0.85 


0.09 


0.09 


2.07 


9.03 


3.53 


1975 


4.54 


2.98 


4.37 


14.30 


6.47 


1.00 


0.16 


0.05 


2.79 


2.20 


4.81 


9.61 


1976 


2.75 


3,02 


12.30 


13.90 


3.91 


1.26 


0.46 


0.08 


0.04 


0,05 


0.67 


1.94 


1977 


1.44 


1.23 


4.07 


11.60 


1.58 


0.30 


0.04 


0,02 


0.61 


4,06 


8.42 


6.70 


1978 


2.56 


1.55 


1.29 


16.40 


5.07 


1.38 


— 


— 


~ 


~ 


— 


~ 



A-6 



TABLE A-U ROSSEAU RIVER AT NO. U 
Station 2EB103 

Discharte in Cubic Metres per Second 



YEAR JAN FEB NAR APR HAY JUN JUL AUG SEP OCT «0V DEC 

1970 „ .. _ 11,70 2.32 0.53 4.30 i.55 0.74 3.17 3.33 3.11 

1971 - ~ - 

1972 - - - 

1973 ~ 

1974 

1975 - - - 

1976 _ _ - 

1977 „ „ „ 6,96 0,75 0,25 0,42 - 2.95 4.74 6.43 



19.10 


3.33 


0.44 


0.12 


0,11 


0.13 


0.16 


0.87 


4.18 


15.60 


4.30 


1.32 


0.41 


0.98 


0.40 


2,98 


2,56 


1,56 


9.51 


3.29 


1.90 


0.48 


2.21 


0.33 


2.04 


4.43 


2.47 


16.30 


3.91 


0.72 


0.22 


0.42 


0.87 


3.99 


5.31 


1.18 


15.20 


4.27 


0.44 


O.IS 


0.11 


0.46 


0.86 


1.52 


5.09 


— 


— 


~ 


0.44 


0.11 


0.06 


0.22 


1.54 


— 



TABLE A-12 BOB CREEK NEAR KINDEN 
SUtion 2HF004 

Discharge in Cubi Kttres per Second 
YEAR JAN FEB HAR APR HAY JUN JUL AUG SEP OCT NOV DEC 



1975 


— 


— 


— 


— 


— 


— 


— 


~ 


0.02 


0.02 


0.08 


0.37 


1976 


0.14 


0.15 


1.32 


0.88 


0.12 


0.02 


0.02 


0.01 


0,04 


0,04 


0,09 


0.07 


1977 


0.03 


0.04 


0.86 


0.38 


0,02 


0.00 


0.00 


0.00 


0,01 


0.14 


0,38 


0.24 


1978 


0.13 


0.05 


0.08 


0.96 


0,24 


0.01 


O.OO 


0.01 


0,00 


0,09 


0,22 


0.15 


1979 


0.12 


0.10 


0,95 


0.73 


0,19 


0.03 


0.01 


0.01 


0,00 


0.22 


0.45 


0.40 


1980 


0.13 


0.01 


0.46 


0,99 


0.07 


0.20 


0.32 


0.04 


0,05 


0,36 


0.28 


0.20 



1981 0.05 1.15 0.39 0.39 0.20 0.05 0,01 0,02 0,18 0.16 0.17 0.11 



A-7 



TABLE A-13 INBI/W RIVER NEAR PEMBROKE 
SUtion 2KC014 



Dischdrse in Cubic lletres ^er Second 
YEAR JAN FEB NAR APR «AY JUN JUL AUG SEP OCT MOV DEC 



1970 


0.81 


1.11 


1.36 


17.10 


10.70 


3.32 


2.82 


0.69 


0.70 


1.44 


2,36 


2.45 


1971 


l,U 


1.38 


2.15 


22.50 


10.40 


3.33 


1.24 


0.65 


0.46 


0.39 


0,95 


2.20 


1772 


1.00 


0.79 


0.87 


15.00 


17.50 


9.48 


4.66 


4.41 


2.63 


5.98 


11,90 


5.41 


1973 


4.06 


3.02 


10.90 


29.70 


10.40 


9.32 


3.60 


1.74 


1.31 


1.64 


2,76 


3,35 


1974 


2.08 


2.22 


4.02 


25.80 


19.80 


5.43 


2,07 


1.39 


0.88 


1.16 


4.19 


3.23 


1975 


2.64 


1.78 


3.36 


17.90 


8.77 


2.45 


0.85 


0.47 


0.51 


0.74 


0.93 


1.81 


1976 


1.28 


1.52 


11.20 


22.40 


9.04 


3.88 


5.42 


1.70 


1.79 


1.95 


2.96 


1.87 


1977 


1.43 


1.20 


8.65 


13.90 


4.02 


1.24 


0.56 


0.45 


0.67 


2.00 


4.24 


5.47 


1978 


3.46 


2.34 


1.82 


23.00 


12.50 


3.11 


0,85 


0.85 


0.64 


0.90 


1.14 


1,14 


1979 


1.27 


1.20 


5.54 


20.70 


11.60 


3.92 


1.39 


0.80 


1.17 


2.28 


3,35 


3.50 


1980 


1.74 


0.85 


4.86 


16.00 


7.45 


4.02 


5.72 


3.24 


1.81 


5.04 


5,25 


4.02 



1981 2.01 10.70 11.20 11.10 7,29 3.20 1.47 1,02 4.24 4.30 3,52 3,27 



I 



APPENDIX B 
SUMMARY OF CLIMATOLOGICAL DATA 



H 



TABLE B-1 SUHHARY OF METEOROLOGICAL DATA FOR BURK'S FALLS 

Precipitation in rtilliBetres and TeBperature in Decrees Celciys 





Year 


Parateter 


Jan 


Feb 


Mar 


Apr 


Hay 


Jun 


Jgl 


Aug 


Sep 


Oct 


Nov 


Dec 




1970 


Precipitation 


63.5 


_ 


70,9 


64,8 


67.3 


85.1 


181.4 


59.9 


108.2 


108.5 


69,9 


96.8 






Hean Hax Dailij Teip 


-8.3 


— 


-- 


10.3 


16.1 


22.2 


24.3 


23.9 


17.9 


13.5 


4.2 


-5.9 






Mean Hin Dai 1^4 TeiP 


-22.1 


— 


— 


-3.1 


4.1 


7.1 


12.2 


9.4 


7.1 


2,3 


-3.4 


-17.7 






Disx 


0.0 


~ 


— 


















lS.d 






Dtin 


0.0 


— 


~ 


















2.8 




1971 


Precipitation 


133.1 


146.8 


92.7 


40.9 


57.9 


75,2 


72.6 


100.8 


57.9 


49.8 


100.3 


108.0 






hean Hax Baib Tetp 


-B.9 


-3.9 


-0.6 


7.3 


18.0 


23.3 


22.6 


21.9 


20.7 


16.6 


2.5 


-1.2 






Mean Hin Daily Teip 


-21.3 


-15.6 


-14.9 


-4,6 


0.9 


8,4 


9.2 


8.7 


7.2 


5.2 


-4.8 


-12.1 






Diax 


2.8 


13.9 


42.0 


















41,7 






Diin 


0.0 


0,0 


0.0 


















1,7 




1972 


Precipitation 


1«.8 


112,8 


78.5 


63.2 


68.1 


112.8 


184.2 


163.1 


66.3 


90.9 


65.3 


190,8 






Hean Hax Daily TetP 


-3.9 


-6.6 


-2,1 


6.5 


19.0 


20,4 


23.9 


21,4 


18.3 


8.7 


1.6 


-3,2 






Mean Min Dai la Teip 


-12.7 


-21.1 


-15.2 


-6,1 


2,4 


6.4 


10.1 


9.8 


6.2 


-1.3 


-4.9 


-14,9 






Diax 


17.8 


2,8 


41.1 


















20.6 






Dain 


0,0 


0.0 


0.0 


















0.0 




1973 


Precipitation 


108.0 


53.8 


101.3 


76,5 


120.9 


148.6 


133.9 


118.6 


77.5 


100,1 


85.3 


96,0 






Mean Max Daily TetP 


-3.9 


-5.1 


5.6 


10.2 


15.4 


22.8 


24,8 


24.3 


18.2 


14.6 


3,0 


-3.8 






Mean Min Daily TetP 


-15.4 


-20.9 


-1.3 


-2,2 


1,7 


9.2 


10,9 


11,9 


4.9 


2.4 


-4.2 


-12,8 






Dsax 


45.0 


8.3 


174.5 


















37.8 






Dtin 


4.4 


0.0 


17.8 


















0.6 




1974 


Precipitation 


137.9 


98.8 


60.7 


94.5 


107.2 


102.6 


101.9 


127.8 


136.9 


126.2 


112.8 


51.0 






Mean Max Daily TeiP 


-4.8 


-6.7 


-0,2 


10.2 


14.1 


21.9 


23.6 


23,8 


16.4 


9.7 


3.9 


-0.5 






Mean Min Bails Teiip 


-16.4 


-21.4 


-11,5 


-3.0 


1,2 


7.3 


9.4 


9.0 


3.1 


-2.3 


-5.2 


-9.6 






Dtax 


32.2 


0.0 


51.7 


















27.2 






B»in 


0.0 


0.0 


2.B 


















0.0 




1975 


Precipitation 


174.5 


123.7 


63.0 


85,9 


26.4 


69.6 


103.4 


119.6 


142.5 


62.5 


109.2 


73.4 






Mean Max Daily Ten? 


-3.6 


-3.4 


-0,3 


5,8 


21.5 


23.7 


25.9 


24.2 


16.2 


12.9 


7.6 


-4.3 






Hean Min Daily leap 


-15.7 


-15.7 


-14,3 


-8.3 


4.4 


7.8 


9.9 


8.4 


3.6 


1.2 


-1.6 


-15.8 






Dmx 


25.0 


23.9 


52,8 


















35.0 






Dtin 


0.0 


0.0 


0.0 


















0.0 




1976 


Precipitation 


78.7 


114.8 


152.9 





121.2 


90.4 


74,2 


62.0 


84.3 


71,9 


113.5 


109.0 






Hean Max Daily TeiP 


-8.4 


-1.2 


2.6 


~ 


15.6 


24.8 


23,4 


23.8 


17.8 


8.9 


0.6 


-7.3 






(tean Min Daily Teip 


-22.4 


-14.8 


-10.0 


~ 


3.2 


9.6 


10.7 


7.9 


5.3 


-1.7 


-7.3 


-21.1 






Dtax 


10.0 


38.3 


134,5 


















1.7 






Dain 


0.0 


0.0 


22.2 


















0.0 




1977 


Precipitation 


80.0 


99.3 


92.7 


89.5 


33.0 


97.1 


86,2 


119,8 


163.6 


34.5 


168.4 


30.0 






Mean Max Daily TetP 


-10.8 


-3.4 


5.5 


12.3 


21.3 


21.1 


24.9 


21.2 


16.5 


11.3 


4,7 


-3.6 






ftean Min Daily fetp 


-23.2 


-13.5 


-6.9 


-3.3 


3,6 


6.9 


11.3 


9.7 


7,5 


-0.2 


-3.9 


-12,4 






Daax 


0.0 


11.2 


188.3 


















16.5 






Dain 


0.0 


0.0 


13.3 


















0.0 




1978 


Precipitation 


91.4 


30.2 


66.3 


74,4 


55.8 


105.0 


85.0 


161.4 


138.0 


87.4 


77.0 


110.4 






Mean Max Daily TetP 


-8.2 


-7.0 


-0,4 


6.1 


20.0 


21.6 


23.9 


23.5 


16,8 


10,2 


3.7 


-3.1 






Mean Min Daily Teap 


-19.8 


-20.8 


-14.7 


-4.2 


5.2 


7,3 


9.8 


10.0 


4.7 


-0.1 


-5.8 


-14.4 






Daax 


4.0 


0.0 


44.5 


















8.5 






toin 


0.0 


0.0 


0.0 


















0*0 




1979 


Precipitation 


188.0 


74.0 


13.2 


— 


— 


71,8 


96,0 


88.2 


100.6 


141.4 


118.2 


108.3 






Mean Max Daily Teip 


-B.8 


-5.5 


~ 


~ 


17.8 


23.0 


26.0 


21.6 


18.6 


11.5 


5.0 


-1.5 






Mtan Min Daily Teap 


-17,5 


-18.2 


— 


~ 


4.6 


8.2 


10,5 


9,8 


6.5 


3.3 


-2,0 


-12,7 






Omx 

MM 


3.0 


15.5 


— 


















27,5 






0.0 


0.0 


-— 


















0.0 




1980 


Precipitation 


133.7 


60.4 


113.3 


101,8 


60.0 


157,2 


92,4 


124.2 


96.8 


123.0 


92.6 


116.8 






Mean Max Daily Teap 


-5.1 


-6.7 


0,3 


10,1 


13.8 


19.0 


23.3 


24.2 


16.8 


8.5 


2.3 


-6.6 






Mean Min Daily Teip 


-16.4 


-18.5 


-8.6 


-0,9 


3.7 


5.7 


11.6 


13.0 


6.1 


-0.4 


-5.9 


-20.7 






Oaax 


22.5 


1.0 


98.5 


















9.5 






Alio 


0.5 


0.0 


2.0 


















M 




1981 


Precipitation 


80.4 


107.4 


67.2 


96.2 


81.4 


109.8 


94.0 


137,2 


165.2 


109.2 


49.6 


118.1 






Mean Max Daily Tetp 


-8.4 


0.4 


2.2 


10.9 


16.3 


21,5 


25.4 


23.5 


16.4 


9.3 


4.7 


-2.7 






Mean Min Daily Teap 


-23.2 


-9.7 


-8.6 


-1,5 


2.4 


8.7 


10.7 


10.0 


6.5 


-1.1 


-6.0 


-11.3 






Bmx 


3.5 


94.7 


114.2 


















9.5 






Dain 


0.0 


17,5 


14.0 


















0.0 

































M 



Year 



TAM.E B-2 SUfiHARY OF «£TE(KOL06ICAL DATA FOR OWLK RIVER 

Preci^iUtion in Hilliietres and Tet^erature in Decrees Celcius 
Paraieter Jan Feb Har Apr May Jun Jul Aua Sep Oct Kov 



Dec 



I 



1970 


Precipitation 


22.1 


42.4 


57.7 


46.7 


88,4 


62,0 


112.3 


24.4 


79.8 


75,2 


50.5 


65,0 




Hean Hax Daily Te«p 


-10.0 


-6.3 


1.0 


11.2 


17.2 


23,2 


25.2 


24.8 


18.3 


13.4 


4.4 


-7.1 




Hean nin Daih Tetp 


-21,2 


-18.6 


-10.8 


-1,1 


4.9 


9,6 


13,9 


U.9 


8.9 


3.7 


-2,1 


-15.4 




Ohx 


1.7 


20.6 


62,2 


















IM 




Diin 


0.0 


0.0 


0.0 


















M 


1971 


Precipitation 


54.6 


121.2 


46.7 


22.4 


43.7 


89,4 


43.4 


61.5 


34.5 


29.5 


66,3 


94.5 




Mean Hax Dai la Te»p 


-9.8 


-4.4 


0.7 


8,3 


19.6 


24.4 


23.9 


23.8 


21,7 


16.2 


2,4 


-3.6 




Hean Hin DailH TeiP 


-20.6 


-15.2 


-11.2 


-2.6 


4.9 


10.2 


11,6 


11.0 


9.6 


5.7 


-3.6 


-13.7 




Dmx 


1.1 


18.9 


54.5 


















37.8 




Diin 


0.0 


0.0 


0.0 


















0.0 


1972 


Precipitation 


3S.6 


82.3 


81.8 


54.6 


36.1 


142.5 


138,4 


123.2 


89.2 


81.8 


72.9 


126.2 




Hean Hax HsUh lev 


-4.0 


-6.7 


-1.6 


7.7 


20.2 


22.1 


25,2 


22.5 


18,6 


8.9 


2.0 


-5.3 




Hean Nin Dai la Teip 


-16.9 


-19.4 


-12.6 


-3.2 


6.1 


10.2 


13.8 


12.6 


3.4 


1.2 


-0.4 


-14.8 




Omx 


33.9 


5.0 


41.7 


















3.3 




Mn 


0.0 


0.0 


0,6 


















0.0 


1973 


Precipitation 


30.7 


43.9 


91.4 


64.8 


74.7 


239.0 


73,9 


63.8 


41.9 


74.9 


64.5 


107.2 




Hean Hax Dailtj leap 


-4.6 


-5.5 


5,4 


10,7 


16.8 


23,6 


27.0 


25.9 


18.8 


14.6 


2.6 


-3.8 




Mean Nin Dai la Teip 


-14.2 


-18.9 


-2.3 


-0.3 


5.0 


12.2 


13.8 


15.4 


8.7 


5.3 


-1.6 


-11.2 




teax 


21.1 


6.1 


168.5 


















33.9 




Da in 


8.9 


0.0 


12.8 


















3.9 


1974 


Precipitation 


62.7 


59.4 


50.3 


92,5 


125.0 


62.2 


85,6 


79.5 


53.8 


51.8 


108.5 


51.1 




Mean Hax Dai la Teip 


-5.7 


-5.8 


0.3 


11.1 


14.9 


22.8 


25,3 


24.7 


17.3 


10,1 


3.5 


-0.7 




Hean Kin Daila Tetp 


-15.6 


-IB.l 


-10.5 


-1.2 


3,6 


11,2 


13,6 


12.9 


6.4 


0,6 


-2.4 


-7.2 




Dtsx 


27.8 


9.4 


62.2 


















26,1 




Diin 


0.6 


0.0 


3,3 


















0,6 


1975 


Precipitation 


63.2 


42.9 


60.5 


44.2 


37.8 


77.7 


62.7 


28.4 


77.2 


93.7 


55.9 


49,3 




Hean Hax Dai la Teip 


-5.1 


-3.6 


0.7 


7.9 


23.8 


24.9 


27.8 


26.3 


17.7 


13.1 


7.1 


-5.6 




Hean Hin Dai la Te^p 


-14.4 


-13.8 


-10.4 


-4.1 


8,1 


12.4 


14.5 


13.3 


7,9 


3.3 


-1.3 


-15.8 




Dmx 


19,4 


18.4 


63.3 


















26.7 




Diin 


1.6 


0.0 


0.6 


















0.0 


1976 


Precipitation 


55.6 


67.8 


116.1 


9,9 


111.3 


76.2 


51.6 


34.3 


64,0 


47.8 


47.0 


39.4 




Mean Max Dai la Tetp 


-8.9 


-1.2 


1.9 


13.7 


16.4 


26.2 


24.7 


24.1 


18.2 


9.7 


0.8 


-7,5 




Mem Hin Dai la Tew> 


-14.3 


-13.7 


-9.9 


0,1 


4.8 


13.0 


12.9 


11.9 


7.6 


0.1 


-5.6 


-19,1 




D»3X 


4.4 


32.2 


127.2 


















0.6 




Diin 


0.0 


0.0 


7.8 


















0.0 


1977 


Precipitation 


45.0 


45.5 


61.3 


83.0 


25.6 


63.4 


58.5 


107.9 


91.3 


50.6 


103.6 


69,0 




Hean Max Dai la Tetp 


-9.1 


-3.8 


6.7 


12.9 


22,7 


21.9 


26.8 


23.2 


17.5 


11.5 


5.1 


-5.0 




Hean Hin Dai la leap 


-19.0 


-12,6 


-4.4 


-1.0 


7,0 


10.3 


13.7 


12.8 


9.5 


1.7 


-1.6 


-12,4 




Okax 


0.0 


13.4 


216.1 


















16.2 




Diin 


0.0 


0.0 


20,5 


















1.1 


1978 


Precipitation 


ao.o 


3.9 


50.9 


30.2 


46,5 


31.8 


43.4 


135.8 


72.9 


47.7 


70,8 


110.8 




Hean Max Dai la Teap 


-8,1 


-6,1 


0.7 


7.5 


21,4 


22.9 


25.6 


24.4 


17.1 


10.8 


5.1 


-3.2 




Hean Hin Daila Tetp 


17.7 


-18.6 


-11.4 


-2.6 


8.2 


10.1 


13.5 


12,4 


6.S 


2.0 


-4.9 


-12.0 




Daax 


7.2 


0.0 


65.7 


















21,0 




Dain 


1.1 


0.0 


0.0 


















0,0 


1979 


Precipitation 


75.3 


22.3 


33.6 


86.9 


94,4 


54.2 


93.3 


55,9 


86.0 


158.0 


74.3 


42.8 




Mean Max Daila Te»p 


-8.3 


-8.7 


4.9 


10.0 


18.3 


23.8 


26.5 


22,5 


19,4 


10.4 


5.3 


-2,4 




Hean Hin Daila Tetp 


-17.2 


-20.9 


-5.7 


-0.5 


6.0 


10,3 


14.1 


11,6 


7.2 


2.9 


-1,8 


-10.6 




Daax 


9.5 


23.5 


166.9 


















33.0 




tain 


0.0 


0.0 


15.0 


















0.0 


1980 


Precipitation 


33.4 


22.0 


70.7 


86.4 


71.3 


131,2 


158.3 


72,6 


110.7 


116.4 


57.5 


53.2 




Hean Max Daila Tea? 


-5.7 


-5.7 


1,8 


11,1 


19.5 


20,6 


25.5 


24.7 


17.3 


9.0 


1.6 


-8.9 




Hean Hin Daila Teap 


-15.7 


-17.1 


-9.6 


0,6 


5.2 


8,1 


13.7 


14.7 


6.9 


1.2 


-5.0 


-20.8 




Oaax 


13.0 


4.0 


107.8 


















3.5 




Siii) 


0.0 


0.0 


2.5 


















0»0 


1981 


Precipitation 


18.7 


74.7 


66,6 


76,3 


79.5 


80.8 


43.1 


129.1 


146.0 


84,3 


43.3 


42.0 




Mean Max Daila Teap 


-10.2 


1.1 


3.3 


11,7 


17,7 


23.2 


26.5 


24,4 


16,8 


9.5 


4.9 






Hean Hin Daila Teap 


-21.8 


-10.4 


-6.2 


O.S 


5.3 


11.4 


13.5 


12.5 


9.1 


-1. 5 


-3.2 


-8.6 




Oaax 


7,5 


105.5 


131.8 


















19.5 




Diin 


0.0 


16.5 


13,0 


















0.5 



B-3 



Year 



TABLE B-3 SUmftRY OF METEOROLOGICAL DATft FOR HUNTSVILLE 

Precipitation in Millisetres snd Teapersture in Decrees Celcius 



P3r3i»eter 



Jan Feb ttsr Aft Mas Jun Jul 



Au3 



3e? Oct Nov Dec 



1970 


Precipitation 




40.7 


41.9 


60.2 


68,3 


106,4 


91.2 


202.7 


65.3 


'7.5 


101.1 


71.1 


111.8 




Mean Ma;; DsiIh Tsiup 


-7.8 


-4.1 


1.0 


10.7 


16.4 


22.9 


24.9 


23.9 


13.0 


13.8 


5.2 


-4.6 




Mean Min Caila 


Tesp 


-19,0 


-16.1 


-11.7 


-i.9 


6.0 


10.7 


14.6 


13.5 


9.9 


5,0 


-0.6 


-14,1 




Oma;? 




0.0 


25.0 


61.0 


















22,8 




Dmin 




0.0 


0,0 


0.6 


















2.8 


1971 


Precipitation 




99.3 


133.7 


63.3 


35.6 


45.0 


75.2 


67.1 


93.5 


31.2 


55,1 


35.9 


114.3 




Mean Max Baily Temp 


-7.0 


-2.7 


0.9 


0,3 


IS. 2 


23.6 


23.1 


22.7 


20.9 


16.9 


3.7 


-0.1 




Mean Min Dai is 


TeRP 


-13.8 


-14.6 


-12.4 


-2.2 


3.3 


11.6 


12,8 


12,0 


11,4 


7,5 


-1.8 


9.4 




Dtax 




5.6 


21.1 


27.2 


















53.3 




Dnin 




0.0 


0.0 


0.0 


















7,2 


1972 


Precipitation 




135.4 


39.4 


64.3 


64.0 


83.3 


80.0 


125.7 


118,9 


79.5 


109.5 


64.3 


128.3 




Mean Max Daily Te»p 


-2.9 


-5.2 


-1.1 


7.2 


19.3 


20,4 


23.9 


22.2 


18,9 


9.4 


2.3 


-2.1 




Mean Min Daily Tetp 


-14,8 


-17.9 


-12.1 


-3.4 


5,9 


9.8 


13.3 


13.0 


9.2 


1,1 


-2.4 


-10.9 




Dnax 




35.0 


6,1 


52,3 


















33.9 




Diin 




0.0 


0,0 


1.1 


















0,6 


1973 


Precipitation 




119,4 


53,3 


76.7 


58.9 


102,6 


235.5 


71,6 


46.2 


54.6 


83,1 


93.0 


90.7 




Mean Max Daila Te»p 


-3.3 


-3.7 


5.9 


10.6 


15,5 


22.9 


25,2 


24.9 


18.9 


14.7 


3.9 


-2,7 




Mean Min Daily 


Teup 


-12.7 


-17.1 


-1.2 


-0.1 


4.9 


12.4 


14.6 


15.7 


9.0 


5.3 


-1.5 


-10,3 




Dnax 




36.1 


117.0 


173.5 


















41.1 




Diin 




5.0 


0,0 


32.8 


















6.7 


1974 


Precipitation 




148.8 


52.6 


55.1 


100.1 


95.8 


50.0 


48.8 


117.6 


129.5 


104.1 


104.4 


43.9 




Mean Max Daily 


Tesf 


-3.9 


-5.1 


0.3 


11.0 


13.9 


21.9 


23.9 


23.9 


16.9 


10.4 


4.4 


0.0 




Mean Min Daily 


TeiP 


-14,1 


-18.1 


-8.8 


-0.7 


4.0 


11.6 


13,2 


13.3 


7.4 


0.6 


-1.5 


-4.6 




Diax 




33.3 


10.6 


50.6 


















31.1 




Dniin 




0,0 


0.0 


2.2 


















0.0 


1975 


Precipitation 




114.0 


91.7 


65.3 


56.1 


57.2 


52.8 


87.1 


93.2 


105.4 


54,6 


100.3 


105.7 




Mean, Max Daily 


Tesip 


-2.7 


-2,5 


0.3 


6.4 


21.4 


23.6 


25.6 


24.3 


16.8 


13.6 


3.6 


-3.4 




Mean Min Dai y 


Tern? 


-11.9 


-12.6 


-11.3 


-4.3 


3.2 


13.6 


14.2 


13.1 


7,7 


4.7 


-0.1 


-12.9 




Dnait 




26.1 


25.0 


60.6 


















37.2 




Diiin 




2.8 


1,7 


2.2 


















0.0 


197A 


Precipitation 




92.2 


102,6 


114.0 


23.1 


53,4 


33,3 


42,2 


66,6 


152.4 


53.3 


120.7 


82.8 




Mean Max Daily 


TeidP 


-6.9 


-0,4 


3.2 


11.7 


15,5 


25,0 


23,9 


23,7 


18.0 


10,1 


1.3 


-6.3 




Mean Min Daily 


TeiP 


-20.1 


-13.3 


-7,3 


0.2 


5,0 


12.9 


13,7 


12.2 


7,7 


2.3 


-4.8 


-18.6 




Daax 




10.4 


54.4 


166.7 


















4.1 




Dnin 




0.0 


0.0 


22.2 


















0.0 


1977 


Precipitation 




115.0 


48.9 


59.4 


55.6 


23,6 


61.8 


70,9 


115.8 


118,0 


33.6 


147.0 


110.2 




Mean Max Daily 


Tesp 


-9.5 


-3.2 


5.8 


12.8 


21.3 


21,7 


25.2 


22.6 


17.7 


11.8 


5.5 


-4.0 




Mean Min Daily Temp 


-20.0 


-12,6 


-5.4 


0.6 


7,0 


10,3 


14,3 


12.3 


9.9 


1.3 


-1.4 


-11.3 




Hmy, 




0.0 


9.6 


185.4 


















11,0 




Dtin 




0.0 


0,0 


5.6 


















0.0 


1978 


Precipitation 




56.9 


13,6 


34.7 


51.8 


41.6 


88,4 


87.0 


145.2 


110.8 


61.6 


82.1 


33.0 




Mean Max Daily 


Teup 


-7,6 


-6,7 


0.1 


6.8 


20.0 


22.1 


24.4 


24.3 


16.3* 


9,5* 


3.4* 


-2,7* 




Mean Min Daily TeiP 


-13.2 


-19.6 


-13.6 


-3.1 


6.9 


9.9 


13.2 


13.2 


5.3* 


0.3* 


-5.2* 


-3.7* 




Diax 




7.5 


0,0 


29.7 


















21.0* 




Di»in 




0.0 


0,0 


0.0 


















0.0* 


1979 


Precipitation 




161.0 


44.4 


91. 2t 


28.0 


116.0 


52.4 


85.4 


94.5 


90.4 


138,3 


88.3 


55. 4 




Mean Max Daily 


Temp 


-7.1 


-7.6 


4.8 


9.4 


17.4 


22,' 


25,9 


22,1 


19.3 


10,6 


6.1 


-0,3 




Mean Min Daily 


leap 


-15.0 


-19,6 


-5.5 


-0,3 


6.3 


9.4* 


13.7 


12.3 


3.7 


3.3 


-0.5 


-3.6 




Diiiax 




4.4 


13.8 


167,1 


















30,0 




Do in 




0.0 


0.0 


9.4 


















1,1 


1980 


Precipitation 




64.9 


48,7 


74,1 


86.7 


60.7 


137.1 


116,6 


127.7 


55.3 


119.3 


40.6 


143.0 




Mean ilax Dai lb- 


TeBP 


-3.6 


-5.6 


1.3 


10.7 


19.1 


19.2 


24.9 


24.7 


17.6 


8.3 


2.9 


-6,0 




Mean Min Daily 


Tea? 


-13.6 


-16.9 


-5,0* 


0.3 


6.0 


7.8 


13.3 


15.3 


3.4 


1,6 


-4,0 


-17,6 




Dnax 




17.2 


1.7 


72,5 


















6.4 




Dain 




0.6 


0.0 


-15,0* 


















0.3 


1981 


Precipitation 




63.2 


64.3 


67.3 


99,0 


64.4* 


39.7 


67,2 


200.7 


99.6* 


48.6* 


37.21 


43.0 




Mean Max Daily 


Teap 


-8,4 


0.7 


2.6 


11.1 


15.0* 


21.8 


25,5 


23.9 


15.6* 


3.4* 


3,8* 


-2.3 




Mean Min Daily 


Teap 


-22.3 


-9,1 


-7.3 


0,6 


3.2* 


11.0 


14.0 


13.2 


7.6* 


0,4* 


-3,5* 


-3.7 




Dma;; 




5,3 


54,9 


68.0 


















6.4 




Dnin 




0.0 


3.9 


4.4 


















0,3 



* Data frot Dwiaht 



M 

TARE B-4 SUGARY OF HETEORaOGICAL iATA FOR MINDEN 

Precipitation in ^1ilUietre5 and leapersture in Beirees Celcius 

Year Paraieter Jan Feb (far Apr May Jun Jul Aua Sep Oct Nov Dec 

1?75 Precipitation 97.8 120,9 47.1 74.9 74.2 91.9 76.7 88.4 106.4 50.0 91.4 96.0 

Sean «3X DailH TeiP -2.5 -2,2 0.7 7.2 22.4 24.6 26.7 25.1 17.3 14.4 9,1 -3.1 

Mean Hin Dails Te»p -12.2 -13.2 -12.1 -6.1 6.1 10.5 11.4 10,2 5.9 2.2 -1.4 -13.3 

D«3x 23.9 28,3 57.8 38.9 

Dtin 2.2 1.6 5.0 0.0 

1976 PreciPiUtion 98.0 134.1 140.7 22.6 80.3 64.5 83.8 110.7 105.7 59.9 63.0 73.7 
Hean Max Daila Teip -6.7 0,2 3,1 13,6 16,5 25,4 24.4 24.2 19.2 9.9 2.1 -5.6 
Mean Kin Daily TeiP -20.7 -12.4 -9,2 -1.2 2.5 10.2 U.l 9.7 5,7 -0.8 -6.3 -18.9 
Dmx 8.9 67.8 142.2 10.0 
Diin 0.0 0,0 9.4 0,0 

1977 PreciPiUtion 81.1 73.4 86,8 33,4 16.0 86.3 40.1 132.7 89.2 95.0 138.5 116.3 
Mean Max Daily Teip -8.9 -4,0 6.0 13.4 22.7 22.8 25.9 22.5 18.1 12.3 6.0 -3,2 
Mean Hin Daily Tea? -20.3 -14.0 -6.0 -0.9 3.4 7.6 11.8 10.4 9,2 0.3 -2.2 -12,2 
Diax 0.0 15.5 201.0 19.6 
Dtin 0.0 0.0 10.1 2.8 

1978 Precipitation 103.4 lO.l 72.2 58.7 81.3 39,5 98.2 126.5 78.3 23.8 77.0 102.3 
Mean Max Daily TeiP -6.9 -5.1 1.1 7.6 20.7 22,5 25.3 25.0 18.7 11.7 5.7 -1,6 
Mean Min Daily Tew -18,0 -20.1 -13,2 -4,0 5.7 7.8 11.0 U.l 5.6 -0.5 -4.8 -U.l 
Dtax 7.7 0.0 70.2 26.5 
Biin 0.0 0.0 0.0 0.0 

1979 Precipitation 84.3 51.4 56.5 99.6 110.2 104.9 78.2 88.8 41.4 156.8 156.7 98.4 
Mean Max Dails TetP -6.3 -6.9 5.5 10.2 17.9 23.3 26.4 22.8 20.3 11.5 6.4 -0,5 
Mean Min Daily TeiP -15.2 -21.7 -5.5 -1.3 5.2 8.8 11.3 10.7 6.3 2.5 -1.6 -9.3 
taax 7.0 26,5 198,0 48.5 
Diin 0.0 0.0 16.0 2.0 

1980 PreciPiUtion 84.6 26.0 109.4 106.2 68.7 102.6 145,6 97,1 71.5 137,4 74,2 142,8 
Mean Max Daily Te»p -3,0 -4,9 2.0 n,4 20.2 19.7 25,1 25.2 18.7 9.9 3.2 -5,1 
Mean Min Daily Tetp -13.3 -17.7 -10.1 -0.5 4.1 5.6 11.6 12.9 6.2 0.1 -5.1 -19.0 
Dmx 37.0 4.5 109.0 13.0 
D>in 0.0 0,0 3.0 0,0 

1981 Precipitation 32.0 112.8 49.7 83.0 85.2 103.1 54.2 91.7 116.7 80.9 61.3 67.0 
Mean Max Daily TeiP -8.0 1.3 4.6 11,0 17.7 22.8 25.6 24.5 17.3 10.2 6,0 -1.3 
Mean Min Daily Teip -22,8 -9,8 -6.9 -0.9 2.9 9,1 11.0 10.3 6.9 -0.9 -5.4 -10.2 
Dmx 10.0 91.9 163.6 18.0 
Oiin 0.0 15.5 28.0 0.0 



B-S 



Year 



Paraieter 



TABLE B-5 SUMMARY OF NETEOfiaOGICAL DATA FOR NORTH BAY 

Precipiistion in Milliutres and Teti^ersture in BHrees Celcius 
Jan Feb Har Apr Hay Jun Jul Aua Sep Oct Nov 



Dec 



I 



1970 


Precipitation 


39.1 


~ 


3A.8 


61.0 


167.4 


125.7 


227.8 


135.9 


154.7 


49.9 


39.1 


82.3 




Mean Max Daily Tetp 


-9,7 


— 


0,4 


10.6 


15.8 


22.9 


25.1 


24.7 


13,2 


13,3 


3.1 


-7.1 




Mean Min Daily TeiP 


-21.1 


— 


-11.1 


-1.3 


4.9 


10.5 


14.9 


13.6 


9.1 


4.3 


-3.8 


~ 




Dftax 


1.7 


0.0 


58.0 


















13.9 




Dtin 


0.0 


0.0 


0.0 


















1.7 


1?71 


Precipitation 


54.9 


99.6 


63.8 


37.6 


53.1 


32.3 


58.4 


38.4 


89.9 


64.3 


72.4 


91.7 




Mean Max Daily Teip 


-9.8 


-4.1 


0.3 


8.5 


18.2 


24,1 


24.1 


22.7 


21.1 


15.7 


2.7 


-1.5 




Mean Min Daily TetP 


-20,4 


-15,2 


-11.7 


-2.0 


4.2 


11.7 


12.2 


11.7 


10.8 


7.1 


-3.0 


-11.1 




Dtax 


1.7 


19.4 


56.1 


















35.4 




Diin 


0.0 


0.0 


1.7 


















3.9 


1972 


Precipitation 


78.5 


87.1 


65.3 


51.3 


49.0 


74.7 


104,9 


210.8 


40.9 


49.1 


45.0 


125.2 




Mean Max Daily Tetp 


-3.9 


-7,1 


-2.0 


7,3 


19.8 


21.3 


24.2 


21.7 


18.2 


8.9 


2.1 


-5.4 




Mean Min Daily Teip 


UA 


-19.6 


-12.8 


-4,1 


6.3 


9.7 


13.2 


12.0 


8.4 


0.6 


-4.1 


-14.8 




Otax 


28.9 


11.7 


42.8 


















6.1 




Diin 


0.0 


0.0 


0.0 


















0.0 


1973 


Precipitation 


Si.l 


31.5 


65.3 


50.5 


109.2 


143.3 


97.8 


93.7 


58.2 


90.4 


47.3 


74.9 




Mean Max Daily Teup 


-4.2 


-6.1 


5.8 


10.6 


16.3 


23,1 


25.8 


24.3 


14.9 


14.8 


2.3 


-4.5 




Mean Min Daily Tetp 


-13.3 


-19.1 


-2.4 


-0.5 


4.3 


12.7 


14.4 


15.2 


4.8 


5.3 


-3.3 


-12.1 




Daax 


30.0 


5.0 


185.9 


















27.8 




Diin 


5.0 


0.0 


16.7 


















1.7 


1974 


Precipitation 


79.2 


63.5 


31.5 


114.6 


104.6 


74.4 


44.5 


85.6 


100.4 


124.0 


48.8 


55.4 




Mean Max Daily Teip 


-7.2 


-7.1 


0.0 


9.4 


14.1 


22.5 


24.1 


24.1 


14.3 


9.3 


3.3 


-0.5 




Mean Min Bails Tetp 


-16.4 


-\9.6 


-10.3 


-1.7 


4.3 


11.9 


13.9 


12.7 


4.4 


0.2 


-2.2 


-7.6 




Dtax 


22.2 


7.2 


50.6 


















29,4 




Dtin 


0.0 


0.0 


1.7 


















0.4 


1975 


Precipitation 


117.9 


83.8 


58.2 


53.6 


74.7 


44.7 


141.5 


25.1 


84.8 


89.2 


131.3 


41.0 




Mean Max Daily TetP 


-5.3 


-3.9 


-0.1 


6.7 


21.7 


24.4 


26.8 


24.4 


14.4 


13.0 


4.7 


-5.2 




Mean Min Daily Teip 


-16.1 


-13.9 


-11,6 


-4.9 


9.0 


12,9 


15.2 


11.3 


7.9 


3.9 


-1.2 


-14.4 




Btax 


12.2 


24.4 


56.1 


















29.4 




Dtin 


0.0 


0.0 


0.0 


















0.0 


1976 


Precipitation 


68.6 


111.3 


95.8 


12.4 


116.6 


76.2 


80.8 


48.5 


99.4 


50.3 


88.9 


56.4 




Mean Max Daily Tetp 


-8.8 


-1.7 


-0.1 


11.2 


16.1 


25.1 


23.8 


23.4 


17.7 


9.2 


0.8 


-9.1 




Mean Min Daily TetP 


-21.2 


-15.1 


-11,2 


-0.1 


5.2 


13.5 


13.7 


11.9 


7,2 


0.8 


-5.4 


-21.3 




Otax 


3.3 


52.8 


103,9 


















0.6 




Dtin 


0.0 


0.0 


8,9 


















0.0 


1977 


Precipitation 


74.9 


76.9 


74,9 


52.4 


40,3 


44.8 


142.1 


141.7 


140,2 


45.8 


95.2 


75.3 




Mean Max Daily Tetp 


-10. 1 


-5.0 


5.9 


12.4 


20.1 


20.0 


24.8 


21.4 


15,4 


11.2 


5.0 


-4.4 




Mean Min Daily TetP 


-19.5 


-13.8 


-4,9 


-0.8 


6.6 


9.2 


14.0 


12.3 


8.1 


1,8 


-2.0 


-10.9 




Dtax 


0.0 


14.0 


192.3 


















12.9 




Diin 


0.0 


0.0 


18,5 


















0.0 


1978 


Precipitation 


6^.2 


16.2 


3i.3 


45.9 


45.5 


45.9 


109,6 


146.1 


134.4 


79.9 


61,2 


106.9 




Mean Max Daily TetP 


-9.8 


-6.9 


-0.1 


7.0 


21.3 


21.7 


24,4 


22.7 


14,5 


10.1 


3,6 


-4.3 




Mean Min Daily Tetp 


-18.6 


-19.0 


-13.2 


-3.0 


8.5 


10,3 


13.5 


12.3 


8.0 


2,3 


-4.4 


-13.4 




Dtax 


1.1 


0.0 


55.1 


















11.5 




Diin 


0.0 


0,0 


0.0 


















0.0 


1979 


Precipitation 


87.2 


37.5 


39.6 


100.7 


30.9 


78.4 


43.4 


118.4 


71,4 


117.2 


94.1 


67.4 




Mean Max Daily Tetp 


-9.2 


-10,1 


3.1 


8.3 


16.5 


21.4 


25.7 


21.2 


17,8 


9.5 


5.1 


-1.9 




Mean Min Daily Tetp 


-17.5 


-20.8 


-7.6 


-1.0 


5.5 


10.1 


14.9 


13.0 


9,1 


3,7 


-1.1 


-10.2 




Dtax 


5.2 


8.0 


130.3 


















32.5 




Dtin 


0.0 


0.0 


2,9 


















2.0 


1980 


Precipitation 


69.0 


47.0 


121.4 


92.0 


39.6 


138.4 


134.8 


85.8 


173,8 


105.5 


80.2 


58.4 




Mean Max Daily Tetp 


-5.7 


-7,0 


1.0 


9.2 


18,8 


17.7 


23.9 


23.8 


15.4 


7.0 


2.2 


-9.2 




Mean Min Daily Tetp 


-14.6 


-16.9 


-9.3 


-0.3 


6.4 


7.3 


14.1 


15.0 


5.6 


-0,7 


-4.1 


-20.7 




Dtax 


20.5 


2.0 


36.0 


















3.8 




Dtin 


0.0 


0.0 


3.0 


















«,0 


1981 


Precipitation 


33.4 


75.4 


77.2 


100.7 


74.8 


126.4 


31,6 


132.1 


145.9 


112.8 


40.8 


43.4 




Mean Max Daily TeiP 


-9.4 


-1.2 


0.6 


9.5 


14.7 


22.0 


25.2 


24.1 


14.5 


9,4 


5.0 


-2.9 




Mean Min Daily TetP 


-20,3 


-10.4 


-8.0 


-1.4 


5.0 


12.2 


13.7 


13,5 


9.1 


0.9 


-3.4 


-9.3 




Dtax 


8.5 


59.1 


77.6 


















'?:8 




Dtin 


0.0 


10,7 


e.7 



















I 



Year 



TABLE B-6 SUHMARY OF METEOROLOGICAL DATA FOR SUDBURY 

Precipitation in MilliMtres and Tetwrature in Decrees Celcius 
Paraaeter Jan Feb Nsr Apr )1au Jun Jul Aud Sep Oct Nov 



M 



Dec 



1970 


Precipitation 


37.1 


31.0 


40.6 


56.1 


132.6 


94.7 


146.1 


70.6 


257.3 


67,1 


66.8 


57.9 




Hean Hax Dail^i Teip 


-10,7 


-7.6 


-1.0 


9,9 


14.2 


22.7 


24.8 


24,3 


18,0 


12.8 


3,1 


-6.6 




Hean Kin Dail» Tetp 


-21,0 


-20,4 


-12.2 


-1.4 


3.7 


U.4 


15.4 


12.6 


9.9 


3,6 


-4.2 


-15.0 




^ 


0.0 


3.9 


35,6 


















10.& 




0.0 


0.0 


0.0 


















0.0 


1971 


Precipitation 


67.1 


90.2 


50,3 


32.0 


62,5 


50.5 


95.0 


55.4 


69,1 


61.5 


71.6 


93,2 




Hean Kax Itailv TeiP 


-10.8 


-4.6 


-0,1 


6.2 


17.7 


24.2 


23.7 


22.3 


21.2 


14,2 


1.7 


-3.7 




t^vn Kin Dail<i Tetp 


-21.2 


-15.4 


-11.8 


-3.6 


4.7 


11.4 


13,1 


10.9 


11.4 


6.4 


-4.8 


-13.1 




tSi 


0.0 
0.0 


18.9 
0.0 


42.2 
1.1 


















26.7 
1.1 


1972 


Precipitation 


72.6 


79.2 


65.8 


60.7 


45.7 


111.5 


78.2 


115.6 


102.1 


36,3 


31.2 


37.6 




Mean Hax Dail<j Tetp 


-5.3 


-7.3 


-3.9 


4,9 


19.4 


20.8 


24.7 


20,9 


16.9 


7,5 


0.3 


-6.9 




Hean Min Dail;* Tetp 


-18.4 


-19.1 


-14.7 


-4,7 


6.4 


8.8 


13.2 


11.6 


6.5 


-0.2 


-4,8 


-16.6 






15.0 


l.l 


16.1 


















M 




Mn 


0.0 


0.0 


0.0 


















m 


1973 


Precipitation 


32.8 


29.7 


78.2 


41.7 


114.3 


83,8 


99.3 


100.8 


27.4 


89.9 


61.0 


49,3 




Hean Kax Daily Teip 


-4.0 


-8.1 


3.9 


9.1 


14.7 


22.2 


24.4 


25.6 


18.7 


13.2 


1.1 


-5.9 




Hean Hin Daily Teap 


-13.6 


-19.3 


-3.9 


-1.7 


4,4 


12.9 


14.1 


16.4 


8.9 


4.7 


-4,4 


-13,8 




fiiax 


31.1 


0.0 


129.8 


















9.4 




Diin 


1.7 


0.0 


17,8 


















0.0 


1974 


Precipitation 


83.1 


76.2 


31.5 


71.1 


85.9 


8S.6 


78.7 


64,8 


104.1 


106.2 


53,8 


67.6 




Hean Hax Dail^i Teap 


-8.1 


-9.2 


-2.4 


7.8 


12.8 


20.9 


24.2 


23,2 


15.1 


8.0 


2,5 


-2.1 




He») Hin Dails Teap 


-16.9 


-20.1 


-12.1 


-2.6 


3.2 


10.9 


13.4 


12,5 


5.1 


-0.3 


-3.6 


-9.0 




Diax 


17.8 


1.1 


27.2 


















12.2 




lain 


0.0 


0.0 


0.0 


















0.0 


1975 


Precipitation 


116.8 


53.1 


65,3 


33.5 


71.9 


55,1 


27,2 


21.8 


75.9 


35.1 


109.7 


43.9 




Hean Hax Daily Teip 


-5,9 


-5.3 


-1,9 


7,1 


21.0 


23.2 


26.7 


24.9 


15.7 


11.4 


6.4 


-5.7 




Hean Kin Dail<< Teap 


-16.7 


-14.2 


-12,8 


-4.8 


9,1 


12.4 


14.4 


13.2 


6.5 


2.7 


-1,4 


-14,3 




Daax 


8.9 


13.9 


26.1 


















26.1 




Dain 


0.0 


0.0 


0.6 


















0,0 


1976 


Precipitation 


52.8 


95.6 


101,3 


17.5 


95.3 


44.2 


48.8 


75.2 


158.8 


57.4 


44,5 


29.2 




Hean Hax Daily Teap 


-9.4 


-3.4 


0.4 


10,1 


15.4 


25.2 


24.2 


24.2 


17.4 


7,3 


-0.8 


-9,3 




Hean Hin Daily Teap 


-22.3 


15.6 


-10.7 


-1.1 


3.9 


12.7 


13.1 


11.8 


6,1 


-0.4 


-7.3 


-21.3 




Oaax 


0,6 


17,3 


19,5 


















0.0 




Oain 


0.0 


0.0 


5.0 


















0.0 


1977 


Precipitation 


50.9 


64.5 


71.2 


36,2 


22.1 


47.5 


156.7 


101.8 


125.9 


58.4 


96.9 


86.3 




Hean Hax Daily Teap 


-11.8 


-4.9 


4,9 


10,8 


21.6 


21.9 


25.1 


20.8 


16.7 


10.5 


3.2 


-6.7 




Hean Kin Daily Teap 


-21.4 


-13.0 


-5.5 


-1.6 


7.2 


10,5 


13,6 


10.5 


3.4 


1.0 


-3.4 


-13.3 




Oaax 


0.0 


10.0 


158.6 


















6,1 




Dain 


0.0 


0.0 


8.4 


















0,0 


1978 


Precipitation 


83.6 


14.3 


42.4 


51,1 


70.7 


72,2 


92,3 


157.8 


160.7 


63.8 


60.6 


100,5 




Hean Max Daily Teip 


-U.l 


-8.0 


-1,4 


5.6 


20,2 


20,8 


23,9 


22.8 


15.6 


8.3 


2,2 


-6.1 




Hean Hin Daily Teap 


-19.5 


-17.8 


-13.2 


-4.1 


7.3 


8.4 


12,6 


12.0 


6,3 


0.4 


-6,5 


-14.9 




Daax 


0.0 


0.0 


35.8 


















1.4 




Dain 


0.0 


0.0 


0.0 


















0.0 


1979 


Precipitation 


68.0 


36.0 


105,1 


176.8 


48.7 


79.1 


80.9 


97.5 


54.0 


96.7 


100,1 


66.1 




He.an Hax Daily Tea? 


-10.8 


-11.5 


1.3 


6.3 


15.9 


21.8 


25.3 


20.7 


18.0 


8.3 


2,7 


-2.4 




Hean Hin Daily Teap 


-18.7 


-22.7 


-8.9 


-2,1 


5.1 


10.1 


13.7 


11.2 


6.7 


1.6 


-4,0 


-11.0 




Daax 


0.1 


10.3 


99,9 


















32.3 




Dain 


0.0 


0.0 


1.7 


















0.3 


1980 


Precipitation 


96.1 


30.7 


78.7 


111.4 


61.9 


64,4 


68,0 


100.0 


148,1 


67,0 


68,0 


56,5 




Hean Hax Daily Teap 


-7.2 


-8.2 


-1.4 


8.2 


18.7 


18,9 


24,7 


24.0 


15.6 


6.8 


1,0 


-8,3 




Hean Hin Daily Teap 


-16.1 


-18.4 


-11.3 


-0.7 


6.1 


7.8 


13.1 


14.1 


6.1 


-1.2 


-6.2 


-20.5 




Daax 


9.8 


1.1 


36.4 


















0,7 




Dun 


0.0 


0.0 


0,9 


















0.0 


1981 


Precipitation 


16.0 


83.6 


37.3 


130,5 


57.0 


112.2 


30.9 


68.7 


177.1 


109,4 


36.9 


41.2 




Hean Hax Daily Teap 


-10.6 


-2.2 


0.8 


9,5 


15.8 


20.5 


26.1 


24.0 


15.3 


3.0 


3.8 


-4.3 




Hean Hin Daily Teap 


-20.5 


-11.5 


-7.7 


-1.6 


4.0 


10.0 


13,6 


13.2 


7.0 


-0,1 


-4.1 


-U.l 




Daax 


0.0 


47.4 


79.4 


















7.3 




toin 


0.0 


0,9 


5.5 


















0.0 
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APPENDIX C 



DETAILS OF DISCHARGE CALCULATIONS 
FOR THE ROUTING MODEL 



The subbasin discharge must be defined as belonging to one of 
the following classes 

- no regulation 

- regulation subject to a rating curve 

- regulation subject to rule curve operation (or historic 
lake levels ) . 

The subbasin continuity is schematized in Figure C-1 for each 
of these three categories. The following subsections define 
these regulating classes in terms of subbasin continuity. 

CI - NO REGULATION 

Where no regulation exists, the discharge from the subbasin 
is simply the sum of subbasin runoff, upstream subbasin 
discharge and external inflow. The subbasin discharge is 
written as 

Qijk = ^ijX + SQij^ + lij^ (CI) 



where 



Q. ., - discharge (m /s) from subbasin k 

R,- ^v ~ subbasin runoff (m /s) from subbasin k 
1 I jc 



3 



SQ . ., - summation of upstream subbasin discharges (m /s) 
which flow directly into subbasin k 




REGULATING LAKE 
LEVEL Y,j ;Y,(j + ,)k 



SUB-BASIN WITH NO REGULATION 
LEGEND 
(^^ SUB- BASIN 
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'ijk 



- external inflow {in /s) to subbasin k 

i - subscript denoting year 

j - subscript denoting month 

k - subscript denoting subbasin. 



The external inflow (lijx) is user supplied and repre- 
sents inflow into the drainage basin. If concurrent drainage 
basins are modeled separately for one watershed, then the 
discharge from the upstream drainage basin becomes the exter- 
nal inflow to the subbasin immediately downstreeim. 



C2 - REGULATION SUBJECT TO 
A RATING CURVE 



The discharge may be regulated by a natural control or by a 
man-made structure at the outlet of the regulating lake. The 
discharge is represented by a rating curve of the form 



Q ■ ■, = a 
i]k 



V + Y 

ijk i(j+l)k 



(C2) 



where 



Y . 
ijk 



i(j+l)k 
a, b 



elevation of the regulating lake (with respect to 
the elevation of the weir crest) at the start of 
month j 

elevation of the regulating lake at the end of 
month j 

coefficients. 



The subbasin continuity expression is written as 



^ijk = ^^ijk - Yi(j+l)k - ^jk ^ PRijx^(AREGI^)(f) + 



R. ., + SQ. ., + I . .. 

ijk ijk i]k 



(C3) 



I 
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where 



e.^ - evaporation coefficient (m/month) for subbasin k 

PR. ., - actual precipitation (m/month) which fails on the 



ijk 



regulating lake in subbasin k 



2 

AREGL. - area of regulating lake (km ) for subbasin K 

f - conversion factor (to convert MCM/ month to m /s). 



To calculate the subbasin discharge (Qijk)/ Equations C2 
and C3 are solved simultaneously. Because of the nonlinear 
relationship of discharge to lake elevation, the chain-of- 
lakes model sequentially calculates the subbasin discharges 
for ten equal time segments within the monthly period. The 
model outputs only the average discharge over the month and 
the month-end lake elevation. 



C3 - REGULATION SUBJECT TO 
RULE CURVE OPERATION 
(OR HISTORIC LAKE LEVELS) 



To describe the imposed regulation, the user can supply 
either 

- rule curves which define the elevation of the regulating 
lake as a function of the calendar date. The chain-of- 
lakes model adjusts the user-supplied rule curve so that it 
corresponds to month-end calendar dates. 

- historic lake levels for each month and year of the 
simulation period. 

Any arbitrary datum may be used in the definition of the rule 
curves or the historic lake levels. 
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Given the above imposed regulation, the continuity expression 
is given as 



QijX - ^^ijk - Yi(j+l)k - «jk ^ PRij^)(AREGL^)(f) + 



(C4) 



^ijk + SQij^ + liji. 



The elevations of the regulating lake CY. .^, Y-i(-i+i)v3 ^"^ 

Equation C4 are defined by the rule curves or by historic 
lake levels. For operative purposes, the model treats 
historic lake levels as if they were rule curve values. 

If the discharge (Qi-;v)» aa calculated from Equation C4, 
is less than zero, then the lake level CY., . ..,] cannot 
be achieved, therefore 



i]k 



and 



x{3 + l)k i:k 3k i:k ^ (AREGL)^(f) ^ ^"^^^ 
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APPENDIX D 



DETAILS OF CALIBRATION PROCEDURES 



The following two calibration techniques are available in the 
chain-of- lakes model 

(i) Historic replication of a specific basin discharge 
record. 

(ii) Improvement of the runoff prediction equations for a 
specific data record. 

These calibration techniques are described in the following 
subsections. 



Dl - HISTORIC REPLICATION 

When calibrating so that the historic basin discharges are 
replicated, the model adjusts the subbasin runoff on a 
month- by- month basis so that the calculated basin discharge 
is equal to measured basin discharge. This calibration 
method is event specific and, therefore, is limited to only 
those months where historic discharge data is available. 
Given the requirement to replicate the historic basin 
discharge, the model distributes the basin-wide calibration 
into the individual subbasins by prorating the subbasin 
runoff, as calculated from the runoff prediction equations. 

To replicate the historic basin discharge, the following 
steps are followed to determine the subbasin runoff. 
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1 - On a month- by -mo nth basis, compute the runoff for each 
subbasin by using the runoff prediction equations. 

The subbasin runoff is represented as R. , , where the 

13 k 

superscript n denotes the iteration (in this case n=0) 
and i = year, j = month, k = subbasin. 



2 - Route the subbasin runoff {R. ., ) throughout the drainage 

network. Determine the calculated basin discharge {Q . . ) 

I jm 

where m represents the furthest downstream subbasin. 

3 - Increment the iteration counter (n) by 1. Given the 

measured basin discharge (QM. .}, the total measured basin 

runoff is approximated by backrouting the measured basin 
discharge with respect to net changes in lake storage 
throughout the basin. The total measured basin runoff is 
calculated as 



RM 



ij ij I ijm 1^1^ ijk J 



For all, i = 1, 2,...NyRS and j = 1, 2,... 12 months 

where N = the total number of subbasins contoured in 
the drainage network. 

NYRS - the total number of years in the simulation 
period 

The quantity in brackets in the above equation is the net 
change in lake storage throughout the basin. 

4 - Compute the ratio: 

RATIO = RM. 

U 



^ R 
^ijk 



k=l 



I 
I 
I 
I 

I 
f 
I 

I 

I 
I 
f 
t 
I 
I 
I 
I 
I 
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5 - Use the ratio calculated in Step 4 to adjust the subbasin 
runoff 



R.-l = (RATIO) (K..0) 



For all, k = 1, 2,...N subbasins 



6 - Route the subbasin runoff (R. ., ) throughout the drainage 
network. Determine the calculated basin discharge 



7 - If additional iteration is necessary, repeat Steps 3 to 6 
inclusive. 



The calculated basin discharge (Q. •„) will implicitly 

converge toward a value equal to the measured basin discharge 

{QM. .). More than one iteration is required only in those 

^ J 

cases in which the total volume of stored water varies from 
iteration to iteration. 



D2 - IMPROVEMENT TO RUNOFF 
PREDICTION EQUATIONS 



Instead of the month by month replication of historic in- 
flows, the improvement to the runoff prediction equations has 
a broader basis for application. With a given measured 
record for basin discharge, revised regression coefficients 
are derived to make the runoff prediction equations basin 
specific. Extrapolation of the data can then be made by 
utilizing these revised regression coefficients in subsequent 
runs of the chain-of- lakes model. 



D-4 



To determine the revised coefficients, this calibration tech- 
nique is performed as a posts imul at ion exercise. If the 
simulation period consists of (NYRS) calendar years, then for 
each month there are (NYRS) pairs of calculated and measured 
basin discharges. The measured basin discharges form the 
basis for calculation of the revised regression coefficients. 
Linear regression analysis is used to estimate these coeffi- 
cients for each month. It is important to note that the 
revised regression coefficients are only applicable for the 
derivation of total basin runoff. The basin-wide calibration 
must then be distributed to the subbasin level. Working on a 
month- by- month basis, the distribution used herein is similar 
to that used in the historic replication technique. 

Before a detailed description of the calibration technique is 
given, it is necessary to review the general form of the 
runoff prediction equations. 

The general form of the runoff prediction equations are as 
follows. 



April to November 






where 



3 
R, .. - subbasin runoff (m /s) 

A_ - area (km ) of the subbasin k 



(Dl) 



k 



\ 



2 
- area (km ) of lakes in subbasin k 
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2 

A - area {km ) of swamps in subbasin k 

P. . - monthly precipitation (m) computed with respect 
•^ to runoff prediction 

a, - a^ - coefficients (vary from month to month) 

i - subscript denoting year 
j, - subscript denoting month 
le - subscript denoting subbasin 



December to March 



^ijk ^ ^1 -^ ^2 V "■ ^3 V "" ^a^^^^ "" ^5°ijk^ (D2) 



where 

'^ijk " degree days ("C) above freezing 
Sjuation Dl can be simplified to 



\jk - ^k -^ ^Ik^jk ^^^^ 



and Equation D2 can be simplified to 



13k 2k '^2k i]k 2k x]k 13k 



Equations Dl, D2, D3 and D4 are similar to Equations 12, 
13, 14 and 15 of Section 3.4 

Note that the above equations are linear with respect to 
variables Pj^,, and P^ ^v ^-i -;v • ^Y extension, the total 

XjK ^J^ IjK 

basin runoff - R. ., is also a function of P and (P D). 
E 1 jk 

k=l 
Let the furthest downstream subbasin be denoted by the 
subscript m. Then the total basin runoff can be 
represented as either 



I 
I 

I 

I 
I 
t 
I 

I 
1 
I 
I 
I 
I 
I 
f 
I 
I 
I 
I 
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E^ijk " ^1 "^ ^l^ijm (D5) 

k=l 



or 



N _ 

Z^ijk " -^2 ^ ^2^ijm ^ ^2^ijm°ijm (D6) 

k=l 



for the summer and winter months respectively. 



If the parameters P. . and P. . D. . are applicable to the 
*^ ijm 13m ijm ^^ 

entire basin then 



^1 Z Ik /. ^2 Z 2k 

k=l k=l 



k=l k=l 



C = ^ Y 

k=l 



A basic premise of the chain-of-lakes model is the limited 
data base. In general. Equations D5 and D6 adequately 
represent the total basin runoff. Given the total measured 
basin runoff (RM. .), then the regression equations 

RM = A' + B'P (D7) 

ij 1 1 ijm 



and 



RM. , = a; + b:p. . + c;p. . d. . (ds) 

13 2 213m 2i3mi]m 

can be derived. Coefficients A,', Ai, B ' , BI, and Cj, make 

12 12 2 

the model basin specific. 
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The following steps outline the methodology of this 
calibration technique, 

1 - For the entire simulation period of (NYRS) calendar 

years, compute the runoff for each subbasin by using 
the runoff prediction equations. 

The subbasin runoff is represented as R.^,,r where the 

1 J K 

superscript n refers to the iteration number (in this 
case n = 0) . 

2 - For the entire simulation period, route the subbasin 

runoff (R. ") throughout the drainage network. E)eter- 

1 J It 

mine the calculated basin discharge (Q; 4™) ' where m 
represents the furthest downstream subbasin. 

3 - Increment the iteration counter (n) by 1. Given the 

measured basin discharge (QM. .), the total measured 

basin runoff is approximated by backrouting the 
measured basin discharge with respect to changes in 
lake storage and for net evaporation throughout the 
basin. The total basin runoff is calculated as 



RM 



ij ^ 13 \ ijni I ijk / 
V k-1 ' 



4 - For each of the 12 months, determine by linear regres- 

sion the coefficients A| , A^, E'^, B^ and C^ of 

Equations D7 and D8 . 

5 - On a month-by-month basis, calculate the estimate of 

total basin runoff (RM..) by using the regression 

^ J 

equations as derived in Step 4. 
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6 - On a month- by-month basis, compute the ratio 



RM. . 
RATIO = ij 

N _ O 

k=l 



7 - On a month- by -mo nth basis, use the ratio (calculated in 
Step 6) to adjust the subbasin runoff 



■*iJK- (RATIO) (R..°) 



For all, k = 1, 2,...N subbasins 



8 - Route the subbasin runoff (R. ., ) throughout the drain- 
age network. Determine the calculated basin discharge 

1 jm 



9 - If additional iteration is necessary, repeat Steps 3 to 
8. 

The derivation of the revised regression coefficients 
(Step 4) is a posts iraulation analysis. Subsequent steps in 
the methodology are performed on a monthly basis. Using 
the revised regression coefficients, an estimate of basin 
runoff is determined. The estimated basin runoff is then 
compared with the calculated basin runoff (iteration n = 
0). The ratio so calculated is then used to distribute 
runoff throughout the subbasins. At this stage, this 
method is similar to that used in the historic replication 
technique. 
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The monthly average for total basin discharge, as computed 
in this calibration technique, will implicitly converge 
toward the corresponding average for measured discharge. 
This convergence is inherent because the linear regression 
analysis will replicate the mean of the measured runoff. 
Basin discharge is equal to the sum of basin runoff and the 
change in volume of stored water. More than one iteration 
is required only for those cases in which the total volume 
of stored water varies from iteration to iteration. 

The results of the calibration (A', Al, B ' , B^, and CL) can 

be read into the program and used in subsequent data 
extension runs. For application of the calibration results 
to data extension, omit Steps 3, 4 and 9 from the above 
methodology . 
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FILE SPECIFICATIONS AND 
DATA INOUT REQUIREMENTS 



The data input to the chain-of-lakes model is described 
herein. Note that all data input is unformatted. 



El - CHAIN-OF-LAKES 

MODEL: FILE STRUCTURE 



Within the chain-of-lakes computer program, the following 
unit numbers are accessed. 



Unit Number 



File Type 



Description 



I 
I 



5 
6 
8 



10 
11 

12 
13 
14 



15 



input 

output 

input 

input 

input 
input 
input 
input 
output 



output 



input control data 

output listing 

temperature data for runoff 
calculations 

precipitation data for runoff 
calculations 

monthly precipitation data 

measured discharge data 

measured month-end lake levels 

external inflow data 

basin discharge data; can be 
used in the chain-of-lakes 
model as external inflow data 
for the consideration of 
downstream basins 

subbasin discharge data, which 
replaces the data output of 
the 'HYDRO' module of the 
Serpent River model 
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E2 - INPUT CONTROL DATA 



1 - TITLE 1 



- enter 3 lines of title description 



2 - TITLE 2 

3 - TITLE 3 

4 - NYRl, MTHl, NYR2 , MTH2 

where NYRl - start year of simulation 

MTHl - initial calendar month for simulation 

(e.g., for January start month, set 

MTHl = 1) 
NYR2 - final year of simulation 
MTH2 - final calendar month for simulation (e.g. 

for December termination, set 

MTH2 =12). 



5 - IPRTl, IPRT2, IPRT3 

where IPRTl - print flag to obtain the monthly summary 

output; can be used to check subbasin 
continuity 
= do not print monthly summaries 
= 1 print monthly summary output 
IPRT2 - print flag to obtain the subbasin 
discharge summary output 
= do not print subbasin discharge 

summaries 
= 1 print subbasin discharge summary 
output 
IPRT3 = print flag to obtain the month-end lake 
level summary output 
= do not print month-end lake level 

summaries 
= 1 print month-end lake level summaries 
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6 - lADJl, IADJ2 

where lADJl - flag to indicate the calibration such that 

there is a month- by-month historical 
replication of basin discharges (see 
Section 5.1) 
= do not adjust data 
>_ adjust data, lADJI represents the 

number of iterations to be performed 
IADJ2 - flag to indicate the calibration tech- 
niques to improve the runoff prediction 
equations (see Section 5.2) 
= do not adjust data 
>_ adjust data, IADJ2 represents the 

number of iterations to be performed 

NOTE: lADJl and IADJ2 cannot both be 
greater than zero. 

7 - NBASIN 

where NBASIN - the total number of subbasins contained 

in the watershed 

REPEAT CARD GROUPS 8-18 INCLUSIVE FOR (NBASIN) SUBBASINS 

8 - NAME 

where NAME - subbasin name or identifier (_<8 

characters) 

9 - NODE, NREG, NINE, NDIS 

where NODE - the node number used to represent this 

subbasin 
NREG - flag to indicate the type of discharge 
regulation for this basin (refer to 
Section 4.2.2) 
= no regulation 



E-4 



= 1 regulation by means of a control point 
at the outlet of the regulating lake 
(discharge calculation follows a weir 
discharge equation (Q = ayt>) 
= 2 imposed regulation by means of rule 
curves or measured lake levels 
NINF - flag to indicate whether external inflow 
data (see E8) is to be input for this 
subbasin 
= do not input data 

= 1 input external inflow data for this 
subbasin 
NDIS - flag to indicate whether measured 

discharge data (See E6 ) is to be input 
for this subbasin 
= do not input data 

- 1 input measured discharge data for this 
subbasin 

10 - NBDWN, NTIES, CNBUP{J), J = 1,NTIES] 

where NBDWN - the node number representing the subbasin 

which is immediately downstream from the 
current subbasin 
= if there is no downstreeim basin (i.e., 
the discharge from this subbasin is 
the outlet from the watershed 
NOTE: only 1 watershed should be 

considered in each run of the 
chain-of- lakes model 
NTIES - the total number of subbasins which 

directly discharge into the current 
subbasin 
= if there is no upstream subbasin 
NBUP - the node number representing the subbasins 
which discharge directly into the 
current subbasin 
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11 - AB, AL, AS 

where AB - area of subbasin (km^) for use in the 

runoff prediction equations 

AL - area of lake within subbasin {km2) for 

use in the runoff prediction equations 

AS - area of swamps within subbasin {km2) for 
use in the runoff prediction equations 

SKIP CARD GROUPS 12-17 INCLUSIVE IF NREG = 
INPUT CARD GROUP 12 IF NREG = 1 
SKIP CARD GROUP 12 IF NREG = 2 

12 - AREGL, Yd), A, B 

where AREGL - area of regulating lake {kni2) foj. this 

subbasin 
y(l) - initial lake level (m) for this regulating 
lake with respect to weir crest 
elevation 
A, B - coefficients for discharge calculation, 
the discharge Q(I) in month 1 is 
calculated as 

B 



Q(I) = A 



Y(I) + Yd + 1 ) 



(m^/s) 



where Yd) and YCl+l) are the initial 
and final regulating lake elvations in 
month I 

SKIP CARD GROUPS 13-17 INCLUSIVE IF NREG = 1 
INPUT CARD GROUPS 13-17 INCLUSIVE IF NREG = 2 

13 - AREGL 

where AREGL - area of regulating lake (km2) for this 

subbasin 
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14 - IRC 

where IRC - flag to indicate whether rule curves or 

measurd lake levels are read in to 
describe the imposed regulation for this 
subbasin 
= input rule curve data 
- 1 input measured lake levels (see E7 ) 

SKIP CARD GROUPS 15-17 INCLUSIVE IF IRC = 1 

15 - NRCPTS 

where NRCPTS - the number of data points which are used 

to describe the rule curve for the reg- 
ulating lake of the current subbasin 

16 - [RCDAYS (1), I = 1, NRCPTS] 

where RCDAYS - are the cumulative calendar days (1-365; 

Jan 1=1) which corresponds to the 
input rule curve elevations (RCCRV) 

17 - [RCCRV (I), 1=1, NRCPTS] 

where RCCRV - are the regulating lake elevations (m) 

which define the rule curve with respect 
to the calendar days (RCDAYS) 
NOTE ; The chain-of-lakes program 

interpolates rule curves for the 
month end calendar dates. 

18 - lEQN 

where lEQN - flag to indicate type of runoff prediction 

equation 
= historic sequential runoff equations 
= 1 ranked runoff equations 



I 
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REPEAT CARD GROUP 19 FOR (12) MONTHS 

19 - Ca(J), J = 1,6] 

where a - coefficients for monthly regression equa- 
tions which are used to describe basin 
runoff [called C(J) in Fortran code] 

- see Tables 2 and 3 

20 - NEVAP 

where NEVAP - the total number of evaporation records to 

be input 

REPEAT CARD GROUPS 21 and 22 FOR (NEVAP) RECORDS 

21 - NEV, [NEVP(I), 1=1, NEV] 

where NEV - the total number of subbasins for which 

this evaporation data is applicable 
NEVP - the node numbers of subbasins for which 
this evaporation data is used in the 
calculation of evaporation loss 

22 - [EVAP(J), J = 1, 12] 

where EVAP - monthly evaporation coefficients 

(m/ month) 

- to calculate the evaporation loss in a 

given month, the evaporation coefficient 
is multiplied by the area of the regu- 
lating lake for those subbasins which 
have regulation control 

- EVAP(l) corresponds to the January evapor- 

ation coefficient 
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23 - ICOEFF 

where ICOEFF - flag to indicate whether the linear 

. regression coefficients (A', A', B', 

B-, CI) are to be input by the user (as 

derived by the calibration technique to 
iit^rove the runoff prediction 
equations) 
- refer to Section D2 in Appendix D 

= do not input regression coefficients 
= 1 input regression coefficients 

INPUT CARD GROUPS 24, 2 5 and 26 IF ICOEFF ^ 1 

24 - [A2(I), I = 1,3; Al ( I ) , 1=4, 11; A2(12)] 

where Al and A2 are monthly regression coefficients 
input by the user to adjust the runoff 

- Al and A2 are equivalent to A' and A' in Equations D7 

and D8 

25 - CB2(I), I = 1, 3; Bl (I) = 4, 11; B2 (12)] 

where Bl and B2 are monthly regression coefficients 
input by the user to adjust the runoff equations 

- Bl and B2 are equivalent to B' and Bl in Equations D7 

and D8 

26 - [C2(I), 1=1, 3; 0, 0, 0, 0, 0, 0, 0, 0,; C2(12)] 

where C2 - monthly regression coefficients input by 

the user to adjust the runoff equations 
- C2 is equivalent to CI in Equation DB 
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E3 - TEMPERATURE DATA 

FOR RUNOFF CALCULATIONS 



The temperature data file is set up so that it is directly- 
compatible with the runoff regression equations. If the 
regression analysis has determined that runoff in a particu- 
lar month is dependent on temperature records for the current 
and previous months, then the input temperature data consists 
of the aggregate of the temperature records for those months. 

1 - NTEMP 

where NTEMP - the total number of separate temperature 

data records to be input 

REPEAT CARD GROUPS 2, 3 AND 4 for (NTEMP) RECORDS 

2 - NBT, [NBTMP(J), J = 1, NBT] 

where NBT - the total number of subbasins for which 

this temperature data set is applicable 
NBTMP - the node numbers of those subbasins which 
utilize this temperature data set in the 
regression equations of runoff 

3 - NYRS 

where NYRS - the total number of years of monthly 

temperature data to be input for this 
particular temperature record 

REPEAT CARD GROUP 4 FOR (NYRS) YEARS 
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4 - NYEAR, CX{J), J = 1, 12] 

where NYEAR - calendar year 

X(J) - monthly temperature data { "C degree days 
above freezing/month) which is used 
directly in the regression equations to 
determine subbasin runoff for a 
particular month 
- X(l) corresponds to January data 



E4 - PRECIPITATION DATA 

FOR RUNOFF CALCULATIONS 



The precipitation data file is set up so that it is directly 
conpatible with the runoff regression equations. If the 
regression analysis has determined that runoff in a particu- 
lar month is dependent on the precipitation in the current 
and previous months, then the input precipitation data 
consists of the aggregate of the precipitation values for 
those months. 

1 - NPPT 

where NPPT - the total number of separate precipitation 

records to be input 

REPEAT CARD GROUPS 2, 3 and 4 FOR (NPPT) RECORDS 

2 - NBPT, [NBPPT(J), J = 1, NBPT] 

where NBPT - the total number of subbasins for which 

this precipitation data set is 
applicable 
NBPPT - the node numbers of those subbasins which 
utilize this precipitation data set in 
the regression equations of runoff 



I 
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3 - NYRS 

where NYRS - the total number of years of monthly 

precipitation data to be input for this 
particular precipitation record 

REPEAT CARD GROUP 4 FOR (NYRS) YEARS 

4 - NYEAR, CX(J), J= 1, 12] 

where NYEAR - calendar year 

X(J) - monthly precipitation data (m/month) which 
is used in regression equations to deter- 
mine subbasin runoff for a particular 
month 
- X(l) corresponds to January data 



E5 - MONTHLY PRECIPITATION DATA 



The monthly precipitation data in this file is equivalent to 
the monthly historic record. 

1 - MPT 

where MPT - the total number of separate precipitation 

data records to be input 

REPEAT CARD GROUPS 2, 3 AND 4 FOR (MPT) RECORDS 

2 - MBPT, CMBPPT(J), J = 1, MBPT] 

where MBPT - the total number of subbasins for which 

this precipitation data set is 
applicable 
MBPPT - the node numbers of those subbasins which 
utilize this precipitation data in the 
water continuity calculations 



I 



E-12 



3 - NYRS 

where NYRS - the total number of years of monthly 

precipitation data to be input for this 
particular precipitation record 

REPEAT CARD GROUP 4 for (NYRS) YEARS 

4 - NYEAR, [X(J), J 1, 12] 

where NYEAR - calendar year 

X(J) - monthly precipitation data (m/month) which 
is used in the water continuity calcula- 
tions for the regulating lake in a 
particular subbasin 
- X{1) corresponds to January data 



E6 - MEASURED DISCHARGE DATA 



For any subbasin in the drainage network, measured discharge 
data can be input and the chain-of- lakes model will give the 
comparative statistics between calculated and measured dis- 
charges. Calibration, based on the measured discharge data, 
is only performed for the farthest downstream subbasin 
(defined by variable NBDWN = 0? See E2). The calibration 
technique is dependent on the value of parameters lADJl and 
IADJ2 (See E2 ) . 



REPEAT CARD GROUPS 1 and 2 FOR EACH SUBBASIN IN 
WHICH VARIABLE NDIS = 1 (See E2 ) 



1 - NC, NYRS 

where NC - the node number representing the subbasin 

for which the following historic dis- 
charge data is applicable 
NYRS - the total number of years of discharge data 
for this subbasin 
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REPEAT CARD GROUP 2 FOR (NYRS) YEARS 

2 - NYEAR, [X(J), J - 1, 12] 

where NYEAR - calendar year 

X(J) - monthly discharge (ni3/s) 

- X(l) corresponds to January data 



E7 - MEASURED MONTH-END 
LAKE LEVELS 



Within the chain-of- lakes model, the discharge from each 
subbasin may be described by the following 

- unregulated 

- regulation as defined by a discharge weir or other natural 

or man-made control structure 

- imposed regulation by either rule curves or the measured 

lake levels. 

This data set contains the measured lake levels for those 
subbasins with imposed regulation. 



REPEAT CARD GROUPS 1, 2 AND 3 FOR EACH SUBBASIN IN WHICH 
VARIABLES NREG = 2 AND IRC = 1 



I - NCODE, NYRS 

where NCODE - the node number representing the subbasin 

for which the following lake level data 
is applicable 
NYRS - the total number of years of lake level 
data to be input for this subbasin 

REPEAT CARD GROUPS 2 AND 3 FOR (NYRS) YEARS 
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2 - NYEAR, EL VI 

where NYEAR - calendar year 

ELVI - the January 1 lake elevation (m) for the 
regulated lake for this subbasin 

3 - CX(J), J = 1, 12] 

where X(J) - month end lake elevations (in) for the 

regulated lake for this subbasin 
- X(l) corresponds to the January month-end 
lake elevations 



E8 - EXTERNAL INFLOW DATA 



An external inflow can be applied to any subbasin. This data 
could represent the channel inflow to a subbasin which is 
gauged or has previously been calculated by the chain-of- 
lakes model. 



REPEAT CARD GROUPS 1 ' AND 2 FOR EACH SUBBASIN IF VARIALE 
NINF = 1 {See E2 ) 



1 - NC, NYRS 

where NC - the node number representing the subbasin 

into which this external inflow is input 
NYRS - the total number of years of inflow data to 
be input for this particular basin 

REPEAT CARD GROUP 2 FOR (NYRS) YEARS 

2 - NYEAR, CX(J), J=l, 12] 

where NYEAR - calendar year 

X(J) - monthly inflow (m3/s) which is used in 

the monthly water continuity calculations 
for the particular subbasin 
- X(l) corresponds to January data 
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EXAMPLE RUN - 
HISTORICAL REPLICATION 
SOUTH BRANCH MUSKOKA RIVER 



In this appendix the historical replication of the South 
Branch Muskoka River for the period 1975 to 1980 is presented 
as a sample run. 

All input and output files are included. Tables F-1 to F-6 
are the input files while Tables F-7 to F-9 are the output 
files. 

F-1 - Control data file for the historical replication at 
Baysville 

F-2 - Temperature data (celcius) with respect to runoff 
calculations 

F-3 - Precipitation data (M/month) with respect to runoff 
calculations 

F-4 - Precipitation data (M/month) as recorded 

F-5 - Measured basin discharge (m3/s) 

F-6 - Historical month-end reservoir levels (tn) 

F-7 - Output file for historical replication at Baysville. 
Individual monthly summaries are only included for 
January and February 197 5 



I 
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F-8 - Output data file to be used as input to the Serpent 
River Model 

F-9 - Summary of calculated discharges from basin, to be used 
as an input file for the downstream basin 



I 
I 
I 

I 
I 
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TABLE F-1 CONTROL DATA FILE FOR 

HISTORICAL REPLICATION AT BAYSVILLE 



SOUTH BRANCH OF THE MUSKOKA RIVER - DETERMINISTIC RUNOFF EQUATIONS 
VERIFICATION RUNS - ACTUAL LAKE LEVELS USED AT KAWAGAMA AND BAYS 
SIMULATION PERIOD 1975 - 1930 
197S 1 1980 12 
1 1 1 





3 
6 
BURNT 




12 




2 

49.45 3.209 1.523 




9.60 






5 




1 120 243 238 365 




29.37 30.73 30.73 29.26 29.37 


JOE 




2 2 




5 I 1 

53.00 4.293 2.507 
5.80 









9 

1 46 74 95 120 151 253 238 365 




421.55 421.35 421,00 421.00 421,90 421.35 




421.35 421.55 421.55 


KAWA 




3 2 




6 
354.80 33.614 4.276 


31.80 




1 




RAGGED 

4 2 




5 




64.15 6.358 1.257 
6,30 









5 

1 151 227 253 365 




431.31 432.67 432.67 431.45 431.31 












-« 
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TABLE F-i CONTROL DATA FILE FOR 
CONT'D HISTORICAL REPLICATION AT BA'fSyiLLE 

TSC 

5 2 

6 2 2 4 
142.38 10.377 3.113 
12.30 

9 

1 46 74 90 120 151 258 363 
417.53 417.48 417.21 117.21 417.38 417.32 
417.82 417.58 

BAY 

6 2 1 

2 3 3 

S91,3S 37.349 24.494 

68.40 

1 


0.244 0.048 O.GOO 0.025 0.00066 0.000 
0.213 0.042 0.000 0.013 0,04900 0.000 
1.033 -0.039 0.0071 0.000 0.04300 0.000 
-0.926 0.329 0.020 0.000 -0.267 0.262 
■0.384 0.051 0.023 0.000 0.213 
0.071 0.119 0.000 0.000 0.11-1 
0.399 0.064 0.000 0.927 0.000 
0.030 0.031 0.000 0.522 0.000 
0.050 0.170 -0.0086 0.000 -0.022 
0.06-1 0.313 -0.017 0.000 -0.021 
0.060 0.296 -0.012 0.000 0.000 
0.161 0.058 0.000 0.036 0.00130 
1 

6 12 3 4 5 6 

0.000 0.000 0.010 0.020 0.080 0.135 0.14S 0.120 0.075 0.045 0.015 0.000 




-0.108 








-0.074 








-0.066 








-0.033 








0.000 








0.000 








0.000 








0.000 








0.145 


120 





075 



I 
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TABLE F-2 TEMPERATURE DATA (CELCIUS) 

WITH RESPECT TO RUNOFF CALCULATIONS 



1 

6 12 3 4 


5 6 






















12 


























1970 


0.0 


0.0 


1.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


51.0 


94.0 


1971 


14,0 


0.0 


0.0 


0,0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


34.0 


91,0 


1972 


44.0 


0.0 


11,0 


0.0 


0.0 


0,0 


0.0 


0,0 


0.0 


0.0 


13.0 


-15.0 


1971 


35.0 


3,0 


33.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


:55.o 


ai.o 


197-1 


37.0 


0.0 


2.0 


0.0 


0.0 


0.0 


0,0 


0,0 


0.0 


0.0 


31.0 


88.0 


1975 


29.0 


2.0 


2.0 


0.0 


0.0 


0.0 


0.0 


0,0 


0.0 


0,0 


61.0 


153. 


1976 


24.0 


0.0 


22.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


3.0 


47.0 


1977 


3.0 


0,0 


6.0 


0.0 


0.0 


0,0 


0.0 


0.0 


0.0 


0,0 


21.0 


16.0 


1978 


9.0 


i.O 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


19,0 


21.3 


1979 


11.5 


0,0 


9,0 


0.0 


0.0 


0,0 


0.0 


0.0 


0.0 


0,0 


26.0 


107. 


1980 


24.0 


0.0 


15.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0,0 


6.0 


49,0 


1931 


6.0 


4.0 


4.0 


0,0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


10.0 


71,0 



TABLE F-3 PRECIPITATION DATA (M/MONTH) 

UITH RESPECT TO RUNOFF CALCULATIONS 



1 

6 12 3 4 5 6 

12 

1970 0,061 0,051 0,054 0.064 0,087 0.099 0.133 0.]20 0.122 0.038 0.080 0.091 

1971 0.106 0.117 0,099 0,050 0,010 0,060 0,062 0,079 0.064 0.060 0.073 0.100 

1972 0.125 0.112 0.096 0.064 0.074 0,082 0,096 0,108 0.108 0.103 0.033 0.096 

1973 0.124 0.036 0.083 0.063 0.081 0.169 0.137 0,118 0.057 0.061 0.077 0.092 

1974 0.120 0.101 0.086 0,078 0.098 0.073 0,065 0,072 0,099 0,117 0.111 0.074 

1975 0.079 0.103 0,090 0.061 0.0:j7 0.035 0.066 0.073 0.095 0.084 0.093 0.103 

1976 0.099 0.097 0,103 0.069 0.041 0.071 0.061 0.064 0.037 0.091 0.106 0.102 

1977 0.099 0.0B2 0.074 0.053 0.040 0.013 0.0:32 0,083 0,102 0,106 0,126 0,129 

1978 0.084 0.036 0.024 0.043 0.047 0.065 0.072 0.107 0.114 0.106 0.094 0,033 

1979 0.122 0,103 0.070 0,096 0.072 0.084 0.085 0.073 0,090 0.103 0.098 0.072 

1980 0.060 0,057 0,063 0,080 0.074 0.099 0.105 0.127 0.113 0,114 0.078 0.092 
1931 0.103 0.064 0.065 0.073 0.730 0.037 0,069 0,119 O.HO 0,160 0,104 0.068 
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TABLE F-4 PRECIPITATION DATA (M/MONTH) 

(RECORDED MONTHLY PRECIPITATION) 



6 12 3 4 5 6 
12 

1970 0.061 0.042 0.060 0.063 0.106 0.091 0.203 0.065 0.097 0.101 0.071 0.112 

1971 O.iOO 0.1;M 0.064 0.036 0.045 0.075 0.067 0.094 0.031 0.053 0.036 0.114 

1972 0.135 0.089 0.064 0.064 0.083 0.080 0.126 0.119 0.079 0.109 0.064 0.123 

1973 0.119 0.033 0.077 0.039 0.103 0.235 0.072 0.046 0.035 0.0S3 0.093 0.091 

1974 0.149 0.053 0.055 0.100 0.096 0.050 0.049 0.113 0.130 0.104 0.104 0.044 

1975 0.114 0.092 0.066 0.056 0.057 0.053 0.037 0.093 0.105 0.0:55 0.101 0.106 

1976 0.092 0.103 0.114 0.023 0.05S 0.083 0.042 0.067 0.152 0.053 

1977 0.115 0.049 0.059 0.036 0.024 0.062 0.071 0,116 0.113 0.0S<> 

1978 0.057 0.014 0.035 0.052 0.042 0.088 0.0B7 0.145 0.111 0.06;^ 

1979 0.161 0.044 0.091 0.023 0.116 0.053 0.036 0.094 0.090 0.133 

1980 0.065 0.049 0.074 0.087 0.061 0.137 0.117 0,128 0.096 0.120 0.041 0.143 

1981 0.063 0.064 0.068 0.089 0.064 0.090 0.067 0,201 0.182 0.100 0.490 0.043 



0.121 


0.033 


0.147 


0.110 


0.0S2 


0.033 


0.088 


0.055 



TABLE F-5 MEASURED BftSIN DISCHARGE 
CUBIC METERS PER SECOND 



6 12 






1970 


20 


n 


1971 


19 


7 


1972 


19 


1 


1973 


29 


.3 


1974 


25 


9 


1973 


20 


.3 


1976 


25 





1977 


15 


.4 


1973 


29 


7 


19 79 


21 


5 


1980 


30 


2 


1981 


20 


9 



22.2 22.3 28.3 43.6 13.3 49.9 20.4 11.2 30.9 21.3 25,6 25.9 

27.3 34.1 47.5 33.6 11,1 5.13 5.79 13.2 7.55 3.36 10,4 13.7 

24.0 31.6 40.1 64.3 18.5 26.1 31.9 16.0 24.7 31.9 21. S 29.2 
40.9 137.3 53.3 34.7 31,9 17,6 16,5 15,9 16,3 10,7 25,5 29.6 
33.3 40.6 57,1 64.7 19.7 7.64 6.83 12.3 26.1 28.4 22.5 2S.7 

23.1 40,1 35,5 31,5 4,62 3,38 5,06 13.3 20.1 10.4 32.0 20.0 

27.6 50.5 69.1 24.0 8.49 18.6 6.31 6.18 11.8 3.17 12.9 22.4 
13.3 11.3 33,7 13,9 4.04 3.04 3.29 14.4 33.7 24.1 27,5 19.3 
19,5 27,7 25.1 33.4 11.3 3.97 6.32 19.7 25.6 14.7 21.4 21.6 

30.7 50.1 37,1 44.7 15.3 5.09 9.52 22,1 23,3 43.1 40.1 30.6 
25.1 23.0 78.5 20.6 22.5 17.0 22.8 24.3 33.9 30.3 24.3 29.3 
37.1 52.6 31.3 22.3 13,6 5.24 3,96 61.3 37.5 14.6 16.7 26.4 



i 
I 



■ABLE F-6 HISTORICAL .MONTH END 
RESERyaiR LEVELS («> 



F-7 



3 6 








1973 


5.13 






5.05 


4.73 


4 


.44 


1976 


3»08 






4.71 


4.24 


4 


.80 


1977 


4.<S7 






4.10 


3.36 


4 


► 45 


1978 


5.35 






4.73 


4,22 


3 


.72 


1979 


5.09 






5.06 


4.53 


4 


.63 


1980 


4»79 






4.73 


4,27 


4 


,77 


6 3 








1973 


5.13 






5.21 


4.87 


4 


95 


1974 


5.13 






5.17 


4.92 


4 


,25 


1975 


5.06 






5.07 


4.99 


4 


40 


1976 


5.16 






5 . 03 


4.32 


4 


31 


1977 


5.03 






4.92 


4.81 


4 


66 


1973 


5.13 






4.89 


4.79 


4 


40 


1979 


5.18 






5.16 


4.77 


4 


68 


1930 


5.13 






5.05 


4.82 


4 


77 



5,09 5.31 5,37 5,31 3.24 5,08 4.41 4.74 5.08 
5.43 5.59 5.43 5.26 5.08 4,94 4.36 4.23 4.47 
5.38 5,35 5.27 5.21 3.21 4.94 4.65 5.12 5,35 
5.37 5.49 5.61 5.53 5.41 5.06 4.68 4.39 5,09 
5.59 5.61 5.58 5,49 5.44 4,92 4.89 4.92 4.79 
5.49 5.47 5.72 5,72 5.31 5.12 4.36 3.02 5.05 

5.18 5.19 5.23 5.22 5.10 5.20 5.04 5.23 5.18 

5.28 5.13 5.16 5.10 5.19 5.19 5.13 5,22 5.06 
4.88 5.13 5.18 5.22 5.22 5.21 5,16 5.21 5.16 
5.18 5.14 5.24 5.17 5.16 5.21 5.19 5.13 5.03 
5.25 5.23 5.22 5.19 5.25 5.27 5.14 5.23 5.13 
5.21 5.19 5.16 5.18 5.27 5,21 5,16 5.24 5,18 

5.29 5.24 4.88 5.22 5.30 5.14 5.20 5.30 5.18 
5.13 5.19 5.24 5.32 5.25 5,30 5,18 5.1? 5.10 
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TABLE F-7 OUTPUT FILE FOR HISTORIC 
REPLICATION AT BAYSVILLE 
(PAGES F-9 TO f-U) 



SOUTH BRANCH HUSKOKA RIVER - HISTORIC SEQUENTIAL RUNOFF EQUATIONS 
VERIFICATION RUNS - ftCTUrtL LAKE LEVELS USED AT KAaAGAHA AND BAYS 
SIMULATION PERIOD 1975 - 1980 



SIMULATION CONTROL DATA 

START - YEAR = 1975 

MONTH = I 

END - YEAR = 1980 

MONTH 5 12 

IPRTl = I : PRIHT MONTHLY SUMMARY 

IPRT2 = 1 : PRINT BASIN DISCHARGE SUMMARY 

IPRT3 = 1 : PRINT LAKE LEVEL SUMMARY 

lADJl = 3 ; NUMBER OF ITERATIONS FOR MOHTH-BY-MONTH BASIN DISCHARGE ADJUSTMENT 
IADJ2 = : NUMBER OF ITERATIONS FOR LINEAR REGRESSION ADJUSTMENT 

NUMBER OF BASINS = 6 • 



I 
I 
I 
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SOUTH BRANCH MUSKOKA RMR - HISTORIC SEQUEHTIAL RUNOFF EQUATIONS 
VERIFICATION RUNS - ACTUAL LAKE LEVELS USED AT KAUAGAMA AND BATS 
SIMULATION PERIOD 1975 - 1980 



BASIN BURNT - NODE NUMBER 1 

BASIN REGULATION = 2 
EXTERNAL INFLOU FLAG = 
MEASURED BASIN DISCHARGE FLAG - 



DOUHSTREAM BASIN = 2 

NUMBER OF UP3TREAH INFLOWS = 

AREA OF BASIN = 49,4500 KM«I2 
AREA OF LAKES = 3.2090 mttl 
AREA OF SUAHPS = 1.5230 K««2 

AREA OF REGULATED LAKE = 9.4000 K«»I2 
RULE CURVE FOR LAKE OPERATION (M) 

JAN FIB MAR APR HAY JUH JUL AUG SEP OCT NOV DEC 
30.099 30.314 30.551 30.780 30.780 30.780 30.780 30.780 29.830 29.387 29.624 29.870 
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SOUTH BRANCH HUSKOKA RIVER - HISTORIC SEQUENTIAL RUNOFF EQUATIONS 
VERIFICATIOf* RUNS - ACTOAL LAKE LEVELS USED AT KAMAGANA AND BAYS 
SIMULATION PERIOD 1975 - 1980 



BASIN JOE - NODE WHBER 2 



BASIN REGULATION = 2 
EXTERNAL IMFLQU FLAG - 
HEASURED BASIN DISCHAR6E FLAG = 



DOUNSTREAd BASIN = 5 

fflJMBER OF UPSTREAM INFLOyS = 1 

UPSTREAM BASINS = 1 

AREA OF BASIN = 53.0000 KM«2 
AREA OF LAKES = 4.2930 KHt«2 
AREA OF SyAMPS = 2.5070 K«t«2 



AREA OF REGULATED LAKE = 5.8000 KMt«2 
RULE CURVE FOR LAKE OPERATION (H) 



JAN FEB MAR APR MAY JUH JUL AUG SEP OCT NOV DEC 
421.550 421.295 421.000 421.900 421.850 421.850 421.350 421,850 421.700 421.530 421.550 421.550 
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SOUTH BRANCH rtlMOKA RIVER - HISTORIC SEQUENTIAL RUNOFF EQUATIONS 
VERIFICATION RUNS - ACTUAL LAKE LEVaS USED AT mmm AND BAYS 
SIMULATION PERIOD 1975 - 1980 



BASIN KAUA - NODE NUHBER 3 

BASIN REGULATION = 2 
EXTERNAL INFLOU FLAG = 
MEASURED BASIN DISCHARGE FLAG - 



DOKHSTREAM BASIN = 4 
NUMBER OF UPSTREAM INaOUS 



AREA OF BASIN - 354.9000 KNtl2 
AREA OF LAKES '■ 33.4140 KM»2 
AREA OF SUAHPS = 4.27A0 K«»I2 



AREA OF REGULATED LAKE ^ 31.3000 K»1«t2 
MONTH-END LAKE LEVaS (H) 



YEAR 


JAN 


FEB 


HAR 


APR 


MAY 


JUN 


JUL 


AUG 


SEP 


OCT 


NOV 


DEC 


1975 


5.05 


4.73 


4.44 


5.09 


5.31 


5.37 


5.31 


5.24 


5.08 


4.41 


4.74 


5.08 


1976 


4.71 


4.24 


4.80 


5.43 


5.59 


5.43 


5.26 


5.08 


4.94 


4.36 


4.28 


4.47 


1977 


4.10 


3.86 


4.45 


5.38 


5.35 


5.27 


5.21 


5.21 


4.94 


4.65 


5.12 


5.35 


1978 


4.73 


4.22 


3.72 


5.37 


5.49 


3.61 


5.53 


5.41 


5.06 


4.48 


4.89 


5.09 


1979 


5.06 


4.53 


4.63 


5.59 


5.61 


5.58 


5.49 


5.44 


4.92 


4.89 


4.92 


4.79 



1980 4.73 4,27 1.77 5.49 5.47 5,72 5.72 5.31 5.12 4.86 5,02 5,05 
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SOUTH BRANCH MUSKOKA RIVER - HISTORIC SEQUENTIAL RUNOFF EDUATIONS 
VERIFICATION RUNS - ACTUAL LAKE LEVELS USED AT KAUAGAiW AND SAYS 
SIMULATION PERIOD 1975 - 1980 



BASIN RAGGED - NODE NUilBER 4 



BASIN REGULATION = 2 
EXTERNAL INFLOW FLAG - 
MEASURED BASIN DISCHARGE FLAG = 



DO«NSTREA« BASIN = 5 
NUMBER OF UPSTREAM im.O«S 



AREA OF BASIH = 64,1500 KN»I2 
AREA OF LAKES = 6.3580 KN>« 
AREA OF SUAHPS = 1.2570 K««2 



AREA OF REGULATED LAKE = 6.3000 KHIt2 
RULE CURVE FOR LAKE OPERATION («) 

m FEB MAR APR MAY JUN JUL AUG SEP OCT NOV DEC 

431,982 432.143 432.320 432,492 432.670 432.670 432.670 432,040 431.500 431.605 431.706 431.810 
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SOUTH BRANCH MUSKOKA RIVER - HISTORIC SEQUENTIAL RUNOFF EQUATIONS 
VERIFICATION RUMS - ACTUAL LAM LEVELS USED AT KAUAGAWA AND BAYS 
SIHULATION PERIOD 1975 - 1980 



BASIN TSC - NODE MU«BER 5 

BASIN REGULATION = 2 
EXTERNAL INFLOW FLAG = 
MEASURED BASIN DISCHARGE FLAG = 



DOWNSTREAM BASIN = A 

ML'MEER OF UPSTREAM INFLOWS = 2 

UPSTREAH BASINS = 2 4 

AREA OF BASIN = 142.3800 K«lt2 
AREA OF LAKES = 10.3770 K««2 
AREA OF SWAMPS = 3.1130 KN»»2 

AREA OF REGULATED LAKE ^ 12.5000 KM»»2 
RULE CURVE FOR LAKE OPERATION (M) 

M FEB HAR APR MAY JUN JUL AUG SEP OCT NOV DEC 

417.513 417.355 417.210 417.880 417.820 417.820 417.820 417.820 417. 78A 417.717 417. iSO 417.580 



SOUTH BRWO MUSKOKA RMR - HISTORIC SEQUENTIAL RUNOFF EQUATIONS 
MERIFICrtTION RUNS - ACTUAL LAKE LEVELS USED AT KAWAGAtW AND BAYS 
SIMULATION PERIOD 1975 - 1980 



BASIN BAY 



NODE NUHBER 6 



BASIN REGULATION = 2 
EXTERNAL INFLOU FLAG = 
MEASURED BASIN DISCHARGE FLAG = 
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DOUNSTREAM BASIN = 
NWBER OF UPSTREAM INFLOUS = 



UPSTREAM BASINS 

AREA OF BASIN = 
AREA OF LAKES = 
AREA OF SyAMPS - 



5 3 

591.3500 K««2 
37.8490 KMtt2 
24.4940 K««2 



AREA OF REGULATED LAKE = i3.400O KM»2 
MONTH-END LAKE LEVELS (M) 



YEAR 


JAN 


FEB 


MAR 


APR 


MAY 


JUN 


JUL 


AUG 


SEP 


OCT 


nm 


DEC 


1973 


5.21 


4,87 


4.95 


5.18 


5.19 


5.28 


5.22 


5,10 


5.20 


5.04 


5.23 


5.18 


1974 


5.17 


4.92 


4.25 


5.28 


5,13 


5.16 


5.10 


5.19 


5.19 


5.13 


5,22 


5.06 


1975 


5.07 


4,99 


4.40 


4.88 


5.13 


5.18 


5.22 


5,22 


5.21 


5.16 


5.21 


5.16 


1976 


5.03 


4,82 


4.81 


5.18 


5.14 


5.24 


5.17 


5.16 


5.21 


5.19 


5,18 


5.08 


1977 


4.92 


4.81 


4.^ 


5.25 


5.23 


5.22 


5.19 


5,25 


5.27 


5.14 


5.23 


5.13 


1978 


4.89 


4.79 


4.40 


5.21 


5.19 


5.16 


5.18 


5.27 


5.21 


5.16 


5.24 


5.18 


1979 


5.U 


4.77 


4.68 


5.29 


5.24 


4.88 


5.22 


5.30 


5.14 


5.20 


5.30 


5.18 


1980 


5.05 


4.82 


4.77 


5.13 


5.19 


5.24 


3.32 


3.25 


5.30 


5.18 


5.19 


5.10 
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SOUTH BRANCH MUSKOKA RIVER - HISTORIC SEQUENTIAL RUNOFF EQUATIONS 
VERIFICATION RUNS - ACTUAL LAKE LEVaS USED AT KAWAGA.1A AMD BAYS 
SIHULATION PERIOD 1975 - 1980 



COEFFICIENTS FOR RUNOFF PREDICTION EQUATIONS - HISTORIC SEQUENTIAL EQUATIONS 



MONTH 


Al 


A2 


A3 


A4 


A5 


M 


JAN 


0.244000 


0.048000 


o.oooooo 


0.025000 


0,000660 


0.000000 


FEB 


0.213000 


0.042000 


0. 000000 


0.013000 


0.049000 


0,000000 


NAR 


1.033000 


-0.039000 


0. 007100 


0.000000 


0,043000 


0.000000 


APR 


-0,926000 


0.529000 


0.020000 


0.000000 


-0,267000 


0.262000 


MAY 


-0.384000 


0.051000 


0.023000 


o.oooooo 


0.213000 


-0,108000 


JUN 


0,071000 


0.119000 


0.000000 


0,000000 


0,114000 


-0.074000 


JUL 


0.399000 


0. 064000 


o.oooooo 


0.927000 


O.OOOOOO 


-0.066000 


AUG 


0.080000 


0.051000 


0.000000 


0.522000 


0.000000 


-0.038000 


SEP 


0.050000 


0,170000 


-O.OOSiOO 


0.000000 


-0,022000 


O.OOOOOO 


OCT 


0.064000 


0.313000 


-0.017000 


0.000000 


-0.024000 


O.OOOOOO 


NOV 


0.060000 


0.296000 


-0.012000 


0.000000 


0.000000 


0,000000 


DEC 


0.161000 


0.058000 


0,000000 


0.036000 


0.001300 


0.000000 
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SOUTH URANCH HUSKOKA RIVER - HISTORIC SEQUENTIAL RUNOFF EQUATIONS 
VERIFICATION RUNS - ACTUAL LAKE LEVELS USED AT KAUA6AM AND BAYS 

3IHULATI0N PERIOD 1975 - 1990 



EVAPORATION COEFFICIENTS (fl/MONTH) 

OATA I APPLICABLE FOR BASINS I JAM FEB MAR APR HAT M JUL AUG SEP OCT MflV DEC 
SET I I 



II 12 3 4 5 6 I 0.0000 0.0000 O.OIOO 0.0200 0.0800 0.1350 0.145O 0.1200 0.0750 0.0450 0.0150 0.00< 



SOUTH BRANCH «USKOKA RIVER - HISTORIC SEQUENTIAL RUNOFF EQUATIONS 
VERIFICATiaN RUNS - ACTUAL LAKE LEVELS UbT;D AT KAWAGArtA AND BAYS 
SIMULATION PERIOD 1975 - 1930 
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SET 1 : TEMPERATURE DATA HRT RUNOFF CALCllATIOfI 



APPLICABLE FOR 6 BASINS 
BASI«S = 12 3 4 5 4 



YEAR 



JAN 



FEB 



MAR 



APR 



MAY 



JUN 



JUL 



AUG 



SEP 



OCT 



NOV 



DEC 



1970 


0.00 


0.00 


1.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


51.00 


94,00 


1971 


14.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


o.oo 


0.00 


0.00 


34.00 


91.00 


1972 


44.00 


0.00 


11.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


13.00 


45.00 


1973 


35.00 


3.00 


33,00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


25.00 


81.00 


1974 


37.00 


0.00 


2.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


31.00 


88.00 


1975 


29.00 


2. GO 


2.00 


0,00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


61.00 


153.00 


1976 


24.00 


0.00 


22.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


3.00 


47.00 


1977 


3.00 


0.00 


6.00 


0.00 


0.00 


0.00 


0.00 


o.oo 


0.00 


0.00 


21.00 


16.00 


1978 


9.00 


1.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0,00 


0.00 


19.00 


21.80 


1979 


11.50 


0.00 


9.00 


0.00 


0.00 


0.00 


0.00 


o.oo 


0.00 


0.00 


26.00 


107.00 


1980 


24.00 


0.00 


15.00 


0.00 


0.00 


0.00 


0.00 


0,00 


0.00 


0.00 


6.00 


49.00 


1981 


6,00 


4,00 


4,00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


10.00 


71.00 
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SOUTH BRANCH NUSKOKA RIVER - HISTORIC SEQUENTIAL RUNOFF EQUATIONS 
VERIFICATION RUNS - ACTUAL LAKE LEVELS USED AT KAUAGAMA AND BAYS 
SIMULATION PERIOD 1?75 - 1980 



SET 1 : PRECIPITATION DATA (M/fWNTH) URT RUNOFF CALCULATION 



YEAR 



APPLICABLE FOR i BASINS 
BASINS = 12 3 4 5 6 



JAN 



FEB 



NAR 



APR 



MAY 



JUN 



JUL 



AUG 



SEP 



OCT 



NOV 



DEC 



1970 
1971 
1972 
1973 
1974 
1975 
1974 
1977 
1978 
1979 
1980 
1981 



0,0610 
0.1040 
0.1250 
0.1240 
0.1200 
0.0790 
0.0990 
0,0990 
0.0840 
0.1220 
0.040O 
0.1030 



0.0510 
0.1170 
0.1120 
0,0840 
0.1010 
0.1030 
0.0970 
0.0820 
0.0340 
O.103O 
0.0570 
0.0440 



0.0540 
0.0990 
0.0940 
0.0830 
0.0840 
0.0900 
0,1030 
0.0740 
0.0240 
0.0700 
0.0430 
0.0450 



0.0440 
0,0500 
0,0440 
0.0480 
0.0780 
0,0410 
0.0490 
0.0580 
0.0430 
0.0940 
0.0800 
0.0730 



0.0870 
0.0400 
0.0740 
0.0810 
0.0980 
0.0570 
0.0410 
0.0400 
0.0470 
0.O72O 
0.0740 
0.730O 



0.0990 
0.0400 
0.0820 
0.1490 
0.0730 
0.0550 
0.0710 
0,0«0 
0.0450 
0,0840 
0.0990 
0.0870 



0.1330 
0.0420 
0.0940 
0.1370 
0.0450 
0.0440 
0.0410 
0.0520 
0.0720 
0.0850 
0.1050 
0.0490 



0.120O 
0.0790 
0.1080 
0.1180 
O.072O 
0.0780 
0.0440 
0.0830 
0.1070 
0.0780 
0.1270 
0.1190 



0.0440 
0.1080 
0.0570 
0.0990 
0.0950 
0.0870 
0.1020 
0.1140 
0.09O0 
0.1130 
0.1400 



0.0880 
0.0400 
0.1030 
0.0410 
0.1170 
0.0840 
0,0910 
0.1040 
0.1040 
0.1080 
0.1140 
0.1400 



0.0800 
0.0730 
0.0830 
0.0770 
0.1110 
0.0930 
0.1040 
0.1240 
0.0940 
0.0980 
0.0780 
0.1040 



0.0910 
0.1000 
0.0960 
0.0920 
0.0740 
0.1030 
0.1020 
0.1290 
0.0830 
0.0720 
0.0920 
0.0680 



SOUTH BRANCH MUSKOKA RIVER - HISTORIC SEWENTIAL RUNOFF EQUATIONS 
VERIFICATION RUNS - ACTUAL LAKE LEVELS USED AT KAUAGA«A AND BAYS 
SIHULATIQN PERIOD 1975 - 1980 



SET 1 ! MONTHLY PRECIPITATION DATA (h/«ONTH) 



APPLICABLE FOR 6 BASIMS 
BASIHS = 12 3 4 5 
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YEAR 



JAN 



FEB 



MAR 



APR 



HAY 



JUN 



JUL 



AUG 



SEP 



OCT 



MOV 



DEC 



1970 


0.0610 


0.0420 


0.0600 


0.0680 


0.1060 


0.0910 


0.2030 


0.0650 


0.0970 


0.1010 


0.0710 


0.1120 


1971 


0.1000 


0.1340 


0.0640 


0.0360 


0,0450 


0.0750 


0.0670 


0.0940 


0.0310 


0.0550 


0.0860 


0.1140 


1972 


0.1350 


0.0890 


0.0640 


0.0640 


0.0830 


O.OSOO 


0.1260 


0.1190 


0.0790 


0.1090 


0.0640 


0.1280 


1973 


0.1190 


0.0530 


0.0770 


0,0590 


0.1030 


0.2350 


0.0720 


0,0460 


0.0550 


0.0830 


0.0930 


0.0910 


1974 


0.1490 


0.0530 


0.05S0 


0.1000 


0,0960 


0.0500 


0.0490 


0.1180 


0.1300 


0.1040 


0.1040 


0.044O 


1975 


0.1140 


0.0920 


0.0660 


0.0560 


0.0570 


0.0530 


0.0870 


0.0930 


0.1050 


0.0550 


0.1010 


0.1060 


1976 


0.0920 


0.1030 


0.1140 


0.0230 


0.0580 


0.0830 


0.0420 


0.0670 


0.1520 


0.0530 


0.1210 


0.0830 


1977 


0.1150 


0.0490 


0.0590 


0.0560 


0.0240 


0.0620 


0.0710 


0.1160 


O.1180 


0.0840 


0.1470 


O.UOO 


1978 


0,0570 


0.0140 


0.0350 


0.0520 


0.0420 


0.0880 


0.0870 


0.1450 


0.1110 


0.0620 


0.0820 


0.0830 


1979 


0.1410 


0.0440 


0.0910 


0.02S0 


0.1160 


0.0530 


0.0860 


0.0940 


0.0900 


0.1380 


0.0880 


0,0550 


1980 


0.0650 


0.0490 


0.0740 


0,0870 


0.0610 


0.1370 


0,1170 


0.1280 


0.0960 


0.1200 


0.0410 


0.1430 


1981 


0.0630 


0.0640 


0.0680 


0.0890 


0.0640 


0.0900 


0.0670 


0.2010 


0.1820 


0.1000 


0.4900 


0.0430 
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SOUTH BRWO MUSKOKft RIVER - HISTORIC SEQUENTIAL RUKOFF EQUATIONS 
VERIFICATION RUNS - ACTUAL LAKE LEVELS USED AT KAiAGAMA ANi BAYS 
SIMULATION PERIOD 1975 - 19B0 





MEASUKIU BASIN DISCHARGE («tl3/S) FOR BASIN BAY 


- NODE NUnSER 4 


AUG 


SEP 


OCT 


NOV 




YEAR 


JAN 


FEB 


MAR 


APR 


MAY 


M 


JUL 


DEC 


1970 


20.2000 


22.2000 


22.3000 


28.3000 


43,4000 


13,8000 


49.9000 


20.4000 


11.2000 


30.9000 


21,3000 


25,4000 


1971 


19,7000 


27.3000 


34.1000 


47,5000 


33,4000 


11.1000 


5,1800 


5.7900 


13.2000 


7.5500 


8,3400 


10.4000 


1972 


19,2000 


24.0000 


31.4000 


40.1000 


44.3000 


18.5000 


24.1000 


31.9000 


14.0000 


24.7000 


31.9000 


21.8000 


1973 


29,5000 


40.9000 


57.3000 


58.8000 


34.7000 


31.9000 


17,4000 


14.5000 


15.9000 


16.8000 


10,7000 


25,5000 


1974 


25.9000 


33.3000 


40.4000 


57.1000 


44,7000 


19.7000 


7.4400 


4.8300 


12.3000 


24.1000 


28,4000 


22.5000 


1975 


20.3000 


23.1000 


40.1000 


35.5000 


31.5000 


4.4200 


3.3800 


5.0400 


13.8000 


20.1000 


10.4000 


32,0000 


1976 


25.0000 


27.4000 


50.5000 


49,1000 


21.0000 


8.1900 


13.4000 


4.3100 


4.1800 


11.8000 


8.1700 


12,9000 


1977 


15.4000 


13,8000 


41.5000 


38.7000 


13.9000 


4.0400 


3.0400 


3.2900 


14.4000 


33.7000 


24.1000 


27.5000 


1978 


29.7000 


19.5000 


27,7000 


25.1000 


53.4000 


11.3000 


3.9700 


4.3200 


19.7000 


25.4000 


14,7000 


21,4000 


1979 


21,5000 


30.7000 


50.1000 


57.4000 


44,7000 


15.8000 


5.0900 


9.5200 


22.1000 


28.3000 


13.1000 


40.1000 


1980 


30.2000 


25.1000 


23.0000 


73,5000 


20.4000 


22.5000 


17,0000 


22.8000 


2^.3000 


33.9000 


30.3000 


24,3000 


1981 


20.9000 


37.1000 


52.4000 


31 .3000 


22.3000 


13.4000 


5.2400 


3.9400 


61.3000 


37.5000 


H.4000 


16.7000 
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SOUTH BRANCH MUSKOKA RMR - HISTORIC SEQUENTIAL RUNOFF EQUATIONS 
VERIFICATION RUNS - ACTUAL LAKE LEVELS USED AT KAUA6AMA AND BAYS 
SIMULATION PERIOD 1975 - 1980 



ITERATION ! SUNHARY OF .iONTHLY SlrtULATIGN - YEAR 1975 , JAN 



BASIN 


NET BASIN 


RUNOFF 


RUNOFF 


BASIN 


UNIT 


■JPSTREAW 


BASIN 


LAKE ELEVATION 


CHANGE n 


mm NODE 


EVAP LOSS 


TEHP 


KKtCIP 


RUNOFF 


RUNOFF 


INFLOW 


DISCHARGE 


INITIAL 


FINAL 


LAKE ELR 






(M) 


DEGREE C 


(M) 


(Mtt3/S) 


(MM) 


(«tl3/S) 


(M»3/S) 


(N) 


(M) 


(H) 


BURNT 




-0.1140 


29.0000 


0.0790 


0.5705 


30.8987 


0.0000 


0.1568 


29.8700 


30.0994 


0.229* 


JOE 




-0.1140 


29.0000 


0.0790 


0.6349 


32.0841 


0.1568 


1.0385 


421.5500 


421.5500 


0.000( 


KAUA 




-0.1140 


29.0000 


0.0790 


3.0947 


23.3619 


0.0000 


5.9916 


5.1800 


5.0500 


-0.130( 


RAGGED 




-0.1140 


29,0000 


0.0790 


0.7969 


33.2714 


O.OOOO 


0.660b 


431.8100 


431.9820 


0.172( 


TSC 




-0.1140 


29,0000 


0.0790 


1.2386 


23.2998 


1.6990 


3.7807 


417.5800 


117.5133 


-0.066; 


BAY 




-0.1140 


29.0000 


O.OTfO 


4.1228 


18.6735 


9.7723 


16.5510 


5.0600 


5.0700 


0.010( 
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SOUTH BRANCH HUSKOKA RIVER - HISTORIC SEQUENTIAL RUHOFF EQUATIONS 
VERIFICATION RUNS - ACTUAL LAKE LEVELS USED AT KAWGAi'lA AND BAYS 
SIMULATION PERIOD 1975 - 1980 



ITERATION 3 ! SUMMARY OF MONTHLY SIMULATION - YEAR 1975 t JAN 



■ BASIN 


NET BASIN 


RUNOFF 


RUNOFF 


BASIN 


UNIT 


UPSTREAM 


BASIN 


LAKE ELEVATION 


CHANGE Ii 


1 NAME NODE 


EVAP LOSS 


TEMP 


PRECIP 


RUNOFF 


RUNOFF 


INFLDU 


DISCHARGE 


INITIAL 


FINAL 


LAKE aE 






(M) 


DEGREE C 


(M) 


(M«3/S) 


(HN) 


(N«3/S) 


(«tS3/S) 


{«) 


(M) 


(H) 


■ BURNT 


I 


-0.1140 


29.0000 


0.0790 


0.7750 


11.9748 


0.0000 


0.3613 


39.8700 


30.0994 


0.229 


JOE 




tL 


-0.1140 


29.0000 


0.0790 


0.8625 


43.5852 


0.3613 


1.4706 


421.5500 


421.5500 


0.000 


■ KAUA 


3 


-0.1140 


29.0000 


0.0790 


4.2040 


31.73i3 


0.0000 


7.1010 


5.1300 


3.0500 


-0.130' 


■ RAGGED 


4 


-0.1140 


29.0000 


0.0790 


1.0825 


45.1981 


0.0000 


0.9461 


431.3100 


431.9820 


0.172' 


TSC 


5 


-0.1140 


29.0000 


0,0790 


1.682& 


31.6520 


2.4167 


4.9424 


417.5800 


417.5133 


-0,066 


« BAY 


6 


-0.1140 


29.0000 


0.0790 


5.4007 


25.3673 


12.0434 


20.3000 


5.0600 


5.0700 


0.010 
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SOUTH BRANCH HUSKOKA RIVER - HISTORIC SEQUENTIAL RUNOFF EQUATIONS 
VERIFICATION RUNS - ACTUAL LAKE LEVaS USEP AT KAUAGA«A AND BAYS 
SIMULATION PERIOD 1975 - 1980 



ITERATION ! SUHMARY OF MONTHLY SIHULATIDN - YEAR 1975 t FEB 



BASIN 


NET BASIN 


RUNOFF 


RUNOFF 


BASIN 


UNIT 


UPSTREAH 


BASIN 


LAKE ELEVATION 


CHANGE Ih 


NAME NODE 


EVAP LOSS 


TEMP 


PRECIP 


RUNOFF 


RUNOFF 


iNaou 


DISCHARGE 


INITIAL 


FINAL 


LAKE ELE\ 






(M) 


DEGREE C 


(«) 


(Mt«3/S) 


(MH) 


(««3/S) 


(M«37S) 


(N) 


(«) 


(H) 


BURNT 


1 


-0.0920 


2.0000 


0.1030 


0.93B6 


45.9m 


0.0000 


0.4540 


30.0994 


30.3135 


0.2141 


JOE 


2 


-0.0720 


2,0000 


0.1030 


1.0266 


44.8572 


0,4540 


2.3133 


421.5500 


421,2946 


-0.255- 


KAUA 


3 


-0.0920 


2.0000 


0.1030 


5.8i39 


39.9832 


0.0000 


11.2796 


5.0500 


4.7300 


-0.320( 


RAGGED 


4 


-0.0920 


2.0000 


0.103O 


1.2475 


17.7995 


0.0000 


1.0890 


431.9820 


432.1425 


0.1404 


TSC 


5 


-0.0920 


2.0000 


0.1030 


2.3500 


39.9290 


3.4023 


7.0474 


417.5133 


417.3544 


-0.158: 


BAY 


4 


-0.0920 


2.0000 


0.1030 


8.34S1 


34.2792 


18.3272 


32.0584 


5.0700 


4.9900 


-0.080( 
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SOUTH BRANCH WSKOKA RIVER - HISTORIC SEQUENTIAL RUNOFF EOUATIONS 
VERIFICATION SUNS - ACTUAL LAKE LEVELS USED AT KAWAGAHA AND BAYS 
SIHULATION PERIOD 1975 - 1980 



ITERATION 3 ! SU«HART OF MONTHLY SIMULATION - YEAR 1975 , FEE 



BASIN 


NET BASIN 


RUNOFF 


RUNOFF 


BASIN 


UNIT 


UPSTREAM 


BASIN 


LAKE aEVATION 


CHANGE li 


NAHE NODE 


EVAP LOSS 


TEMP 


PfJECIP 


RUNOFF 


RUNOFF 


iNao« 


DISCHARGE 


INITIAL 


FINAL 


LAKE ELE 






(«) 


DEGREE C 


(«) 


CMIH/S) 


(MN) 


(MW3/S) 


(MI13/S) 


(N) 


(N) 


(H) 


BURNT 


1 


-0.0920 


2.0000 


0.1030 


0.5247 


25.W93 


0.0000 


0.0401 


30.0994 


30.3135 


0.214 


JOE 


It 


-0.0920 


2.0000 


0.1030 


0.5739 


26.1940 


0.0401 


l.'!467 


421.5500 


421.2946 


-0.255: 


KAUA 


3 


-0.0920 


2.0000 


0.1030 


3.2780 


22.3513 


0.0000 


8.6937 


5.0500 


4.7300 


-0.320' 


RAGGED 


4 


-0.0920 


2.0000 


0.1030 


0.7086 


26.7207 


0.0000 


0.5300 


431.9820 


432.1425 


0,160. 


TSC 


5 


-0.0920 


2.0000 


0.1030 


1.3137 


22.3210 


1.9767 


1.5857 


417.5133 


417.3546 


-0.158: 


BAY 


& 


-0.0920 


2.0000 


0.1030 


4.9574 


20.2807 


13.2795 


23.1000 


5.0700 


4.9900 


-o.oeo 
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SOUTH BRANCH MUSKOKA RIVER - HISTORIC SEQUENTIAL RUNOFF EQUATIONS 
VERIFICATION RUMS - ACTUAL LAKE LEVELS USED AT KAUAGAMA AND BATS 
SIMULATION PERIOD 1975 - 1980 



ITERATION 3 ! DISCHARGE (««3/S) SUMMARY FOR BASIN - BURNT f NODE NUMBER 1 



YEAR TYPE M FEB MAR APR MAY JUN JUL AU6 SEP OCT NOV DEC MEAN 



1975 CALC 


0.361 


0.040 


0.811 


0.936 


1.182 


0.116 


0.209 


0.083 


3.846 


2.008 


0.000 


0.712 


0.859 


1976 CALC 


0.236 


0.113 


1.871 


1.814 


0.705 


0.265 


0.938 


0.042 


3.634 


1.713 


0.000 


0.000 


0.944 


1977 CALC 


0.000 


0.000 


0.905 


1.104 


0.329 


O.OOO 


0.090 


0.139 


3.860 


2.713 


0.678 


0.549 


0.864 


1978 CALC 


0.181 


0.000 


0.157 


0.539 


1.956 


0.381 


0.185 


0.276 


3.942 


2.376 


0.064 


0.288 


0.862 


1979 CALC 


0.478 


0.000 


1.935 


2.052 


1.407 


0.090 


0.684 


0.380 


3.828 


2.896 


1.180 


0.772 


1.308 


19S0 CALC 


0.757 


0.000 


0.575 


2.388 


0.631 


1.044 


0.987 


0.612 


4,336 


2.788 


0.571 


0.341 


1.253 



MEAN CALC 0.336 O.026 1.042 1.472 1.035 0.316 0.516 0.255 3.908 2.416 0.416 0.444 
SDEV CALC 0.263 0.046 0.715 0.719 0.597 0.382 0.403 0.215 0.233 0.473 0.479 0.291 
SKiy CALC 0.586 1.901 0.336 -0.011 0.511 1.307 0.197 0.931 1.354 -0.642 0.770 -0.444 
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SOUTH BRANCH MUSKOKA RIVER - HISTORIC SEQUENTIAL RUNOFF EQUATIONS 










VERIFICATION RUNS - ACTUAL LAKE LEVaS USEB AT KAlWGAm AND BAYS 












SlrtHAlIOM PERIOB 1975 - 1980 






























YEAR TYPE 


ITERATION 3 : DISCHAfil^ <Nlt3/S) SUifilARY Fi]R BASIN - 


JOE i MDE NUMBER 2 


NOV 


BEG 




JAN FEB MAR APR MAY JUN 


JUL AUG 


SEP OCT 


MEAN 




1975 CALC 


1.471 1.447 3.039 0.368 2.698 0.401 


0.511 0.215 


4.466 2.715 


0.731 


2.330 


1.743 


1976 CALC 


1.244 1.593 5.209 2.777 1.671 0.655 


2.041 0.168 


3.973 2.152 


0.361 


0.667 


1.876 


1977 CALC 


0.410 1.030 4.411 1.232 0.928 0>119 


0.303 0.294 


4.478 4.146 


2.116 


1.951 


1,785 




1978 CALC 


1.196 1.189 2.029 0.076 4.398 0,393 


0.465 0.533 


4.656 3.493 


0.957 


1.455 


1.778 




1979 CALC 


1.631 1.344 5.426 3.226 3.064 0.341 


1.495 0.821 


4.451 4.436 


3.253 


2.507 


2.667 


1980 CALC 


2.402 1.356 2.686 3.848 l.SOl 2.191 


2,088 1.273 


5.478 4,247 


2.072 


1.144 


2.549 



I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 



MEAN CALC 1.392 1.327 3.800 2.004 2.377 0.767 1.150 0.551 4.584 3.537 1.582 1.726 

SDEV CALC 0.649 0.197 1.412 1.490 1.267 0.717 0.823 0.429 0.494 0.920 1.088 0.678 

SKEU CALC 0.094 -0.314 -0.006 -0.050 0.672 1.792 0.203 1.081 1.202 -0.708 0.526 -0.500 
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SOUTH BRANCH MUSKOKA RIVER - HISTORIC SEQl^NTIAL RUNOFF EQUATIONS 
yERIFICATION RUMS - ACTUAL LAKE LEVELS USEP AT KAWAGAftA AND BAYS 
SIHULATION PERIOD 1973 - 1980 



IT 


ERATION 3 


: DISCHAR6E (MW3/S) SIMARY FOR BASIN - 


KAUA ) 


MODE NUMBER 3 










TEAR TYPE 


JAN 

7.101 


FEB 
8.694 


«AR 
8.542 


APR 
5.905 


HAY "^ 


JUN 


JUL 


AUG 


SEP 


OCT 


NOV 
0.404 


DEC 
6.000 


HE AN 


1975 CALC 


10.957 


2.088 


2.201 


1.898 


3.549 


10.270 


5.634 


1974 CALC 


9.494 


10.398 


3.304 


13.500 


4.507 


5.383 


7.403 


3.244 


1.746 


7.552 


3.074 


1.661 


6.567 


1977 CALC 


6.521 


5.167 


4.426 


5.723 


3.842 


2.445 


1.598 


1.144 


4.909 


10.156 


3.225 


5.445 


4.567 


1978 CALC 


12.254 


9.474 


8.361 


4.723 


8.150 


2.925 


2.380 


3.195 


6.509 


9.324 


2.936 


4.318 


6.213 


1979 CALC 


4.244 


10.454 


11.427 


10.483 


14.058 


2.790 


5.494 


3.935 


7.944 


7.474 


12.635 


12.066 


8.7S4 


1980 CALC 


8.557 


9,475 


1,364 


15.462 


7.581 


5.336 


6.3^ 


9,632 


7.099 


9,991 


7,253 


6.090 


7.889 



MEAN CALC 8.362 8.943 7.221 9.299 8.516 3.578 4.233 3.845 5.294 9.161 4.921 5.930 
SDEV CALC 2.278 1.964 3.633 4.523 3.544 1.407 2,470 3.011 2.341 1.244 4.372 3.427 
SKEU CALC 1.084 -1.864 -0.782 0.419 0.487 0,858 0.197 1.827 -0.542 -0.708 1.294 1.110 



I 
I 
I 
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SOUTH BRANCH MUSKOKA RIVER - HISTORIC SEQUENTIAL RUHOFF EQUATIONS 














VERIFICATION RUNS - ACTUAL LAKE LEVaS USED AT KAUAGAHA AND BAYS 
















SIHULATION PERIOD 1975 - 1980 














































YEAR TYPE 


iimmm z 


I DISCHARGE (!1tt3/S) SUMHARY FOR BASIN - 1 


RAGGED ; 


NODE NUHBER 4 


OCT 


NOV 


DEC 




JAN 


FEB 


HAR 


APR 


MAY M 


JUL 


AUG 


SEP 


HEAft 


1975 CALC 


0.946 


0.530 


1.241 


0.943 


1.852 0.565 


0.461 


1.709 


1,620 


0.204 


0.603 


1.736 


1.033 


1976 CALC 


0.826 


0.633 


2.713 


1.840 


0.983 0,676 


1.611 


1.733 


1.323 


0.000 


0.040 


0.572 


1.079 


1977 CALC 


0,063 


0.089 


2.695 


1,005 


0,432 0.298 


0.34^ 


1.721 


1.620 


1.01? 


1.448 


1.514 


1.021 




1978 CALC 


0.852 


0.302 


0.800 


0.207 


3.359 0,873 


0.431 


1.842 


1.734 


0.660 


0.806 


1.234 


1,092 




1979 CALC 


0.982 


0.48^ 


2.731 


2,4^9 


1.998 0.474 


1.177 


2.187 


1.613 


1.137 


2.219 


1.934 


1.618 




1980 CALC 


1.643 


0.503 


1.236 


2.510 


0,810 1.737 


1.565 


2.448 


2.214 


1.053 


1.536 


1.106 


1.530 



I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 



ICAN CALC 0.386 0.424 1,903 1.496 1.572 0.771 0.932 1.940 1.488 0,679 1.110 1.350 
SDEV CALC 0.502 0.196 0.902 0.927 1.066 0.511 0.590 0.308 0.292 0.480 0.777 0.489 
SKEW CALC -0.268 -1.127 -0.134 -0.159 0.911 1.713 0,220 1.171 1.156 -0,629 0.066 -0,574 
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SOUTH BRANCH HIiSKOKA RIVER - HISTORIC SEaUEMTIAL RUWF EQUATIONS 
VERIFICATION RUNS - ACTUAL LAKE LEVELS USEB AT KAUAGAHA AND BAYS 
SIMULATION PERIOD 1975 - 1930 



ITERATION 3 '. DISCHARGE (««3/S) SUHHARY FOR MSIN - TSC i NOBE NUKBER 5 



YEAR TYPE JAN FEB MAR APR MAY JUH JUL AUG SEP OCT MOV DEC MEAN 



1975 CALC 


4.942 


4.586 


7.653 


5.030 


9.345 


1.338 


1.549 


2.287 


6.942 


4.299 


3.464 


8.340 


5.023 


1974 CALC 


4,415 


4.724 


15.074 


11.114 


5.715 


2.453 


5.766 


2.271 


5.412 


2.763 


1.559 


3.084 


5.363 


1977 CALC 


l.d31 


2.739 


13.351 


6.211 


3.328 


O.B30 


1.009 


2.416 


6.940 


8.320 


7.513 


6.937 


5.102 


1978 CALC 


4.312 


3.413 


5.176 


2.781 


15.543 


3.203 


1.413 


3.012 


7.493 


6.512 


4.323 


5.572 


5.233 


1979 CALC 


5.267 


4.041 


14.707 


12.267 


10.1S8 


1.596 


4.347 


4.176 


6.938 


8.959 


11.079 


8.Ki7 


7,699 


1980 CALC 


7.489 


4.043 


7.811 


14.001 


5.042 


7.025 


6.012 


5.410 


9.910 


8.503 


7.709 


5.373 


7.3A1 



MEAN CALC 4.676 3.925 10.629 8.5A8 8,188 2.825 3.334 3.262 7.273 6.559 5.941 6.360 

SDEV CALC 1.887 0.744 4.250 4.500 4.448 2,209 2,292 1.277 1,470 2.536 3,463 2.136 

SKEU CALC -0.259 -0.729 -0.137 -0.077 0,838 1,764 0.200 1.173 1.116 -0.729 0.299 -0.404 



I 
I 
I 
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SOUTH BRANCH HUSKOKA RIVER - HISTORIC SEQUENTIAL RUNOFF EQUATIONS 
VERIFICATION RUNS - ACTUAL LAKE LEVELS USED AT KAWAGAHA AND BAYS 
SIMULATION PERIOD 1975 - 1980 





ITERATION ] 


; DISCHARGE (tt«3/S) SUMHARY FOR BASIN - 


BAY 1 


NODE NUMBER 6 










YEAR TYPE 


JAN 


FEB 


HAR 


APR 


NAY 


M 


JUL 


AUG 


SEP 


OCT 


NOV 


DEC 


NEAN 


1975 CALC 


20.300 


23.100 


40.100 


35.404 


31.500 


4.420 


3.380 


5,040 


13.800 


20.100 


10.457 


32.000 


19.985 


MEAS 


20.300 


23.100 


40.100 


35.500 


31.500 


^420 


3.380 


5,040 


13.800 


20.100 


10.400 


32.000 


19.988 


1974 CALC 


25.000 


27.593 


50.500 


69.100 


24.000 


8.490 


18.400 


4.310 


4.130 


12.050 


9.010 


13.700 


22.544 


«EAS 


25.000 


27.400 


50.500 


49.100 


24.000 


8.490 


18.400 


6,310 


4,180 


11.800 


8.170 


12.900 


22.387 


1977 CALC 


16.145 


14.105 


41.500 


38.444 


13.700 


4.058 


3.040 


3.290 


14.400 


33.700 


24.100 


27.500 


19.517 


HEAS 


15.400 


13.800 


41.500 


38.700 


13.900 


4.040 


3.040 


3.290 


14.400 


33.700 


2M0O 


27.500 


19.448 


1978 CALC 


29.409 


19.480 


27.700 


22.920 


53.400 


11.300 


3.970 


4.320 


19.700 


25.400 


14.700 


21,400 


21.342 


NEAS 


29.700 


19.500 


27.700 


25.100 


53.400 


11.300 


3.970 


6.320 


19.700 


25.600 


14.700 


21.400 


21.533 


1979 CALC 


21.500 


30.459 


50.474 


57.400 


44,700 


15,800 


5.090 


9.520 


22.100 


28,300 


43.100 


40.100 


30.712 


HEAS 


21.500 


30.700 


50.100 


57.400 


44.700 


15.800 


5.090 


9.520 


22.100 


28,300 


43.100 


40.100 


30.701 


1980 CALC 


30.130 


24.782 


23.000 


78.500 


20.400 


22,500 


17.000 


22.800 


24.300 


33.900 


30.124 


24.300 


29.328 


HEAS 


30.200 


25.100 


23.000 


78.500 


20.400 


22.500 


17.000 


22.800 


24.300 


33.900 


30.300 


24.300 


29.375 


MEAN CALC 


23.784 


23.253 


38.879 


50.295 


31.350 


11.128 


8.513 


8.883 


16.747 


25.608 


21.915 


26.500 




NEAN HEAS 


23.483 


23.300 


38.817 


30.717 


31.350 


11.125 


8.513 


8.883 


14.747 


25.567 


21. m 


24.367 




SDEV CALC 


5.497 


5.849 


11.444 


21.500 


15.090 


7,078 


7.245 


7.114 


4.437 


3.434 


13.197 


9.065 




SDEV HEAS 


5.749 


6.020 


11.370 


20.913 


15.090 


7,081 


7.215 


7.114 


4.437 


3.514 


13.397 


9.294 




SKEW CALC 


-0.075 


-0.540 


-0.451 


0.120 


0,537 


0.783 


0,950 


2.024 


-0.424 


-0.770 


o./8y 


0.178 




SKEy HEAS 


-o.m 


-0.580 


-0,445 


0.193 


0.537 


0.781 


0,950 


2.024 


-0.621 


H).7?5 


0.740 


0.074 





STD ERR 



0.387 0,251 0.188 1.098 O.OOO 0.009 0.000 0.000 0.000 0.125 0.430 0.400 
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SOUTH BRANCH HIJSKOKA RIVER - HISTORIC SEQUENTIAL RUTOFF EQUATIONS 
VERIFICATION RUNS - ACTUAL LAKE LEVELS USED AT KAyAGAHA AND BAYS 
SIMIATION PERIOD 1973 - 1930 



ITERATION 3 : MONTH-END LAKE LEVEL («) SUW1ARY FOR BASIN - BURNT ! NODE NUMBER 1 



YEAR TYPE JAN FEB HAR APR MAY JUN JUL AUG SEP OCT NOV DEC 



1975 CALC 


30.099 


30.314 


30,551 


30.780 


30.780 


30.780 


30.780 


30.780 


29.830 


29.387 


29.609 


29.870 


1976 CALC 


30.099 


30.314 


30.551 


30.780 


30.780 


30.780 


30,780 


30,780 


29.830 


29.387 


29.524 


29.729 


1977 CALC 


29.884 


30.001 


30.551 


30.780 


30.780 


30.775 


30.780 


30.780 


29.830 


29.387 


29.624 


29.870 


1978 CALC 


30.099 


30.236 


30,551 


30.780 


30.780 


30.780 


30.730 


30.780 


29.830 


29.387 


29.624 


29.870 


1979 CALC 


30.099 


30.285 


30.5S1 


30.780 


30.780 


30.780 


30.780 


30.780 


29.830 


29.387 


29.621 


29.870 


1980 CALC 


30.099 


30.289 


30.^1 


30.780 


30.780 


30.780 


30.780 


30.780 


2?.830 


2?. 387 


29.624 


29.870 



MEAN CALC 30,064 30.240 30.551 30.780 30.780 30.779 30.780 30.790 29.830 29.387 29.605 29.846 
SDEV CALC 0.088 0,120 0.000 0,000 0,000 0.002 0.000 0.000 0.000 0.000 0.040 0.058 
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SOUTH BRANCH MUSKOKA RIVER - HISTORIC SEQUENTIAL RUNOFF EQUATIONS 
VERIFICATION RUNS - ACTUAL LAKE LEVELS USED AT KAUAGAtIA AND FAYS 
SKWLATION PERIOB 1975 - 1980 



ITERATION 3 ! W»aH-END LAKE LEVa (tt) SUilMARY FOR BASIN - JOE ! NODE NUHKR 2 



YEAR TYPE JAN FEB HAR APR MAY JUN JUL AUG SEP OCT NOV DEC 



1975 CALC 


421.550 


421.295 


421.000 


421.900 


421.850 


421.850 


421.850 


421.850 


421.700 


421.550 


421.550 


421.550 




1976 CALC 


421.550 


421.295 


421.000 


421.900 


421.850 


421.850 


421.850 


421.350 


421.700 


121.550 


421.550 


421.550 




1977 CALC 


421.550 


421.295 


421.000 


421.900 


421. S50 


421. 8S0 


421.850 


421.850 


421.700 


421.550 


421.550 


42! .550 




1978 CALC 


421.550 


421.295 


421.000 


421.900 


421.850 


421.850 


421.850 


121.850 


121.700 


421.550 


421.550 


421.550 




1979 CALC 


421.550 


421.295 


421.000 


421.900 


421.850 


421.850 


421.850 


421.850 


421.700 


421.550 


421.550 


421.550 




1980 CALC 


421.550 


421.295 


421.000 


421.900 


421.850 


421.850 


421.850 


421.850 


421.700 


421.550 


421.550 


421.550 





MEAN CALC 421.550 421.295 421.000 421.900 421.850 421.850 421.850 421.850 421.700 421.550 421.550 421.550 
SDEV CALC 0.000 0.000 0.000 0.000 0.000 0.000 0,000 0.000 0,000 0.000 0.000 0.000 
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SOUTH BR/tfO MUSKOKA RIVER - HISTORIC SEQUENTIAL RUNOFF EQUATIONS 
VERIFICATION RUNS - ACTUAL LAKE LEVELS USED AT KA«AGA«A AMD BAYS 
SIMULATION PERIOD 1?75 - 1980 



ITERATION 3 ! MONTH-END LAKE LEVEL (fi) STOARY FOR BASIN - KAWA ) NODE NUMBER 3 



YEAR TYPE JAN FEB MAR APR MAY JUH JUL AUG SEP OCT MOV DEC 



1975 CALC 


5.050 


4.730 


4.440 


5.090 


5.310 


5.370 


5.310 


5.240 


5.080 


4.410 


4.740 


5.080 


197i CALC 


i.710 


4.240 


4.800 


5.430 


5.390 


5.430 


3.240 


5.080 


4.940 


4.340 


4.280 


4.470 


1977 CALC 


4.100 


3.840 


4.450 


5.234 


5.350 


5.270 


5.210 


5.210 


4.940 


4.430 


5.120 


5.350 


1978 CALC 


4.730 


4.220 


3.720 


4.273 


5.490 


5.410 


5.530 


5.410 


5.040 


4.480 


4.890 


5.090 


1979 CALC 


3.060 


4.530 


4.430 


5.590 


5.410 


5.580 


5.490 


5.440 


4.920 


4.890 


4.920 


4.790 


1980 CALC 


4,730 


4.270 


4.770 


5.490 


5.470 


3.720 


5.720 


5.310 


5,120 


4.840 


5.020 


5.050 



MEAN CALC 4,730 4.308 4,443 5.185 5.170 5.497 5,420 5.282 5.010 4.442 4.828 4.972 
SDEV CALC 0.349 0.297 0.397 0.482 0.122 0.148 0,194 0.134 0.084 0,221 0.297 0.303 
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SOUTH HIANCH NUSKOKA RIVER - HISTORIC SEQUENTIAL RUNOFF EQUATIOWS 
VERIFICATION RUNS - ACTUAL LAKE LEVELS USED AT KAWA6AHA AND BAYS 
SlttULATION PERIOD 1973 ■ 1980 



ITERATION 3 ! WNTH-END LAKE LEVa («) SUWWRY FOR BASIN - RA66ED 5 NODE WMBER 4 



YEAR TYPE JAN FEB MAR APR HAY M JUL AUG SEP OCT NOV DEC 



1975 CALC 


431.982 


432.143 


432.320 


432.492 


432. A70 


432.470 


432.470 


132.040 


431.500 


431.405 


431,706 


431.810 


1976 CALC 


431.982 


432.143 


432.320 


432.492 


432.470 


432.470 


432.470 


432.0'tO 


431.500 


431.355 


431.704 


431.810 


1977 CALC 


431.982 


432.143 


432.320 


432.492 


432.670 


432.470 


432.470 


432.040 


431.500 


431.405 


431.706 


431.810 


1978 CALC 


431.982 


432.143 


432.320 


432.492 


432. i70 


432.470 


432.470 


432.040 


431.500 


431.405 


431.704 


431.310 


1979 CALC 


431.982 


432.143 


432.320 


432.492 


432.670 


432.470 


432.470 


432.040 


431.500 


431.405 


431.706 


431.810 


1980 CALC 


431.992 


432.143 


432.320 


432.492 


432.470 


432.470 


432.470 


432.040 


131.500 


431,405 


431.704 


431.810 



MEAN CALC 431.982 432.143 432.320 432.492 432.470 432.470 432.470 432.040 431.500 431,594 431.704 431.810 
SDEV CALC 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0,000 0.000 0.020 0.000 O.OOO 



I 
I 
I 
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SOUTH BRANCH NUSKOKA RIVER - HISTORIC SEQUENTIAL RUNOFF EQUATIONS 
VERIFICATION RUNS - ACTUAL LAKE LEVELS USED AT KAyAGAItt AKD BAYS 
SmULATION PERIOD 1975 - 1980 



ITERATIW* 3 ; MONTH-EWI LAKE LEVa (H) SU««ARY FOR BASIN - TSC i NODE iftJHBER 5 



YEAR TYPE JAN FEB MAR APR HAY M JUL AUG SEP OCT NOV DEC 



1975 CALC 417.513 417.355 417.210 417. S80 417.820 417.820 417.820 417.820 417.784 417.717 417.450 417.580 

1976 CALC 417.513 417.355 417.210 417.880 417.820 417.820 117.820 417.320 417.784 417.717 417.450 417.580 

1977 CALC 417.513 417.355 417.210 417.880 417.820 417.820 417.820 417.820 417,764 417.717 417.450 417. 5B0 

1978 CALC 417.313 417.355 417.210 417,^0 417.820 417.820 417.820 417.820 417.784 117.717 417.450 417.580 

1979 CALC 417.513 417.355 417.210 417.880 417.820 417.820 417.820 417.820 417.784 417.717 417.450 417.580 

1980 CALC 417.513 417,355 117,210 417,880 417.820 417,820 417.820 417.820 417.784 417.717 417,450 417.580 



MEAN CALC 417.513 417.355 417.210 417.880 417.820 417.820 417,820 417.820 417.784 417.717 417.450 417.580 
SDEV CALC O.OOO 0.000 0,000 0,000 0,000 0,000 0,000 0.000 0,000 0,000 0.000 0.000 



F-36 



SOUTH mmn huskoka river - historic sequential rumoff equations 

VERIFICATIOH RUNS - ACTUAL LAKE LEVELS USED AT KAyAGAMA AND BAYS 
SIMULATION PERIOD 1973 - 1980 



ITERATION 3 ! MONTH-END LAKE LEVEL («> SUKMARY FOR 8ASIN - BAY ', NODE NtMER & 



YEAR TYPE JAM FEB HAft APR HAY JUN JUL AU5 SEP OCT NOV DEC 



1975 CALC 


5.070 


4.990 


4,400 


4.880 


5.130 


5.180 


5.220 


5.220 


5.210 


5.140 


3.210 


5.160 


1976 CALC 


5.030 


4.320 


4.810 


5.180 


5.140 


5.240 


5.170 


S.UO 


5.210 


5.1W 


5,180 


5.080 


1977 CALC 


4,920 


4.810 


4.&&0 


5.250 


5.230 


5.220 


5.190 


5.250 


5.270 


5.140 


5.230 


5.130 


1978 CALC 


4.890 


4.790 


4.400 


5.210 


5.190 


5.160 


5,180 


5.270 


5.210 


5.140 


5.240 


5.180 


1979 CALC 


5,U0 


4.770 


4.680 


5.290 


5.240 


4.880 


5.220 


5.300 


5.140 


5.200 


5.300 


5.180 


1980 CALC 


5.050 


4.820 


4.770 


5.130 


5,190 


5.240 


5.320 


5.250 


5.300 


5.180 


5.190 


5.100 


MEAN CALC 


S.020 


4.333 


4.i20 


5.157 


5.187 


5.153 


5.217 


5.242 


5.223 


5.172 


5,225 


5.138 



SDEV CALC 0.100 0.079 0.179 0.144 0,045 0,138 0.055 0.048 0.056 0.022 0.043 0.042 



TABLE F-3 OUTPUT DATA FILE TO BE USED AS 

INPUT TO THE SERPENT RIVER rtOOEL 



F-37 



6 

1 
2 


BURNT 















2 


JOE 














5 


1 


1 












3 
6 


KAUA 















4 


RAGGED 












5 

















5 


TSC 














6 


*> 


2 4 












6 


B^Y 

















2 


5 3 












1975 


1 


0.3613 


1.4706 


7.1010 


0.9461 


4.9424 


20.3000 


1973 


2 


0.0401 


1.4467 


3.6937 


0.5300 


4,5857 


23.1000 


1973 


3 


0.8106 


3.0337 


8.5419 


1.2406 


7.6528 


40.1000 


197-5 


4 


0.9364 


0.3677 


5.9051 


0.9429 


5.0299 


35.1061 


1975 


5 


1.1818 


2,6978 


10.9373 


1.3515 


9.3452 


31.5000 


1973 


6 


0,1162 


0,4006 


2.0831 


0.5654 


1.3332 


4.6200 


1975 


7 


0.2094 


0.5109 


2.2014 


0.4605 


1.5487 


3.3300 


1975 


3 


0.0850 


0.2130 


1,3973 


1.7090 


2.2363 


5.0600 


1973 


9 


3. 3456 


4.4661 


3.5487 


1.6198 


6.9413 


13.8000 


1975 


10 


2,00«5 


2.7419 


10.2693 


0.2043 


4.2939 


20.1000 


1975 


11 


0.0000 


0.7311 


0.4042 


0.6080 


3.4638 


10.4567 


1973 


12 


0.7117 


2.3303 


5.9997 


1.7365 


3.3395 


32,0000 


197i 


1 


0.2357 


1.2440 


9.4943 


0.S257 


4.4150 


25.0000 


1976 


-> 


0.1133 


1.5931 


10.3982 


0,6331 


4.7244 


27,5932 


1976 


3 


1.S707 


5.2085 


8.3061 


2.7128 


15.0740 


50,5000 


1976 


4 


i.3in 


2.7773 


13.4993 


1.3405 


11 .11 14 


69.1000 


1976 


5 


0.7048 


1.6708 


6.5063 


0.9831 


5.7150 


24.0000 


1976 


6 


0.2617 


0,6554 


5,3330 


0.6759 


2,4323 


8.1900 


1976 


7 


0.9376 


2.0412 


7.4027 


1.6107 


5.7657 


13.6000 


1976 


3 


0.0423 


0,1632 


3.2644 


1.7325 


2,2710 


6.3100 


1976 


9 


3.6340 


3.9733 


1.7664 


1,3232 


5.4121 


6.1800 


1976 


10 


1.71 '25 


2.1524 


7.5524 


0,0000 


2.7627 


12.0197 


1976 


U 


0.0000 


0.3614 


3.0741 


0.0398 


1 . 5589 


9.0099 


1976 


12 


0.0000 


0.6670 


1.6614 


0.3723 


3,0339 


13,6996 


1977 


1 


0.0000 


0.4099 


6.5211 


0.0634 


1.6313 


16.1651 


1977 


•y 


0.0000 


1.0299 


5,1669 


0,0338 


2.7333 


11.104S 


1977 


3 


0.9045 


4. 4111 


4.6263 


2.6949 


13.3507 


41.5000 


1977 


4 


1.1013 


1.2321 


5.7233 


1,0052 


6.2112 


38.1135 


1977 


5 


0.3238 


0.9276 


3.3417 


0.4324 


3.3279 


ii.9000 


1977 


6 


0.0000 


0,1133 


2.4453 


0.2973 


0.3299 


1.0333 


1977 


7 


0.0905 


0.3031 


1.5977 


0.3459 


1.0095 


3.04 00 


1977 


3 


0.1392 


0.2942 


1.1410 


1,7206 


2.4163 


3.2900 


1977 





3.3602 


4.4777 


4.9089 


1,6203 


6,9404 


14.4000 


1977 


10 


2.7127 


4.1435 


10.1558 


1.0136 


3.3199 


33.7000 


1977 


11 


0.6782 


2.1163 


3.2243 


1.4'?35 


7.5134 


24.1000 


1977 


12 


0.5437 


1,9512 


5.4447 


1,5133 


6.9369 


27.5000 



I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 



F-38 



TABLE F-8 OUTPUT DftlA FILE TO BE USED AS 
COHT'D INPUT TO THE SERPENT RIVER MODEL 



1978 


1 


0.1808 


1.1959 


12.2544 


0.351? 


4,3120 


29,6036 


1973 




0.0000 


1.1390 


9.4736 


0.301? 


3.4173 


19.4802 


1975 


3 


0.1575 


2.0286 


3.3607 


0.7998 


5.1762 


27,7000 


1973 


4 


0.5333 


0.0755 


4.7228 


0.2063 


2.7312 


22.9701 


1978 


5 


1.9564 


4.3982 


3.1503 


3,3595 


13.5425 


53.4000 


1978 


6 


0.3311 


0.8926 


2,9232 


0.3726 


3.2030 


11.3000 


1978 


7 


0.1349 


0.4646 


2,3301 


0,4307 


1.4426 


3.9700 


1978 


3 


0.2757 


0.5327 


3.1949 


1.3421 


3.0120 


6.3200 


1978 


9 


3.9'U3 


4.6561 


6.5095 


1.7338 


7.4934 


19.7000 


1973 


10 


2.3760 


3.4931 


9.3260 


0,6593 


6.5120 


25.6000 


1978 


11 


0.0643 


0.9571 


2.9359 


0,8062 


4 . 3228 


14,7000 


1973 


12 


0.2332 


1.4555 


4.3134 


1.2343 


5.5722 


21.4000 


1979 


1 


0.4782 


1.6308 


6.2437 


0.9813 


5.2665 


21 .5000 


1979 


1 


0.0000 


1.34 40 


10.4536 


0.'t364 


4.0409 


30.4393 


1979 


3 


1.9352 


5.4263 


11.6271 


2.7313 


14,7070 


50.4762 


1979 


^1 


2,0513 


3.2257 


10.4331 


2.4685 


12.2674 


57.4000 


1979 


5 


1.4071 


3.0641 


14.0583 


1.9978 


10.1573 


44.7000 


1979 


& 


0.0897 


0.3410 


2.7904 


0.4741 


1.5956 


15.9000 


1979 


7 


0.6841 


1.4954 


5.4942 


1.1772 


4.3467 


5.0900 


1979 


8 


0.3797 


0.8244 


3.9346 


2.1366 


4.1753 


9.5200 


1979 


9 


3.8232 


4.4514 


7.9438 


1.6176 


6,9380 


22.1000 


1979 


10 


2.8958 


4.4353 


7.6737 


1.1370 


3.9592 


23.3000 


1979 


11 


1.1801 


3.2534 


12.6345 


2.2183 


11.0786 


4 3.1000 


1979 


12 


0.7713 


2.5072 


12,0633 


1.93 40 


3.3369 


40.1000 


1980 


1 


0.7574 


2.4024 


8.5570 


1 .6428 


7.4894 


30.1299 


1930 


2 


0.0000 


1.3564 


9.4749 


0,5031 


4,0454 


24.7317 


1980 


3 


0.5754 


2.6855 


1.3657 


1.2357 


7,3110 


23.0000 


1930 


A 


2.3384 


3.8482 


15.4622 


2.5095 


14.0010 


73.5000 


1980 


5 


0.6305 


1.5015 


7.3841 


0.3103 


5.0421 


20.6000 


1<'80 


6 


1.0439 


2.1907 


5.3357 


1.7374 


7.0254 


22.5000 


1980 


7 


0.9368 


2.0876 


6.5221 


1 .5648 


6.0121 


17.0000 


1930 


3 


0.6119 


1.2729 


9.6313 


2.4^30 


5.4097 


22.3000 


1980 


9 


4.3357 


5.4781 


7.0995 


2.2137 


9.90<?9 


24.3000 


1930 


10 


2.7832 


4.2473 


9.9912 


1,0529 


3.5035 


33.9000 


1980 


11 


0.5712 


2.0721 


7.2S34 


1.5361 


7.7036 


30. 1244 


1980 


12 


o.3no 


1.4438 


6.0901 


1,1064 


5,3733 


24.3000 
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TABLE F-9 ajMHARY OF Ct^LC'JLATa DISCHARGE FRGM 3ASIN 

TO BE USED AS INPUT FILE FOR DOUNSTREAil BASINS 



o 6 BAY 3ASIN DISCHARGE 

1973 20.300 23.100 40.100 33.406 31,500 4.620 3.380 5.060 13.300 20.100 10.457 32.000 

1976 25.000 27.593 50.500 49.100 24.000 8.490 13.600 6.310 6.130 12.050 9.010 13.700 

1977 16.165 14.105 41.500 38.444 13.900 4. 056 3.040 3.290 14.400 33.700 24.100 27.500 
1973 29,609 19,480 27,700 22.920 53.400 11,300 3.970 6,320 19,700 25.600 14.700 21.400 

1979 21.500 30.459 50.476 57.400 44.700 15.300 5.090 9.520 22.100 28.300 43.100 40.100 

1980 30,130 24,782 23.000 79.500 20,600 22.500 17,000 22,800 24,300 33,900 30,124 24,300 





DATE DUE 
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